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Notations and conventions 



All cultural products contain a mixture 
of two elements: conventions and inven- 
tions. 

J. G. Cawelti 



The conventions will mainly follow the book of Peskin and Schroeder Q147Q . We use the metric 
tensor 



with Greek indices running over 0,1,2,3- Repeated indices are summed in all cases. Light italic 
roman type will be used for four- vectors, while boldface italic will be used for three- vectors. 
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Light-cone vectors 

Light-cone vectors will be indicated as 

a^ = [a^, Oj^] = 
The dot-product in light-cone components is 



a a 



a-b = a^b +b a'^ — Uj-bj 
= a+b^ +b^a^ -a'b' 
= a^b^ + b^a^ — a^bx — ayby 

The two-dimensional transverse parts of the vectors will be written in boldface with an index T 
and Latin indices will be used to denote the two transverse components only. Note that 



[0, 0, a 
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Notations and conventions 



We introduce the projector on the transverse subspace 

/O o\ 
0-100 
0-10 

yo oy 



§1 



We define the antisymmetric tensor so that 



-0123 



-1. 



and we define the transverse part of the antisymmetric tensor as 



03mv 



Dirac matrices 



Dirac matrices will be often expressed in the chiral or Weyl representation, i.e 
and we will make use of the Dirac structure 



1 
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Introduction 

The most difficult part of a trip is to 
cross the doorway. 

P. Terentius Varro 



In this thesis I will discuss three different ways to observe the transverse spin of quarks inside 
the nucleons. Before embarking on such an undertaking, I would like to spend a few pages on 
explaining what makes this problem so interesting to justify investing years of research on it. This 
introduction is meant especially for nonexperts, since I will review notions well known to the 
experts in the field. 



1.1 The structure of matter 

When we talk about quarks inside nucleons we are referring to the best paradigm we currently have 
to describe the elementary structure of matter. The comprehension of this elementary structure is 
a question that has allured philosophers and scientists since the historical origins of philosophical 
thought. It is striking to observe that as early as six hundred years BC, Greek philosophers were 
already wondering: are there fundamental elements in nature, what are they and how do they 
interact? Today, after more than two millennia, we learned a lot about the structure of matter, but 
some of the most important questions still elude our comprehension. We are still engaged in one 
of the oldest quests of human mind. 



Since 1803, when Dalton suggested his atomic hypothesis [85], we have gradually realized 



that almost all matter on earth is made up of atoms. Atoms contain electrons - identified for 
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the first time by J. J. Thomson in 1897 [ |166| , |167| ] - and nuclei - introduced for the first time 
by E. Rutherford in 1911 D153| , pS]]. The efforts to explain precisely the structure of atoms and 
the electromagnetic interaction binding together electrons and nuclei lead to two of the major 
achievements of physics in the last century: Quantum Mechanics and Quantum Electrodynamics 
(QED). 

In the meantime, more investigations were carried out to grasp the structure of the nucleus 
inside atoms. The smallest known nucleus was identified with a single particle [ ]154| ], the proton, 
while a second constituent of heavier nuclei, the neutron, was eventually observed by J. Chad- 
wick [[751] ■ Since they are the constituents of the nucleus, protons and neutrons are referred to as 
nucleons. They are kept together by the nuclear force, of which at the moment we have only an 
incomplete understanding. 

Although the electrons are responsible for the chemical properties of atoms, they account for a 
very small fraction of the mass of the atom. The mass of an electron is about 0.5 1 1 MeV, while the 
mass of a proton is about 938 MeV. Therefore, nucleons make up for more than 99.9% of ordinary 
atomic matter. If we want to understand matter, we cannot set aside the problem of explaining the 
structure of nucleons. Nucleons belong to the more general class of hadrons, of which they are 
the most abundant specimen. At first, hadrons were classified as elementary particles, i.e. without 
any internal substructure. Very soon this appeared to be an unsatisfactory hypothesis, in particular 
since there are so many of them (several dozens). Nowadays, the study of the structure of hadrons 
represents a field of research on its own, often designated with the name of hadronic physics. An 
up-to-date review of the field can be found in Refs. ffS] and [M 



1.2 Hadrons and deep inelastic scattering 



To interpret the information available on the properties of hadrons in 1964, M. Gell-Mann [ |99D 
and G. Zweig [ ]172| ] independently suggested that hadrons are composed of smaller constituents, 
the quarks, having spin 1/2, a fractional electric charge and a new degree of freedom, called 
flavor. This model is often referred to as constituent quark model. Gell-Mann himself seemed 
not to believe in the existence of quarks as real entities, but rather regarded them as convenient 



concepts 0990. One of the reasons to be skeptical about the real existence of quarks was that they 
have a charge that is just a fraction of the electron charge, while the charge of all other elementary 
particles is an integer multiple of that. 

The quark model aimed at describing the mass, charge and spin of the hadrons. For instance, 
the proton has a mass of about 1 GeV, a charge +e (the same as the electron, but with opposite 
sign) and spin equal to 1/2. According to the model, a proton with its spin, for instance, in the up 
direction is made of two quarks with flavor up and charge 2e/3 plus one quark with flavor down 
and charge — e/3. Two of the quarks have spin 1/2 in the up direction and one has spin 1/2 in the 
down direction. Each of the three quarks carries about one-third of the mass of the proton. 

Looking at the "extrinsic" properties of hadrons - like their mass, charge and spin - was not 
enough to unravel the details of their structure. To glance at the inside of hadrons, physicists 
resorted to deep inelastic scattering (DIS) experiments, as in the pioneering experiments led by 



1.3 Quantum Chromodynamics and confinement 



3 



Friedman, Kendall and Taylor at the Stanford Linear Accelerator Center (SLAC) In scat 



tering experiments, a focused beam of particles is dispersed by the interaction with a target. The 
way this dispersion takes place yields information on the structure of the target. For instance, the 
existence of the nucleus was suggested by Rutherford as an explanation to the scattering experi- 
ments of Geiger and Marsden [ p8| ] . 

Basically, particle accelerators as the one at SLAC are exploited as microscopes of extremely 
high resolution. The experiments at SLAC scattered electrons off hydrogen. The interaction pro- 
ceeds via the exchange of a virtual photon with high energy and momentum. A measure of the 
resolution of the experiment is given by the four-momentum squared of the virtual photon, Q^, or 
rather by the associated wavelength h/Q. The SLAC experiments reached a maximum of 7.4 
GeV^, corresponding to a resolution of the order of 1/10 of the proton size. 

The results of the SLAC experiments indicated that the scattering data did not exhibit a (strong) 
dependence on 2^. They depended rather on the variable that was later to be named x-Bjorken, Xg, 
in honor of J. Bjorken.Q This property, called scaling, was predicted by Bjorken himself p6[ ] and 
explained by R. Feynman [^^], who introduced the parton model: the proton was pictured to 
be a collection of almost /ree point-like constituents off which the electrons scatter incoherently. 
The constituents of the proton were initially called partons, but it soon became clear that they had 
a lot in common with the quarks of Gell-Mann and Zweig. 

Feynman's partons have spin 1/2, fractional electric charge and flavor, but they have a very 
small mass compared to Gell-Mann's quarks, a few MeV against about 300 MeV. Consequently, 
we call them current quarks, to distinguish them from the constituent quarks of the quark model. 
But this is not the only difference between the two models. In the constituent quark model, the 
proton is made up just of three quarks, while in the parton model it turns out that there is a huge 
number of quark- antiquark pairs, together with a huge number of electrically neutral particles, later 
to be identified as gluons. 

Deep inelastic scattering experiments are performed in some of the world's largest experimen- 
tal facilities for high energy physics, such as CERN, SLAC, DESY, BNL. They usually employ 
beams of electrons, positrons, muons or, more rarely, neutrinos. They scatter off different kinds 
of fixed targets or off a beam of protons, and they operate at different energies and kinematic 
coverage. 



1.3 Quantum Chromodynamics and confinement 

The parton model raised a profound question. We experience that matter, at least in the normal 
conditions on earth, is composed of hadrons - it is the so-called hadronic matter. If quarks are 
the hypothetical constituents of hadrons, they must be bound extremely tight to explain why we 
have never directly observed a single isolated quark, nor a different state of matter other than 
hadronic. This essential feature of quark dynamics is known as confinement. Yet, deep inelastic 
scattering suggests that in the interaction with a high- 2^ virtual photon, quarks behave as if they 

'We will properly define in Chap. ^ 
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Figure 1.1. The running of the electromagnetic coupling constant. The gray area roughly indicates where 
perturbation theory is not trustworthy anymore. Note the huge scales necessary to achieve a appreciable 
increase of the coupling. 



were essentially free. This property is known as asymptotic freedom. The question is then: how is it 
possible to devise a theory to reconcile these two opposing properties, confinement and asymptotic 
freedom? 

A first attempt to implement confinement was done by postulating that quarks have a color 
charge and that all detectable objects have to be colorless (cf. Refs. |142| , |143| and references therein). 
By virtue of this assumption, it is impossible to see an isolated, colored quark. On the other hand, 
such a point of view is not suited to describe asymptotic freedom. There was the need of a theory 
that could describe the binding of colored quarks as a dynamical mechanism. 

To shape a new theory of color interactions, it seemed natural to follow the example of Quan- 
tum Electrodynamics, the quantum field theory of electromagnetic interactions. QED is in essence 
a perturbative theory, which works because electromagnetic interactions are weak. In fact, elec- 
trons can be separated quite easily from atoms and observed as free particles. A measure of the 
strength of the electromagnetic interaction is given by the value of the electromagnetic coupling 
constant a ~ 1/137. From renormalization of QED, it is known that in reality a is not constant, 
but it has to increase as the momentum exchange of the interaction increases, or equivalently as 
the interaction takes place over shorter distances. However, the increase of the coupling is so 
weak (e.g. a ~ 1/135 at = 1000 GeV^) that perturbative QED works brilliantly for any elec- 
trodynamics experiment we might do. Fig. |LT] shows approximately the way the electromagnetic 
coupling constant changes with Q^. 

A breakthrough in the comprehension of quark interactions came in 1973, when D. Gross and 
F. Wilczek Q102| , |10l| ] and D. Politzer Q1490 showed that non-Abelian quantum field theories can 
display the crucial properties of asymptotic freedom, i.e. the interaction they describe is weak 
at high momentum transfer (or long distances). This discovery prompted the birth of Quantum 
Chromodynamics (QCD), a non-Abelian field theory of color interactions. 

The difference between QED and QCD can be likened to the difference between the attraction 
forces of two opposite magnetic poles and of two ends of a spring. In the first case, we know that 
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Figure 1.2. The running of the strong coupling constant. The gray area roughly indicates where perturbation 
theory is not trustworthy anymore. 



increasing the distance between the magnets, the attraction diminishes, while, in the second case, 
separating the two ends the force will increase more and more. 

In field theory language, the electromagnetic coupling constant is reduced at large distances due 
to the effect of vacuum polarization, which is responsible for a screening of the bare electric charge. 
On the contrary, the color coupling constant is reduced at short distances because the vacuum 
polarization induces an antiscreening of the charge, or equivalently an enhancement of the charge 
at large distances. The reason for this different behavior is that gluons, the mediator of the color 
interaction, carry color charge themselves, while photons, the mediator of the electromagnetic 
interaction, are chargeless. Ref. |147| (p. 541) and Ref. ^ (p. 5) present enlightening discussions 
on antiscreening and asymptotic freedom. 

QCD is able to justify asymptotic freedom, but what about confinement? At the moment, we 
know that QCD is not in contradiction with confinement and might in fact explain it, but we are 
not able to demonstrate this statement. As in the case of QED, the strength of color interaction is 
measured by the strong coupling constant Us, which has a value of about 0.1 17 at = 8.3 x 10^ 
GeV^. But at lower energy scales, e.g. 2^ ~ 1 GeV^, the coupling constant grows and becomes 
of the order of 1. The running of the strong coupling constant is illustrated in Fig. [0[ We might 
deduce that the increase of the coupling constant is a sign of the onset of confinement. In reality, 
we can only conclude that at low energies we enter a regime where perturbation theory cannot be 
trusted. Therefore, even if QCD is in principle a consistent theory at any energy scale, we cannot 
use standard techniques to draw conclusions about its behavior in the nonperturbative regime. To a 
certain extent, we cannot be sure that QCD is the correct theory in this regime: maybe it is simply 
an asymptotic approximation of a more profound theory. 

In practice, we have to make a distinction between two major branches of QCD: perturbative 
and nonperturbative, or short-distance and long-distance. Perturbative QCD is relatively well un- 
derstood. It is essentially similar to QED, it is based on Feynman-diagram approach, although it 
often requires larger sets of diagrams to attain the desired accuracy. The theory contains pointlike 
and almost massless fermions (the so-called current quarks) and massless bosons to carry their 
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Figure 1.3. One of the greatest achievements of perturbative QCD: world data on the Q dependence of the 



unpolarized structure function compared to theoretical expectations (solid curves) [138]. 



interactions (the gluons). Probably one of the most important achievements of perturbative QCD 
is the study of the way deep inelastic scattering data change with Q^. The striking agreement be- 
tween theory and experiments is shown in Fig. [TT3|. On the other hand, nonperturbative QCD is 
poorly understood and it is a challenging playground for fundamental physics. At the moment, 
our understanding of this theory relies on lattice calculations, effective chiral field theories, and 
phenomenological models. 

We know that nonperturbative QCD should display confinement as a fundamental property, 
at least under normal conditions. Lattice calculations already provide strong evidence that the 
quark-quark interaction potential increases linearly, and is therefore a confining potential [^. 
It is essential to understand from first principles why this occurs, and it is desirable to explore 
how it is possible to achieve a deconfined phase of QCD, maybe under extreme conditions (e.g. 
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neutron stars). We also know that in QCD chiral symmetry is approximately valid. On the other 
hand, the existence of pions, which are nearly massless, suggests that chiral symmetry should be 
spontaneously broken, with pions being Goldstone bosons. Nonperturbative QCD should be able 
to explain this feature. Nonperturbative QCD should also explain the transition between massless 
current quarks and constituent quarks. 

Finally and more generally, nonperturbative QCD should lead to a reliable quantitative descrip- 
tion of the structure of hadrons and of hadronic phenomena. The question at the heart of hadronic 
physics is: what is the structure of hadrons in terms of their quark and gluon constituents? There- 
fore, we might define hadronic physics as the branch of physics that deals with understanding 
QCD, and in particular nonperturbative QCD. 



1.4 Spin physics and the transversity distribution 

One of the key questions in understanding the structure of hadrons is: where does the spin of the 
nucleons come from? In the constituent quark model, the spin of the quarks adds up to yield the 
total spin of the proton. Deep inelastic scattering experiments, however, show the importance of 
other contributions, such as the spin of the gluons and the orbital angular momenta of quarks and 
gluons. A measure of the quark spin contribution is given by the distribution function g^, often 
denoted as /^q and usually called the helicity distribution. In a frame of reference where the hadron 
is moving with a very large speed {infinite momentum frame) and if the direction of its spin is 
longitudinal to its motion, the helicity distribution describes the number of quarks with their spin 
aligned with that of the hadron minus the number of quarks with opposite spin, it is therefore a 
measure of the longitudinal spin of the quarks in the hadron. The quark helicity distribution has 



been measured with a good precision, as shown in Fig. |L4|on page Naively, if the spin of the 
hadron is entirely due to the quark spin as in the constituent quark model, we expect to have a 
net balance of one quark spinning in the direction of the proton and thus accounting for the whole 
proton spin. In reality, it turns out that (the integral of) the helicity distribution accounts for only 
about 30% of the proton spin! We expect thus that the missing spin is provided by the gluon 
spin and by the orbital angular momentum of quarks and gluons. These two quantities have not 
been measured yet. Even worse, we don't know if it is possible to measure the orbital angular 



momentum directly [ |123| , |103| ,p3l]. Tab. jlTTJon the following page shows a list of all polarized 



deep inelastic scattering experiments, together with their typical energies and the characteristics 
of their beams and targets. The kinematic coverage of each experiment is indicated in the table 
by its average (GeV^) and Xg range. Columns and Pj give the average or typical beam and 
target polarizations as quoted by each experimental group. The column labeled ^ is an estimate 
of the total nucleon luminosity (# of nucleons/cm^ times # of beam particles/s) in units of 10^^ 
nucleons/cm^/s for each experiment. 

So far we talked about the longitudinal spin of the quarks inside the proton, but what about the 
transverse spin? The observable we have to take into consideration is the transversity distribution. 
In the infinite momentum frame with the proton spin transverse to the direction of motion, the 
transversity distribution describes the number of quarks with their spin aligned with that of the 
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Table 1.1. Summary of polarized deep inelastic measurements 



Lab 


Exp. 


Year 


Beam 


\W; / 






Target 




if 


SLAC 


E80 


75 


10-16 GeV e 


2 


0.1 - 0.5 


85% 


T-T-biitannl 

X X 1. ' \A I, iX J. 1. V/ X 


50% 


400 




E130 


80 


16-23 GeV 


5 


0.1 - 0.6 


81% 


l-T-hntannl 

X X I / IX L iX J. 1. V/ X 


58% 


400 




E142 


92 


19-26 GeV 


2 


03 - 6 


39% 


^He 


35% 


2000 






yj 


lU VJC V c 




0^ OS 

— yJ.O 


85^!^ 

OJ /V 




1\J/V 




















25% 


1000 




E154 


95 


48 GeV 


5 


0.01 -0.7 


82% 


^He 


38% 


3000 




E155 


97 


48 GeV 


5 


0.01 -0.9 


81% 




90% 


1000 
















LiD 


22% 


1000 






99 


30 GeV p 




02 - 9 

\j .\j jl. \j,y 


83% 


NH3 


75% 


1000 

J- \J\J\J 
















LiD 


11% 


1000 


CERN 


EMC 


85 


100-200 GeV u+ 


11 


0.01 -0.7 


79% 




78% 


0.3 




SMC 


92 


100 GeV u+ 


4.6 


006 - 6 


82% 


D-hiitatiol 

i ' I / IX L tX J. 1. V/ X 


35% 


0.3 






93 


190 GeV u+ 


10 


003 - 7 


80% 


H-butanol 


86% 


0.6 






94-95 








81 % 


l_y ULiUClllUi 


50% 


6 






96 








77% 


NH3 


89% 


0.6 


DESY 


HFRMFS 


95 


28 GeV c+ 


2 5 


02 - 6 


55% 




46% 


1 






96-97 








55% 


H 


88% 


0.1 






98 


28 GeV 






55% 


D 


85% 


0.2 






99-00 


28 GeV e+ 






55% 


D 


85% 


0.2 






01-? 


28 GeV 






55% 


H 


85%* 


0.2 


CERN 


COMPASS 


01-? 


190 GeVjU+ 


10 


0.005 - 0.6 


80% 


NH3 


90% 


3 
















LiD 


40% 


3 


BNL 


RHIC 


02-? 


200 GeV p-p 


~ 100 


0.05 - 0.6 


70% 


Collider 


70% 


2 


DESY 


ZEUS/Hl 


02-? 


28 X 800 GeV e-p 


22 


0.00006 - 0.6 


50% 


Collider 




0.2 



* Transversely polarized target 
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Figure 1.4. The helicity distribution of the proton, as a function of the fractional momentum x carried by 
quarks /^. 



hadron minus the number of quarks with opposite spin, it is therefore a measure of the transverse 
spin of the quarks in the hadron. The transversity distribution looks very similar to the helicity 
distribution, as at first sight they seem related by rotational symmetry. However, we cannot forget 
the fact that the interpretation of helicity and transversity holds true only in the infinite momentum 
frame, where the direction of motion of the hadron breaks rotational symmetry. 

In the rest frame of the proton, there is a probability to find quark spins aligned with the proton's 
spin. This probability is obviously the same no matter what the spin orientation is. If we now boost 
the proton to a large speed in the direction of its spin, the alignment probability will correspond to 
the helicity distribution. If we boost the proton in a direction transverse to its spin, the alignment 
probability will correspond to the transversity distribution. In a nonrelativistic situation, Galilean 
boosts will not affect the spin distribution and we would still expect helicity and transversity to 
be equal to each other and to the spin distribution in the rest frame. But in a relativistic context, 
Lorentz boosts can affect the spin distribution and can cause helicity and transversity to be different 
from each other and from the rest-frame distribution. The way this difference arises depends on 
the inner structure of the nucleon. 



The transverse spin of quarks is thus another missing piece in the proton spin puzzle. It can 
give new information on the dynamics of quarks inside hadrons, complementary to the helicity 
distribution. In spite of this, the transversity distribution escaped notice until 1979, when it was 
introduced by J. Ralston and D. Soper p51| ]. In the last decade, it has been evaluated in mod- 
els [ |1 16| , |120| , |148| , p0| ] and lattice computations [p4|]. At this point, an experimental measurement 
will be needed to put all these calculations on test, but unfortunately the transversity distribution is 
an elusive object to measure. Today, looking for a practical way to observe transversity is still an 
open problem. Experimental collaborations are planning its measurement at last 01O5| , ^ [73l], and 
some of them will resort to the methods discussed in this thesis. 
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1. Introduction 



1.5 Outline of the thesis 

The goal of the thesis is to discuss three different ways to observe the quark transversity distribution 
in deep inelastic scattering. 

To start with, in chapter ^ I will review the formalism of deep inelastic scattering. I will 
introduce the parton distribution functions and I will devote a particular attention to the transversity 
distribution function. I will discuss totally inclusive DIS, where only the scattered electron is 
detected, and I will show that it is not possible to measure the transversity distribution in this kind 
of process. 

In chapter^, I will turn the attention to one-particle inclusive DIS, where one of the outgoing 
hadrons is detected in coincidence with the scattered electron. I will introduce fragmentation 
functions to describe the production of hadronic fragments. In particular, I will demonstrate that 
the presence of the transverse momentum of the outgoing hadron allows the introduction of the 
Collins fragmentation function. In the cross section of one-particle inclusive DIS, I will show the 
occurrence of the product of the transversity distribution and the Collins function. Therefore, this 
suggests a first way to observe the transversity distribution of the quarks. 

In chapter ^ I will examine the more complex case of two-particle inclusive DIS, where two 
of the outgoing hadrons are detected together with the scattered electron. I will discuss how the 
presence of the relative transverse momentum between the two hadrons permits the definition of 
a new function, H^, to be connected to the transversity distribution. In the same chapter, I will 
study what happens when we assume that the two hadrons are produced only in s and p waves. In 
addition to the usual s wave contributions, I will distinguish the pure p-wave contributions and the 
sp interference contributions. This will lead to the introduction of two new fragmentation functions 
that can be observed in connection with the transversity. They are two distinct components of the 
function and they generate the second and third way I will consider to access the transversity 
distribution. 

In chapter]^, I will analyze the formalism needed to deal with spin-one hadrons in deep inelastic 
scattering. In the first part of the chapter I will focus on spin-one targets, while in the second part 
I will study the production of spin-one hadrons in the final states. This process has something in 
common with one-particle inclusive DIS (because the production of a single hadron is addressed), 
but also with two-particle inclusive DIS (because the spin-one hadron has to decay into two hadrons 
to yield information on its polarization). In particular, I will clarify the connection between spin- 
one fragmentation functions and the pure p-wave sector of the analysis of two-particle production. 

The various options described to measure the transversity distribution all involve the class of 
T-odd fragmentation functions.0 To attempt some quantitative assessments on the magnitude of T- 
odd fragmentation functions, in chapter ^ I will present a model calculation of the Collins function 
and of some of the measurable quantities in which it appears. 



^In fact, in the thesis I will not deal with the well known case of spin-half production, which involves a T-even 
fragmentation function. 



2 

Distribution functions 
and transversity 



You have never given me a transverse look. 

A. Chekhov 



In this chapter, we will introduce the concept of parton distribution functions. In order to do this, 
first of all we will review the general formalism of polarized deep inelastic scattering, starting 
from the simplest case of inclusive processes. This subject is covered in detail in books (e.g. 
Refs. IBIKT^, |133), PhD theses (e.g. Refs. |135] and |l^ and reports (e.g. Refs. ^ and gl]). 
Nevertheless, it is useful to examine the formalism from the point of view we will adopt throughout 
the remaining chapters. In the analysis of distribution functions, we will include beam and target 
polarization and partonic transverse momentum. We will limit the analysis to leading order in 1/2 
and we will only briefly mention Us corrections. 

A particular attention in this chapter and in the rest of the thesis will be reserved to the quark 
transversity distribution. The quark transversity distribution h^ Ql ISQ - also called transverse spin 
distribution [|S^] - was first introduced by Ralston and Soper [ |151[ ] and it is an essential compo- 
nent in the description of the nucleon spin. It is a chiral-odd object describing the difference of 
probabilities to find in a transversely polarized hadron a quark with spin aligned or antialigned 
to the spin of the hadron. The transversity distribution has been upstaged for many years by the 
helicity distribution, g^, which is easier to measure. However, some experimental collaborations 
are planning to measure it in the next years [ ]105| , |S3| ,[73|], possibly using one of the techniques we 
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2. Distribution functions and transversity 




P3, 



Figure 2.1. Inclusive deep inelastic scattering. 



will outline in the thesis. 



2.1 Inclusive deep inelastic scattering 

In deep inelastic scattering, an electron scatters off a nucleon via a large momentum transfer, the 
nucleon is destroyed and many hadrons are formed as a consequence of the collision. In inclusive 
events, only the scattered electron is detected, while the hadronic final states are unobserved. A 



schematic view of the process is provided by Fig. 



2.1.1 Kinematics 

In electron-nucleon scattering, an electron with momentum I scatters off a nucleon with momentum 
P, mass M and spin 5, via the exchange of a virtual photon with momentum q. The electron final 
momentum is 

We define the invariants 



{p+iy 



-{I -I 



!\2 



(2.1) 



and we introduce the variables 



IPq 



PI' 



(2.2) 



In deep inelastic scattering it is required that Q^,P ■ q ^ M^. Usually, the Bjorken limit is assumed 
(Q^,P ■ q ^ °°, Xg fixed). In particular, represents the hard scale of the process. In this thesis, 
only the leading terms in an expansion in 1 /2 will be retained. In agreement with the working 
redefinition of twist proposed by Jaffe, we will very often identify the expression "leading twist" 
with the expression "leading order in 1/2" [ |119| ]. 



2.1 Inclusive deep inelastic scattering 
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When working in the Bjorken limit, the vectors P and q can be conveniently parametrized (in 
light-cone coordinates) as 



AQ 



Q AQ 















1P+' 'J 







2x,P^ 



-x,P' 







(2.3) 



(2.4) 



This parametrization holds in any frame of reference where the virtual photon direction is antipar- 
allel to the z axis. Any frame fulfilling this requirement will be simply called collinear. The pa- 
rameter A specifies uniquely a specific collinear frame of reference. For instance, for A = Mx^/Q 
we select the nucleon rest frame, where P is purely timelike, while for A = 1 we select the so-called 
infinite momentum frame, where q is purely spacelike. 

In a 1 /2 expansion, it turns out that the plus component of P plays a dominant role. This 
statement holds regardless of the value of A, i.e. in any collinear frame. In the nucleon rest frame 
P^ is of the order of 1, while in the infinite momentum frame it is of the order of Q. However, if 
we take for instance the scalar combination P ■q,we see that the component P^q^ is of the order 
of Q^, whereas P q^ is of the order of 1, independently of the frame. Therefore, we can say that 
the plus component of the nucleon's momentum is the relevant or dominant one, although only in 
the infinite momentum frame it is truly dominant. 

We are going to define a process as soft if the relevant component of all momenta remains 
the same. In contrast, in a hard process, such as the interaction with the hard momentum q, the 
relevant component has to change. Similarly, when we describe a momentum as soft with respect 
to another, we mean that their relevant component is the same. 

In the Bjorken limit, the electron and the proton can be considered to be massless, 2P ■ I ^ s 
and = sx^y. The lepton momenta can be parametrized as 



Q A{l-y)Q 



{i-y)Q AQ 

Ay\fl ' y\fl' 



y 

1 

y 











(2.5a) 



(2.5b) 



This parametrization implies that we chose the y axis of our system as pointing in the direction of 
the vector product {—q x /'). Normally, transverse vectors and azimuthal angles will be defined as 
lying on a plane perpendicular to the direction of the virtual photon (see Fig. ^]2|on the following 
page). 



2.1.2 The hadronic tensor 

The cross section for polarized electron-nucleon scattering can be written in a general way as the 
contraction between a leptonic and a hadronic tensor 

d? G ex?' 
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Figure 2.2. Description of the vectors involved in totally inclusive deep inelastic scattering and of the 
azimuthal angle (j)^. 

where the vector S denotes the spin of the nucleon and its azimuthal angle, Xg denotes the 
helicity of the electrons and a = e^/An. Fig. |2]^ illustrates the definition of the scattering plane, 



the z axis of our coUinear frame and the azimuthal angle 

Considering the lepton to be longitudinally polarized, in the massless limit the leptonic tensor 
is given by [ |1410 



(2.7) 



= -Q^Snv + 2 {l/y + I'^lv) + 2i Xe e^ypa in 



The leptonic tensor contains all the information on the leptonic probe, which can be described 
by means of perturbative QED, while the information on the hadronic target is contained in the 
hadronic tensor 

2MW^\ci^P^S) = ^^j^^^{2nf5^'\q + P-P^)H^\P^S^P^l (2.8) 

H^''{P,S,P^) = (/',5|7^(0) |^)( jr|7^(0) \P,S). (2.9) 

The state ^ symbolizes any final state, with total momentum P^ . It is integrated over since in 
inclusive processes the final state goes undetected. By Fourier transforming the delta function and 
translating one of the current operators, we can rewrite the hadronic tensor as 

2MW^\q,P,S) = ^ j d^t, e^?-^ {P,S\J^{^)J\0) \P,S). (2.10) 



2.1 Inclusive deep inelastic scattering 



15 



In general, the structure of the hadronic tensor cannot be specified further, because this would 
require an understanding of its inner dynamics. At most, it can be parametrized in terms of struc- 
ture functions. However, the phenomenology of DIS taught us that at sufficiently high we can 
assume that the scattering of the electron takes place off a quark of mass m inside the nucleon. 
The final state 3^ can be split in a quark with momentum k plus a state X with momentum P^. 
Considering the electron-quark interaction at tree level only, the hadronic tensor can be written as 

+ V0(O) \X){X\ wm \P,S) rl {ft-m\^ 7,5) , 

where k is the momentum of the struck quark, the index q denotes the quark flavor and eq is the 
fractional charge of the quark. Note that, for simplicity, we omitted the flavor indices on the quark 
fields. The integration over the phase space of the final-state quark can be replaced by a four- 
dimensional integral with an on-shell condition, so that the hadronic tensor can be rewritten as 

2MWMV(,,P,5) = £4£ / / d'k5ik'-m')eik^-m) 

: 5(4) (p+q-k-P^) [{P,S\ WiiO) \X){X\ WjiO) \P,S) yjl (Hm^iT^j 



+ {P, S\ yfjiO) \X){X\ YtiO) \P, S) it^ - m\, yf.^ 

(2.12) 

Next, we Fourier transform the Dirac delta function and we introduce the momentum p = k — qio 



obtain 



X ((P,S| |X><X| I///0) |P,i>r| (^+i + m)uyjj 

+ (PA Vj(o) \x){x\ w,(o) \p,s) yl (V+rf -m)„r,';) ■ 

Finally, we use part of the exponential to perform a translation of the field operators and we use 
completeness to eliminate the unobserved X states, so that 

+ {P, S\ yfj ) v^,(0) |P, s)rl{^ + 4- r[^) . 
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Figure 2.3. The handbag diagram, describing the hadronic tensor at tree level. 



The hadronic tensor can be written in a more compact way by introducing the quark-quark 
correlation function 4> and the antiquark-antiquark correlation function 4> 



q J 

x^x[^{p,P,S)y^{j/, + 4 + m)f+^{p,P,S)f{j/> + i-m)y^] 

where 



(2.15) 



{my 



d'^ e-^-^(/',5|v/,(^)V/,(0)|P,5> 



X 



[P,S\ V7,(0) \X){X\ v/,(0) |P,5> 5(4) (P-P-P: 



(2.16a) 



4> 



'jApM = J^4j^'^ e-^-^(/',5| v//^) V/,(0) |/',5> 

/ j^^bi {p^sw^mx){xmoms) 5^'\p-p-p^). 



(2.16b) 



As the quark fields should carry a flavor index that we omitted, also the correlation functions are 
flavor dependent and they should be indicated more appropriately as and 4>^. For simplicity, 4> 
will be omitted henceforth. It can be accounted for simply by extending the summation over quarks 
to a summation over quarks and antiquarks. A graphical representation of the hadronic tensor at 
tree level in the parton model is given by the so-called handbag diagram, depicted in Fig. P3 . 
We parametrize the quark momentum p in the following way 



2xP+ 



(2.17) 



In our approach, we assume that neither the virtuality of the quark, p^, nor its transverse momen- 
tum squared, jpy- p, can be large in comparison with the hard scale Q^. Under these conditions, the 
quark momentum is soft with respect to the hadron momentum and its relevant component is xP^. 



2.1 Inclusive deep inelastic scattering 
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In Eq. ( [Z.lSp , neglecting terms which are \/Q suppressed, we can use an approximate expression 
for the delta function and 



2MW^^(^,/',5)^£4 j d2p^dp-dx2^5(x-Xg)Tr[4>(p,/',5)7^(i^ + ^ + m)7^ 



= E4^Tr 



Pq 



(2.18) 



where we introduced the integrated correlation function 



4>..(x,5)= j d^p^dp- ^/p,P,S) 



p+=xP^ 



(2.19) 



Notice that there is a contradiction between the fact that we assumed the transverse momentum 
of the quark to be small in comparison to the hard scale, yet we are integrating over the entire 
space of pj. Indeed, when dealing with transverse momentum of perturbative origin (i.e. arising 
from the radiation of gluons, see next section), which typically falls down as l/lpj]'^, we have to 
impose a cut-off on the maximum value the transverse momentum can reach. This cut-off depends 
on the scale Q^. On the other hand, the transverse momentum of nonperturbative origin, usually 
called intrinsic transverse momentum, is supposed to fall off very rapidly so that there is effectively 
almost no intrinsic transverse momentum above a typical scale of 1 GeV^. 

Finally, from the outgoing quark momentum, p + q, we can select only the minus component 
and obtain the final form for the hadronic tensor at leading twist 



1 



2MW^^(^,i',5) ^ £4 -Tr [^{xj,,S) 7^7+7^] 



(2.20) 



A few words to justify the last approximation are in order. The dominance of the minus component 
is most easily seen in the infinite momentum frame, where p + q^ is of the order of Q, while 
p^ + q'^ = 0, and Pj and m are of the order of 1 . However, if we perform a 1/2 expansion of the 
full expression, including the correlation function ^» [starting from Eq. ( [2.31[ )], we would be able 
to check that in any coUinear frame the dominant terms arise only from the combination of plus 
component in the correlation function and minus components in the outgoing quark momentum. 



2.1.3 One gluon additions 

Up to now, we took into consideration only quark-quark correlation functions at tree level. The 
addition of a gluon can either lead to the introduction of a quark-gluon-quark correlation function 
or can give rise to perturbative corrections to the photon-quark scattering In this thesis, we 
will not take quark-gluon-quark correlation functions into account, since they start contributing 
only at the twist- three level, and we will not examine perturbative corrections, since they give 
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Figure 2.4. One gluon addition: gluon coming from the soft blob. 



only origin to a logarithmic scale dependence of the quark-quark functions. Nevertheless, for 
completeness we will now give a sketchy view of these very important issues. 

When the gluons come directly from the soft blob, as in the diagram of Fig. Ol longitudinally 
polarized gluons (A+) are the dominant ones, while transversely polarized gluons {Aj.) are subject 
to a 1/2 suppression. If we choose a physical gauge where A+ = 0, then this kind of diagrams 
contribute only at twist three and higher and they require the introduction of quark-gluon-quark 
correlators [Q. On the other hand, in a different gauge the contributions of longitudinal gluons 
is present and is not necessarily suppressed by any power of \/Q. Then, we have to sum all the 
contributions with an arbitrary number of longitudinal gluons. The result of this summation can 
be cast in the form of a gauge link to be inserted in the definition of the quark-quark correlation 
function 

^/P^P^^^ = ji^J d'^e-^-^(/',5|v/,(^)^(^,0;path)v/,(0)|P,S> (2.21) 
where the gauge link is a path-ordered exponential 

=Sf(^,0;path) = ^exp (-i^.^^ ds^A^{s)^ . (2.22) 



with a straight path along the light-cone minus direction 

In A+ = gauges the link is equal to unity (although some subleties have been recently an- 
alyzed in Ref. We will henceforth neglect it and trade off manifest color gauge invariance 
for a lighter notation. Finally, we mention that recently it has been suggested that the gauge link 



could play an important role in the context of T-odd distribution functions [ |78| ] (see Sec. on 
page |2^). In particular, much care should be taken when including partonic transverse momentum. 
In this case, the gauge link path cannot simply run along the light-cone but has to have a transverse 
component and it might not be reducible to unity anymore [ |124| , [78| ]. 

Now we take a brief look to perturbative corrections to the quark-quark correlation function 
(see Refs. ^ and |147D . They are of two kinds: virtual gluon loop diagrams (Fig. ^3| ) and real gluon 
bremsstrahlung diagrams (Fig. Each virtual diagram contains ultraviolet and infrared diver- 
gences. The ultraviolet divergences can be cured using standard renormalization techniques. The 
infrared divergences cancel with analogous divergences in the real gluon emission diagrams 



2.1 Inclusive deep inelastic scattering 
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Figure 2.6. One gluon addition: real loop diagrams. 



The remaining parts of the diagrams give actual a, corrections to the tree level result of the 
previous section. In particular, coUinear divergences give origin to the leading-log part of the 
evolution equations, by which the parton distribution functions (see Sec. ^]2| on page ^1]) acquire 
a dependence on the scale [|TOO|,pl|,p|]. 



2.1.4 Leading twist part and connection with helicity formalism 

To identify the leading twist contributions to the cross section, it is convenient to define the pro- 
jectors 

^- = \r^r. (2.23) 

Before the interaction with the virtual photon, the relevant components of the quark fields are the 
plus components, i//^ = Xjr. They are usually referred to as the good components!^ Vice versa, 
after the interaction with the virtual photon, the relevant components of the outgoing quark fields 
are the minus components, = i//. Therefore, the leading twist part of the hadronic tensor 



'In the infinite momentum frame the good components are truly dominant and we can avoid the distinction between 
"quark" and "good quark". 
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in Eq. ( |2.18| ) can be projected out in the following way (see Fig. p7l)pl 



2P-q 



The differential cross section defined in Eq. 

d^cj 
dxg dy d 



L^vi^j'j^e) Tr 



2P-q 
1) takes the form 



Y^Y 



^4 



(2.24) 



2Pq 



Y^Y 



7' 



^4 



(2.25) 



where we explicitly showed Dirac indices (repeated indices are summed over). In Sec. [2.2. 1| on 
page we will see in detail how the insertion of the projectors effectively reduces the four- 
dimensional Dirac space into a two-dimensional subspace. Chiral-right and chiral-left good quark 
spinors can be used as a basis in this space. Therefore, it is possible to replace the Dirac indices 
with chirality indices (of good fields). By doing this, we put particular evidence on the connection 
with the helicity/chirality formalism (see e.g. Refs. |112| and [T^). In fact, the cross section can be 
conveniently rewritten as 



dxj^dyd^g 



do"i\xX 



dy 



(2.26) 



^Note that (1//+) = xj/ ^_ and (v/_) = v^^^^. 



2.2 The correlation function <I> 
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where the elementary electron-quark cross section is 




(2.27) 



where 



A{y) = l-y + ^-, C{y)=y(l-^y (2.28) 

Finally, we define the matrix F = ^<^^^>7+j , i.e. the transpose of the leading-twist part of 
the correlation function, and we observe that 



(2.29) 

Thus, the transpose of the correlation function describes the forward scattering of a good antiquark 
off a hadron, or equivalently the forward scattering of an antiquark off a hadron in the infinite 
momentum frame. As any scattering matrix, for any antiquark-hadron state |a) the expectation 
value (^a \ M |a) must be positive. In mathematical terms, this means that the matrix is positive 
semidefinite, i.e. the determinant of all the principal minors of the matrix has to be positive or zero. 
This property will prove to be essential in deriving bounds on the components of the correlation 
function, i.e. the parton distribution functions. 



2.2 The correlation function O 



As shown in Eq. ( [2.16aD , the quark-quark distribution correlation function, can be expressed in 
terms of bilocal operators. At leading order in l/Q, it contains all the relevant information about 
the nonperturbative dynamics of the quarks inside the hadron. Due to its nonperturbative nature, 
it is not possible to calculate it from first principles, as we don't know how the hadronic states are 
built up from the elementary quark and gluon fields. 

When considering subleading orders ma l/Q expansion, quark-gluon-quark correlation func- 
tions have to be considered, as we briefly mentioned in Sec. [2.1.3| on page |17[ In this case, the 
general structure of the hadronic tensor becomes richer. In the rest of the thesis, as the analysis will 
be concerned only with leading order terms, we will not deal with quark-gluon-quark correlation 
functions. 
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To get more insight into the information contained in the correlation function, which is a Dirac 
matrix, we can decompose it in a general way on a basis of Dirac structures. Each term of the 
decomposition can be a combination of the Lorentz vectors p and P, the Lorentz pseudovector S 
(in case of spin-half hadrons) and the Dirac structures 

1, 75, 7^75, ic^^''75- 

The spin vector can only appear linearly in the decomposition (cf. Eq. ( [2.43P ). Moreover, each 
term of the full expression has to satisfy the conditions of Hermiticity and parity invariance 

Hermiticity: ^{p,P,S) = f^\p,P,S)f, (2.30a) 

parity: ^{p.P.S) = f^{p,P, -5) f (2.30b) 

where = 5^^p^ and so forth for the other vectors. The most general decomposition of the 
correlation function ^» imposing Hermiticity and parity invariance is [ |151| , |141[ ] 

^{p,P,S)=MA,l+A2f+A^^ + ^a^yP^'p''+iA,p■Sr5 

+MAjr^ +A, ^IJ>y, +A, ^]^r,+iA,a,yr,S^P' (2.31) 

+ \A a vS^'n^'+lA vP^'n''+A ^i'^P''^^^'' P'' 

+ ^\O^HvJs'^ P +^11 ^2 ^l^vh^ P ^^^12 ' 

where the amplitudes A ■ are dimensionless real scalar functions A - = A-(p ■ P, p^) . 

The correlation function can be separated in a T-even part and a T-odd part, according to the 
definition 

^T-evenb,^,5) = '^yV ^T-.^eniP ^ P . S) ijV , (2.32a) 

^T-oMiP:P^S) = -i7y ^T-odd(P'A^)i7y- (2.32b) 

Thus, the terms containing the amplitudes A^, A^ and can be classified as T-odd. 

At leading twist, we are interested in the projection 4>(jCg, 5)7+. After inserting the general 
decomposition of Eq. ( |2.31[ ) into Eq. ( |2.19| ), we can project out the leading-twist components and 
obtain the general expression 



j^+^{x,s)r+ = {Mx)+s^g^ix)r,+h^{x)r,^^) (2.33) 

where we introduced the parton distribution functions 
/i {x) = J d^pj dp^ d{2p-P)5{p^+ x^M^ + p^ - 2xp ■ P) [A^ + xA^] , (2.34a) 
g,{x)=J d^pj. dp^ d{2p ■ P) 5 {p\ + xV + p2 _ 2xp .P) -A^-(^-^ {A, + Mg) 



h^{x) = j d^pj^dp^d{2p-P)5{pl+x^M^ + p^-2xp-P) 



Ag MjQ-F^^Ajj 



(2.34b) 
. (2.34c) 



The correlation function ■ 
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The function /j is usually referred to as the unpolarized parton distribution, and it is sometimes 
denoted also as simply f ox q (where q stands for the quark flavor). The function is the parton 
helicity distribution and it can be denoted also as A/ or A^. Finally, the function h-^ is known as 
the parton transversity distribution; in the literature it is sometimes denoted as 8q, IS.jq or A^./, 
although in the original paper of Ralston and Soper QISIQ it was called hj. In this thesis, we will 



follow the nomenclature suggested by Jaffe and Ji [ |116| ] and later extended in Ref. |141| , because it 
allows a harmonious connection with the most general cases of transverse momentum dependent 
distribution functions, as we shall see later. A thorough discussion on the different naming schemes 
is presented in Ref. [pi]] . 



The individual distribution functions can be isolated by means of the projection 

<I>[r] ^ lTr(4)r), (2.35) 

where T stands for a specific Dirac structure. In particular, we see that 

/iW=^t^^'W, (2.36a) 

^iW=^t^^^'W, (2.36b) 

/z^(jc)=4>['^''^^](x). (2.36c) 



2.2.1 Correlation function in helicity formalism 



We will now examine how it is possible to write the correlation function as a matrix in the chirality 
space of the good quark fields ® the spin space of the hadron. We will work out the steps in a 
meticulous way, even if we will incur the risk of introducing some redundant steps. 

The correlation function is a 4 x 4 Dirac matrix. However, due to the presence of the projec- 
tor on the good components of the quark fields, the leading-twist part spans only a 2 x 2 Dirac 
subspace. This is evident if we express the Dirac structures of Eq. ( P.33| ) in the chiral or Weyl 
representation. Using this representation, the correlation function reads 



^+4>(x,5)7^ 



lf,{x)+Si^g,{x) {S,-iSy)h,{x)\ 





\{S, + iSy)h,{x) f,{x)-Si^g,{x)j 



(2.37) 



As shown by this explicit form, it seems that the four-dimensional Dirac space can be reduced to 
a two-dimensional space, retaining only the nonzero part of the correlation function. The relevant 
part of the Dirac space is the one corresponding to good quark fields. To show this explicitly, 
we introduce the chiral projectors ^j^jj^ = (1 ± 7g)/2 and define the good chiral-right and good 
chiral-left quark spinors, i.e. the normalized projections 



(2.38) 
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Then, we can define a new matrix in the chirality space of the good quark fields 

(^+4>(x,5)r+)^,^^ ^ (^+4>(x,5)r+)^.^.<,, (2.39) 
Any contraction with bad quark fields vanishes. Explicit computation of the matrix elements yields 

(^+4.(^,5)7+)^^ = u^^^[0^^^{x,S)y+) = f,{x) +S^g,ix), (2.40a) 

(^+4>(x,5)r+)^^ = <^ (^+4.(^,5) 7+) ^..^^^^ = Mx) -S^g,{x), (2.40b) 

(^+4>(x,5)r+) =M%(^+4>(x,5)7+) M^^ = (5,-iS,)/zi(x), (2.40c) 

(^+4>(x,S)r+) = m;J^+4>(x,5)7+) = (5, + iS,)/zi(x). (2.40d) 

The correlation matrix in the good quark chirality space is then 

^+4>(x,S)7+) = ^ ^^^^ ^ -'^ M . (2.41) 

^^1^1 V(5. + i5,)^(x) /i(x)-5^^i(..); 

As we could have expected, this result corresponds simply to taking the full Dirac matrix in Weyl 



representation, Eq. ( [2.37[ ), and stripping off the zeros. From the matrix representation in the chi- 
rality space it is clear why the function is defined to be chiral odd. 

The correlation function is a matrix in the parton chirality space and depends on the target spin. 
By introducing the helicity density matrix of the target 

P(\^, =^(1+^-^). .,=^1 ■ " ^ -'h (2-42) 

we can obtain the correlation function from the trace of the helicity density matrix and a new 
matrix in the quark chirality space ® the target spin space: 

(^+4>(x,5)7+) =p{S)^^ (^+$(x)7+)^;^'. (2.43) 
We will refer to the last term of this relation as the matrix representation of the correlation function 




or, more simply, as the correlation matrix. Fig. shows pictorially the position of the spin 
indices. 

Starting from Eq. ( [2.33P and using the relation 



A'lAi 



^f/+5A + ^-.^x + S,^, = p(5)^ ^ I 'I (2.44) 



The correlation function ■ 
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Figure 2.8. Illustration of the position of the indices of the correlation matrix. 



we can cast the correlation function in the matrix form 



A'lAi 



(/l (^) + ^1 (■^) 75)^+ K (-^) (r. - i7y) 75 
^iW (7x + i7y) 75.^+ (/iW -^lW75) ^4 



(2.45) 



Finally, by expressing the Dirac structures in Weyl representation and reducing the Dirac space as 
done before, we obtain the matrix representation of the correlation function 



(a 

















2h^{x) 







2h^{x) 


f,{x)-g,{x) 


\ 








/iW+^iW 



\ 



(2.46) 



where the inner blocks are in the hadron helicity space (indices AjAj), while the outer matrix is in 
the quark chirality space (indices x[X\)- 

The form of the correlation matrix can also be established directly from angular momentum 
conservation (requiring A!^+x[ = ^\+Xi) and the conditions of Hermiticity and parity invariance. 
In matrix language, the condition of parity invariance consists in Ql 190 



A',A, 



X\ X\ 



(2.47) 



The most general form of the correlation matrix complying with the previous conditions corre- 
sponds to Eq. ( [2.46p . 

As mentioned at the end of Sec. [2.1.4| on page [1^, with transposing the quark chirality indices 
of the correlation matrix we obtain the scattering matrix 



X\X\ 



/iW+<?iW 








lh^{x) 



















/iW-<?iW 





lh^{x) 








/iW+,?iW 



(2.48) 



Note that because of the inversion of the quark indices, the lower left block has x[ = X\ = ^ 
and vice versa for the upper right block. Since this matrix must be positive semidefinite, we can 
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Figure 2.9. Probabilistic interpretation of the unpolarized distribution function /[ (a), and of the helicity 
distribution function g ^ (b). 

readily obtain the positivity conditions 



/iW>0, (2.49a) 
\8i{x)\<f,{x), (2.49b) 
|/ii(x)|<K/iW+^iW)- (2-49C) 



The last relation is known as the Soffer bound [ [159^ . 



2.3 The transversity distribution function 

In the previous section, we wrote the forward antiquark-nucleon scattering matrix, F, in the helicity 
basis of the hadron and of the good quark (to be more precise, we used the chirality basis for 
the quark). Each entry, with indices Xi-^'i^ describes the product of the amplitude for the 
scattering of an antiquark with helicity (chirality) x[ off a hadron with helicity Aj going to anything 
times the conjugate of the amplitude for antiquark with helicity Xi off a hadron with helicity A'^ 
going to anything. 

F{xf^^: o^l^f d'P^ {X\ iw^),^ \PA\TiX\ {w+),> \PA,). (2.50) 



In this basis, the probabilistic interpretation of the functions and is manifest, since they oc- 
cupy the diagonal elements of the matrix and they are therefore connected to squares of probability 
amplitudes (see Fig. p!P| ) 

/i W = \ {F{x)il + F{x)u) g,{x) = \ (F(4|-F(4i) (2.51) 

On the other hand, the transversity distribution is off-diagonal in the helicity basis. This means 



2.3 The transversity distribution function 
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(b) 

Figure 2.10. Probabilistic interpretation of the transversity distribution function h^, in the helicity basis (a) 
and in the transversity basis (b). 

that it does not describe the square of a probability amplitude, but rather the interference between 
two different amplitudes [see Fig. |2.10| (a)] 



(2.52) 



The transversity distribution recovers a probability interpretation if we choose the so-called 
transversity basis, instead of the helicity basis, for both quark and hadron [ |1 19| , |1 ISfl . The transver- 
sity basis is formed by the "transverse up" and "transverse down" states. They can be expressed in 
terms of chirality eigenstates 



The same relation holds between the hadron transversity and helicity states. 
In the new basis, the scattering matrix takes the form 



(2.53) 



















+ h^{x) 







/i(^) 


-h^{x) 




-h^{x) 














-/jj(x) 


/i(^) 


-hiix) 























+ h^{x) 



(2.54) 



and clearly the transversity distribution function can be defined as [Fig. |2.10| (b)] 

^W = i(FWjj-F(x)n). 



(2.55) 



The transversity distribution has been the object of several model calculations, using e.g. the 
bag model [ ]1 16| ], the spectator model [ |120| ], the chiral soli ton model 01480 and others pO|]. The 
integral of /ij - also known as the tensor charge of the nucleon - has been evaluated in lattice 
QCD [ p4] ] . A recent review on the transversity distribution is presented in Ref . 

The transversity distribution evolves with the energy scale in a different way as compared to 
the helicity distribution, without mixing with gluons [PS1,|5^,PS|,|132|,|109[]. 
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2. Distribution functions and transversity 



Transverse direction 



Thomas precession orbit 





'uark "spin' 



Longitudinal direction 



(a) (b) 

Figure 2.11. Thomas precession effect on transversity (a) and helicity (b). 

2.3.1 Transversity and Thomas precession 

In the literature, it is common to find the statement that the difference between the helicity distri- 
bution and the transversity distribution is connected to relativistic effects, since boosts and rotation 
do not commute 0114| , |119| , |117| , |1130 . Relativistic effects influence observable quantities depending 
on the dynamics of the system and they can therefore give important information about its struc- 
ture. Therefore, the difference between helicity and transversity distributions can shed light on the 
structure of the nucleon and its spin. 

This statement can be understood even in the framework of purely classical relativistic mechan- 
ics. We will show how relativistic effects can change an observable in a toy model with just a little 
bit of dynamical complexity. Of course, the model is not meant to describe a nucleon. We will 
consider "spin" merely as a pseudovector attached to the quark. "Helicity" will be the projection 
of the spin along the momentum of the quark and "transversity" will be the projection transverse 
to the momentum of the quark. Note that these quantities are only reminescent of the real helicity 
and transversity. 

Suppose we have a system constituted by a quark revolving in a circular orbit with its spin 
aligned in the direction of the orbital axis. To have an analogy of the helicity distribution, we 
boost the system to a velocity v along the direction of the orbital axis, and we measure what is the 
probability of finding the helicity of the quark aligned along the axis direction. For the transversity, 
we boost the system to a velocity v transverse to the axis direction, and measure the probability of 
finding the transversity of the quark aligned along the axis direction. 

To properly deal with this situation in a relativistic way, we have to take into account Thomas 
precession, an effect which occurs whenever a frame of reference (in our case joined to the quark) 
is moving with a velocity v with respect to the observer and is at the same time subject to an 
acceleration a Qlll| , |165[ ]. Thomas precession causes the quark frame of reference (and the spin of 
the quark with it) to precess with an angular velocity 



V X a 



(2.56) 



2.4 Inclusion of transverse momentum 
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In the infinite momentum frame, the speed of revolution of the quark is negligible with respect 
to the overall velocity of the system, so that the quark velocity is approximately v. However we 
cannot neglect its centripetal acceleration a. 

Let us first analyze the situation with transversity [Fig. [2.1 1| (a)]. The cross product of the 
velocity of the quark and the centripetal acceleration is pointing in the transverse direction. Thomas 
precession will not influence the orientation of the spin of the quark. The transversity distribution 
is just 1. 

The situation is dramatically different for helicity [Fig. [2.1 1| (b)] . The cross product of the quark 
velocity and its acceleration is still pointing in a transverse direction, which is now orthogonal to 
the spin of the quark. This causes the spin of the quark to precess around a transverse axis. The 
net helicity of the quark will be zero and so will be the helicity distribution. 

In a simpler system, for instance if the quark would not have any orbital angular momentum, 
helicity and transversity distributions would be the same. In conclusion, from this example we 
see that, due to relativistic effects, the difference between two apparently similar observables can 
reveal something important about the structure of a system. 



2.4 Inclusion of transverse momentum 



So far we have been concerned with the integrated correlation function, defined in Eq. ( [2.19[ ), 
the only relevant one in totally inclusive deep inelastic scattering. In the next chapters, we will 
analyze also semi-inclusive scatterings, where we will need to consider the transverse-momentum 
dependent version of the correlation function, i.e. 



4>..(x,Py,5) = J dp ^j.ip^P.S) 



p+=xP+ 



(2.57) 



Starting from the general decomposition presented in Eq. ( P.31[ ), the leading order part of the 
transverse-momentum dependent correlation function becomes 



Mx,pj-) +ihf{x,pl) ^+S^^g,^{x,p'T) r5+Sj^htt^{x,p'T) 75 ^ 



M 



M 



75 



(2.58) 



The definition of the parton distribution functions in terms of the amplitudes A-, introduced in 
Eq. ( |2.33[ ), can be found elsewhere [ ]135[ , |164| , [120| ]. 
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2. Distribution functions and transversity 



For any transverse-momentum dependent distribution function, it will turn out to be convenient 
to define the notation 



2M 



f{X:P 



(iL)" 

V2M2 J 



(2.59a) 
(2.59b) 



for n integer. We also need to introduce the function 



h^{x,p\)=h^.j^{x,p\)+h^^^\x,p\ 



(2.60) 



The connection with the integrated distribution functions defined in Eq. ( |2.34D is 



fi{x)= J d^pj f^{x,p\), 
g^{x)= J d^pjg^j^{x,pj^), 
h^{x)= j (fp-j^h^ix.p^). 



(2.61a) 
(2.61b) 
(2.61c) 



Note that the distribution functions hj- and are T-odd. At first, this class of functions 



was supposed to vanish due to time-reversal invariance [ |151| ]. Sivers [ |157| , |158| ] was the first one to 
consider an observable arising from the T-odd distribution function f^, since then called the Sivers 
function. The complete analysis of leading-twist T-odd distribution functions and the introduction 
of hj- were carried out by Boer and Mulders . A proof of the nonexistence of T-odd distribution 
functions was suggested by Collins in Ref. ^ but recently it has been repudiated by the same 
author [ [TSD after Brodsky, Hwang and Schmidt JT^ explicitly obtained a nonzero Sivers function 
in the context of a simple model. The question at the moment awaits clarification, but it is likely 
that T-odd distribution functions will stir a lot of interest in the near future, together with T-odd 
fragmentation functions, of which we shall abundantly speak in the next chapters. 

As done in the previous section, we can express the transverse momentum dependent corre- 
lation function as a matrix in the parton chirality space ® target helicity space. To simplify the 
formulae, it is useful to identify the T-odd functions as imaginary parts of some of the T-even func- 
tions, which become then complex scalar functions. The following redefinitions are required:^ 



8^T + '^fi 



IT 



IL- 



(2.62) 



^From a rigorous point of view, it would be better to introduce new functions, e.g. g^j and hfj^, but this would 
overload the notation. 
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The resulting correlation matrix is p^,p3|] 







M 




M 





/i -<? 



IL 



\ 



l/'rl 



^e'^/'/i^-; 
M 1^ 

M2 ^ 



i^e'^p/ij, 
M 



f\-S\L 

\Pt 



M 



fl+SlL 



J 



(2.63) 

where for sake of brevity we did not explicitly indicate the x and dependence of the distribution 
functions and where 0p is the azimuthal angle of the transverse momentum vector. 

The distribution matrix is clearly Hermitean. Notice that by introducing the transverse momen- 
tum of the quark, the angular momentum conservation requirement becomes less constraining and 
we can have non zero contributions in all the entries of the scattering matrix. To be more specific, 
when an exponential e^' '^p appears in the matrix, we have to take into account /' units of angular 
momentum in the final state. The condition of angular momentum conservation becomes then 
Aj + + ^' = -'^i + Zi • Also the condition of parity conservation is influenced by the presence of 
orbital angular momentum and becomes 



F(..,p^)^i^i = (-l)'V(x,/,^)-^i-^; 



(2.64) 



The positivity of the matrix is not hampered by the introduction of the transverse momentum 
dependence, since 



F(x,p^,5)^i^; 

X\X\ 



^^-Ake-P-^(P,A;|(v/+)t^(^)(V/+^^^ 



XT / • 



Bounds to insure positivity of any matrix element can be obtained by looking at the one-dimensional 
and two-dimensional subspaces and at the eigenvalues of the full matrix.0 The one-dimensional 
subspaces give the trivial bounds 



/l(x,/>2)>0, 



\g^^{x,p\)\<f^{x,pl) 



(2.66) 



Cf. Ref. |13q for an earlier discussion on positivity bounds for transverse momentum dependent structure func- 
tions. 
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2. Distribution functions and transversity 



From the two-dimensional subspaces we get 

|^i|<^(/i+^il)</i. 



IT 
.(1) 



1 



< 



Pt 

.2 



(/i+^il) (/i-^il)<(/| 



(1/2) 



(2.67a) 
(2.67b) 
(2.67c) 
(2.67d) 



where, once again, we did not explicitly indicate the x and pj dependence to avoid too heavy a no- 
tation. Besides the Soffer bound of Eq. ( [2.67aD , now extended to include the transverse momentum 
dependence, new bounds for the distribution functions are found. 

The positivity bounds can be sharpened even further by imposing the positivity of the eigen- 
values of the correlation matrix. The complete analysis has been accomplished in Ref. ^(see also 
Ref. ph. 



2.5 Cross section and spin asymmetries 



The cross section for inclusive deep inelastic scattering at leading twist is expressed by Eq. ( [2. 26] ). 
Extracting the target helicity density matrix as done in Eq. ( p.43[ ) and using the matrix F, the 
equation becomes 



d^a 



dxg dy d 



F(jc)^i^i 



dcj^'?\ziz; 1 

dy ) 2' 



(2.68) 



Inserting the expressions of Eq. ( |2.27p , ( [2. 42] ) and ( [2.48[ ) into the previous expression leads to the 
result 



d^a 



2a^ 



dxg dy d sxgy^ 



y 



l-y + '-]f,{xj,)+XeS^(y 



2 



8Mb) 



(2.69) 



where the index q denotes the quark flavor. 

The transversity distribution does not appear in the cross section for totally inclusive deep in- 
elastic scattering at leading twist. The reason is that it is a chiral odd object and in any observable 
it must be connected to another chiral odd "probe". In inclusive deep inelastic scattering, what 
probes the structure of the correlation function is the elementary photon-quark scattering, which 
conserves chirality. We will see in the next chapters how semi-inclusive deep inelastic scattering 
provides the necessary chiral odd partners for the transversity distribution (i.e. chiral odd fragmen- 
tation functions). 

The r.h.s. of Eq. ( |2.69D is independent of the azimuthal angle ^g, so that we can integrate the 
cross section over this angle. The result is 



dxg dy sx^y^ 



y 



r 

2 



sMb) 



(2.70) 



2.6 Summary 
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If the spin of the target is oriented along the direction of the electron beam, it will have its 
longitudinal component oriented along the — z direction, that is will be negative. On the contrary, 
orienting the spin of the target in the opposite direction will produce a positive 5^. Summing the 
cross sections obtained with opposite polarization, we isolate the unpolarized part of the cross 
section 

d^a,, ^ \ (d2a_ + d^a^) - ^ (l -3^ + ^) I^/f (2.71) 

The first subscript indicates the polarization of the beam, while the second stands for the polariza- 
tion of the target. The letter U stands for unpolarized. The right arrow means polarization along 
the beam direction, the left arrow means the opposite. Subtracting the cross section we obtain the 
polarized part 

d^cT,, = \ (d2a_ - d^cT^) ^ ^ I A.I (, - ^) Zelg\{x,). (2.72) 

where now the subscript L specifies that longitudinal polarizations of beam and target are required. 
We can define the longitudinal double spin asymmetry 



Note that, since the beam direction does not exactly correspond to the virtual photon direction, 
the degree of longitudinal polarization of the target will be somewhat smaller than the degree of 



polarization along the beam direction, while a small transverse polarization will arise 0145| , |1280 
These effects are anyway \/Q suppressed. 



2.6 Summary 

In this chapter, we introduced the hadronic tensor, containing the information on the structure of 
hadronic targets in deep inelastic scattering [cf. Eq. (|2.6p]. We studied the hadronic tensor in the 
framework of the parton model at tree level and we came to the introduction of a quark-quark 
correlation function [cf. Eqs. ( |2.19[ ) and ( [2.20[ )]. The correlation function can be parametrized in 



terms of parton distribution functions. In particular, at leading order in 1/2 (leading twist) we 
introduced the unpolarized distribution function, /j, the helicity distribution function, g^, and the 
transversity distribution function, [Eq. ( |2.33| )]. 

We demonstrated how the leading twist part of the correlation function can be cast in the form 
of a forward scattering matrix [cf. Eq. ( [2.48[ )]. Exploiting the positivity of this matrix, we derived 
relations among the distribution functions [Eq. ( [2.49D ]. We also discussed the probabilistic inter- 
pretation of the distribution function, with a particular emphasis on the transversity distribution. 

We repeated the analysis of the correlation function introducing partonic transverse momen- 
tum. In this case the decomposition of the correlation function contains eight distribution func- 
tions [Eq. ( [2.58[ )]. We have seen that each entry of the corresponding scattering matrix is nonzero 
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2. Distribution functions and transversity 



[Eq. ( |2.63| )], indicating that the full quark spin structure in a polarized nucleon is accessible if 
transverse momentum is included. The connection with the helicity formalism and consequently 
the extraction of positivity bounds on transverse momentum dependent distribution functions are 
among the original results of our work. 

Finally, we expressed the cross section for inclusive deep inelastic scattering at leading order 
in 1/2 in terms of distribution functions [Eq. ( |2.69| )]. We showed that the helicity distribution 
can be accessed by measuring the longitudinal double spin asymmetry, but we concluded that the 
transversity distribution is not accessible in inclusive deep inelastic scattering. 

In the next chapter, we are going to see how the transversity distribution can be measured 
in one-particle inclusive deep inelastic scattering, in combination with the Collins fragmentation 
function. 
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Fragmentation functions 
and the Collins function 



Gather up the fragments that remain, that 
nothing be lost. 

John 6:12 



In the previous chapter we analyzed totally inclusive deep inelastic scattering as a means to probe 
the quark structure of the nucleons. We concluded that some aspects of this structure - notably 
the transversity distribution and the transverse momentum distribution - are not accessible in this 
kind of measurement. It is therefore desirable to turn the attention to a more complex case, that 
of one-particle inclusive deep inelastic scattering, where one of the fragments produced in the 
collision is detected. In this case, we will need to introduce some new nonperturbative objects, the 
fragmentation functions . 



3.1 One-particle inclusive deep inelastic scattering 

In one-particle inclusive deep inelastic scattering a high energy electron collides on a target nucleon 
via the exchange of a photon with a high virtuality. The target breaks up and several hadrons are 
produced. One of the produced hadrons is detected in coincidence with the scattered electron (see 
Fig. ^7T] on the next page). As a result of the hardness of the collision, the final state should 
consist of two well separated "clusters" of particles, one is represented by the debris of the target, 
broken by the collision, and the other is represented by the hadrons formed and ejected by the hard 
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3. Fragmentation functions and the Collins function 




Figure 3.1. One-particle inclusive deep inelastic scattering. 



interaction with the virtual photon [ ]168[ ]. The former are called target fragments and the latter 
current fragments. We will take into consideration only events where the tagged final state hadron 
belongs to the current fragments. 



3.1.1 Kinematics 

As for inclusive scattering, we denote with I and I' the momenta of the electron before and after 
the scattering, with P, M and S the momentum, mass and spin of the target nucleon, with q the 
momentum of the exchanged photon. Moreover, we introduce the momentum and the mass 
of the outgoing hadron. 
We define the variables 



IPq 



y- 



pV 



PPh 

Pq 



(3.1) 



We will assume that the following conditions hold: Q^,P ■ q,P- P^ ^ M^,M^, and Xg.Zj^ fixed. 

For a distinction between the current and the target fragments, the rapidity separation should 
be taken into account. This separation is quite clear for high values of z^. For low values of z^, 
separating the two clusters becomes arduous. In events with higher W^, the limit of Zf^ can be 
pushed lower [ |140| ]. 



In the frame of reference defined in Sec. |2.1.1| on page |T2|, the target momentum and the virtual 
photon momentum have no transverse components. The outgoing hadron's momentum can be 
parametrized as 

'z,Q A{Ml + \P^^\^) ^ 



h 



(3.2) 



To avoid the introduction of further hard scales, it is required that Pf^ 

The frame of reference we adopted is the more natural one from the experimental point of view, 
as the longitudinal direction corresponds (up to order 1/2) to the beam direction and the hadron's 
transverse momentum corresponds to what it is actually measured in the lab. However, as it will 
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become clear later, in order to preserve a symmetry between the distribution and fragmentation 
functions, it is convenient to use a different frame of reference where the target and outgoing 
hadron momenta are coUinear, while the photon acquires a transverse component. To distinguish 
between the two frames of reference when needed, from now on we will use the subscript T when 
denoting a transverse component in the new frame, while we will use the subscript _L to denote 
transverse components in the former frame. 
Li the T frame, the external momenta are 

f2 



h 



AQ_ 

'z,Q AMI 



2P^ 



P^, 



■ Q A{Q^-\q^\^) 



p 



(3.3a) 
(3.3b) 
(3.3c) 



with P^ f» Q^zJ {2P^Xg). The connection with the transverse momentum components of the 
photon in the T frame and of the outgoing hadron in the ± frame is [ |141| ] 



It 



(3.4) 



3.1.2 The hadronic tensor 



The cross section for one-particle inclusive electron-nucleon scattering can be written as 



2E, d^a 



d^P^ dxg dy d 0^ Isx^Q- 



■L^y{lj',Xe)2MW^''iq,P,S,Pf 



or equivalently as 



a' 



■L^y{lj'X)1MW^\q,P,S,P^). 



(3.5) 



(3.6) 



To obtain the previous formula, we made use of the relation d^PjJlEj^ dz^ d^Pj^^/lzj^. 
The hadronic tensor for one-particle inclusive scattering is defined as 



H^\P,S,P^,,P^) = (/',5|7^(0) IP,, JT'Xp,, ^'\riO) \P,S). 



(3.8) 



By integrating over the momentum of the final-state hadron and summing over all possible hadrons, 
we recover the hadronic tensor for totally inclusive scattering 



d^P 



1 



2^,^ 2MW^^{q,P,S,P,) = ^Zj^ 2K5^'\q + P-P,^)H^^{P,S,P, 



(3.9) 



= 2MW^''{q,P,S), 
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Figure 3.2. The bull diagram, describing the hadronic tensor at tree level. 

where now the state ^ indicates the sum of the states and h, and = P^, +Pj^- Note that 
we did not include any dependence of the hadronic tensor on the polarization of the final state 
hadron. The reason is that this polarization is usually measured through the decay of the final state 
hadron into other hadrons. In this sense, this case falls within the context of two-particle (or even 
three-particle) inclusive scattering, as we will show in the next chapter. 

In the spirit of the parton model, the virtual photon strikes a quark inside the nucleon. In the 
case of current fragmentation, the tagged final state hadron comes from the fragmentation of the 
struck quark. The scattering process can then be factorized in two soft hadronic parts connected 



by a hard scattering part, as shown in Fig. p^ . 

We need to introduce a parametrization for the vectors 



2xP+ 
P^ z{k^+M 
z ' 2P 



k 



(3.10a) 



(3.10b) 



Without expliciting including the antiquark contributions, the hadronic tensor can be written at tree 
level as 

2MW^\q,P,S,P,)=Y,e] J d^ 5^ (p + ^ - fc) Tr ($(;?,/',5) 7^ A(fc,/',) /) , (3.11) 



where 4> is the correlation function defined in Eq. ( ]2.16a| ) and A, 
1 



{my 



d^^ e*^(0|i/A,(^)|P,>(/',|V7,(0)|0> 



d-^Pv 



(2;r)3 2P0 



Y / ' 



(3.12) 
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is a new correlation function we need to introduce to describe the fragmentation process [^T]]. 
Neglecting terms which are 1 / Q suppressed, we obtain the compact expression 



where, as we shall do very often, we used the shorthand notation 



(3.13) 



and where we introduced the integrated correlation functions 

^{x,Pj^,S) = I dp-^{p,P,S) 

A(z,M = ^/ dk+A{k,P, 



(3.1 



) 



(3.15a) 
(3.15b) 



Often in experimental situations the transverse momentum of the outgoing hadron is not mea- 
sured. This corresponds to integrating over the outgoing hadron's transverse momentum, P^^^. As 
we already pointed out, this transverse momentum squared has to be small compared to 2^. Ex- 
perimentally, this can be insured by imposing a cut-off on the data. Due to the 5 distribution in 
Eq. ( [3.14D , this condition implies also {jpj + kjY <^ 2^- Considerations as the one discussed in 
Sec. [2.1.2| on page [T3|, can be applied to the perturbative and intrinsic components of the transverse 
momenta. 

The integration of the cross section yields 
d^a 



05" 



where 



dxj^dy dz^^d(^^ 4z,^sxj^Q 



4>(x,5) = J dp d^Pj.^{p,P,S) 
A{z)^ydk+d%Aik,P,^) 



p+=xP+ 



(3.16) 

(3.17a) 
(3.17b) 

(3.17c) 



3.1.3 Leading twist part and connection with helicity formalism 



Using the projectors and defined in Eq. ( [2.23D , we can isolate the leading order part of 
the hadronic tensor, analogously to what was done in Sec. p. 1.4] on page [T^ 



2MWt'''{q,P,S,PjJ^4zj^y 



4z,J^ 



Tr(^+4>(xg,/>^,S) ^_ 7^ ^-A{zM ^+ f) 
Tr (^^^{xs,p^,S) 7+ ^ ^_ ^-HzM r ^ ^4 



(3.18) 
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The differential cross section becomes 



r T 



777 ^uv (^5^ j'^e 



2 



Y+Y^ 



ji 



Y^Y 



^4 



(3.19) 



As we already discussed in Sec. [2.2. 1| on page ^ the restriction to the leading order allows us to 
reduce the four dimensional Dirac space to the two-dimensional subspace of good fields. Writing 
all the components of the cross section in the chirality space of the good fields, we obtain 



dV 



dxj^dydz^d^^d^P^^ 



= PiS) 



a,a; 



x'lXi \ dy J 



y, yi V 



do'"^\XiXvX2X2 



Xix[ ^ 



X2X2 



(3.20) 



The elementary electron-quark scattering matrix is 



do'''i\XiXvX2X2 cc^el 



dy 



Y Y^ 



X\X2 /'y^y^ 



X'2X'\ 





/ A{y) + XeC{y) 








-Biy) \ 
















sxj^y^ 
























A{y)-XeC{y) ) 



(3.21) 



where 



A{y) = \-y + - 



B{y) = {\-y), C{y)=y(\ 



(3.22) 



The internal blocks have indices X\ X\ and the outer matrix has indices X2 Xi- Notice the difference 
between Eq. ( [3.2ip and Eq. ( [2.27[ ): in the latter, the presence of the Kronecker delta was identifying 
two of the chirality indices, thus reproducing a true scattering matrix. Here, we don't have just a 
quark line connecting the two scattering amplitudes, but rather the correlation function A. For this 
reason, we cannot identify the chirality indices of the outgoing quark. Strictly speaking, this is not 
a scattering matrix anymore, but a scattering amplitude times the conjugate of a different scattering 
amplitude [|T9p. However, for conciseness we follow the notation of, e.g., Ref. |117 . 



Im 
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The leading order part of the correlation function corresponds to the transition probability 
density for the process of a quark decaying into hadrons 



Y 

(3.23) 

Analogously to the distribution case, the fragmentation correlation matrix is positive semidefinite, 
allowing us to set bounds on the fragmentation functions. 



3.2 The correlation function A 

While the distribution correlation function describe the confinement of partons inside hadrons, the 
fragmentation correlation function describes the way a virtual parton "decays" into a hadron plus 
something else, i.e. q* — > hY . This process is referred to as hadronization. It is a clear manifestation 
of color confinement: the asymptotic physical states detected in experiment must be color neutral, 
so that quarks have to evolve into hadrons.0 Nowadays, we can count on reliable phenomenological 
descriptions of hadronization, such as the Lund model. On the other side, understanding it from 
first principles, as well as including spin degrees of freedom, is very difficult. 

As on the distribution side the quark-quark correlation function is sufficient to describe the 
dynamics at leading order in l/Q, also for fragmentation it is sufficient to consider quark-quark 
correlation functions. 

The procedure for generating a complete decomposition of the correlation functions closely 



follows what has been done on the distribution side in Sec. on page gl]. It is necessary to 
combine the Lorentz vectors k and Pj^ with a basis of structures in Dirac space, and impose the 
condition of Hermiticity and parity invariance. The outcome is 

A(fc, P,) = 5i 1 + f^+B^f+^ o^,P;^k\ (3.24) 

h 

The amplitudes 5- are dimensionless real scalar functions 5- = B -{k ■ Pj^,k^) . The T-even and T-odd 
part of the correlation function A can be defined in analogy to Eqs. ( |2.32[ ). According to those 
definitions, the last term can be classified as T-odd. 

At leading twist, we are interested in the projection A(z,A;j,) 7^. The insertion of the de- 
composition given in Eq. ( |3.24| ) into Eq. ( |3.15b| ) and the subsequent projection of the leading-twist 



'Note that on the way to the final state hadrons, the color carried by the initial quark can be neutralized without 
breaking factorization, for instance via soft gluon contributions. 
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component leads to 



where we introduced the parton fragmentation functions 

1 r { 



2k -R 



H Bn 

z 



' " i / d(2^ . P,) 5 (4 + ^ + - ^ ) [-^4] ■ 



The fragmentation function Hf^ is known with the name of Collins function [[77|]. 

The individual fragmentation functions can be isolated by means of the projection^ 

A[r] =Tr(Ar), 

where F stands for a specific Dirac structure. In particular, we see that 

D^{z,z^kj)=A^^'\z^r)^ 



(3.25) 



(3.26a) 
(3.26b) 



(3.27) 

(3.28a) 
(3.28b) 



As we have done with the distribution functions, it will be helpful to introduce the notation 



D('/'\z,z'kl)^^:^Diz.z'k 



2M. 



I), 



D{z,z'klj, 



(3.29a) 
(3.29b) 



for n integer. 



3.2.1 Correlation function in helicity formalism 



Expressing the Dirac matrices of Eq. ( |3.25D in the chiral or Weyl representation, as done in 
Sec. |2.2.1| on page |23|, we get for the leading twist part of the correlation function the expres- 
sion 



(^_A(z,A:^)r- 







ki 2 





D,{z.z^kl] 







0\ 



ieh\^H^{z,z^kl) 



-iQ-'^kLl^H^{z,z^k\) 



D,{z,z^kj-) 










0/ 



(3.30) 



•^The absence of the factor 1/2 in Eq. (3.27) as compared to Eq. (2.35) is due to the absence of an averaging over 
initial states. 
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Figure 3.3. Illustration of the position of the indices of the correlation matrix. 



Restricting ourselves to the subspace of good quark fields and using the chirality basis, we can 
rewrite the correlation function as 



h 



D,{z.z^k] 



(3.31) 



Fig. shows diagrammatic ally the position of the indices of the correlation function. Besides 
being Hermitean, the matrix fulfills the properties of angular momentum conservation and parity 
invariance. Because of the presence of factors e''*^*, we have to take into account / units of angular 
momentum in the initial state, therefore the condition of angular momentum conservation is X2 — 
X2 + i ^rid the condition of parity invariance is 



D{z,kj 



-lyDizX 



Positivity of the correlation matrix implies the bounds 

H^^'){z.z^kl) <D['/'\z.z%). 



(3.32) 



(3.33a) 
(3.33b) 



3.3 Cross section and asymmetries 



In order to obtain the cross section for one-particle inclusive deep inelastic scattering, the expres- 
sions for the spin density matrix of the target, Eq. ( |2.42| ), for the distribution correlation matrix, 
Eq. ( [2.48D , and for the fragmentation correlation matrix Eq. ( |3.31D have to be inserted in the for- 
mula for the cross section 



d^a 



Ip(5) 



a,a; 



XiX[ V dy / 



(3.34) 

Instead of presenting the full cross section, we turn directly to sum and differences of polarized 
cross sections. As in the previous chapter, we will use the symbols to indicate polarization along 
the beam direction and ^ opposite to it. We will also use ] to indicate transverse polarization in the 
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Figure 3.4. Description of the vectors and angles involved in one-particle inclusive deep inelastic scattering. 



direction specified by the angle and [ opposite to it. The subscript U will denote unpolarized 
particles, while L and T will denote longitudinally and transversely polarized particles. The first 
subscript describes always the beam polarization and the second subscript the target polarization. 
Fig. gives a pictorial description of the vectors and angles involved in one-particle inclusive 
deep inelastic scattering. The unpolarized cross section reads [BTp 



d^(T, 



uu 



la' 



dxg dj dz;^ d 05 d^P/^^ sxj^y^ 



My)I,ely\f^{x,,pl)D\{z,,zlkl) 



MM, 



(3.35) 



In the following equations, we will omit to explicitly indicate the variables in which the cross 
section is differential and the variables the distribution and fragmentation functions depend on. 
We define the following polarized cross sections differences 



LL 



d^a 



LT 



d^o- 



d^(7 



d^a, 



UL 



d^CT, 



u- 



dV, 



UT 



^ d^CT, 



d^a, 



(3.36a) 
(3.36b) 



for which we obtain the following expressions in terms of distribution and fragmentation func- 
tions []57|] 



2^2 



d'^LL - E -—T \K\ \S,\ C{y) .y [g\D]\ , 



(3.37a) 



3.3 Cross section and asymmetries 



45 



d^(7, 



q ''■^By 

t/r ^ E TTrf \^t\ { Biy) sin + 0^) ^ 



MM, 



IL 1 



g1 



11 



h 



f±qp.q 
JlT 1 



X J 



Pt PhA 

M 

^{Pt ■ Ph^fikj ■P.J- 2{pj ■ P,J {pj. ■kj.)-pi {kj ■ 



+ B{y) sin {3(j),^-(l)s) 

,2 , 



(3.37b) 
(3.37c) 



2M^M, 



IT I 



(3.37d) 



3.3.1 Transverse momentum measurements 

Due to the presence of an observable transverse momentum, i.e. that of the final hadron, one- 
particle inclusive deep inelastic scattering gives the possibility of extracting some information on 
the transverse momentum of partons. It is convenient to introduce the transverse momentum of 
the hadron with respect to the quark, Kj^ = —zkj. In experiments where jets can be identified, it 
corresponds to the transverse momentum of the detected hadron with respect to the jet axis. The 
simplest quantity to be measured is 



{Pld^^z.)) 



uu 



Jd^.d^P,^ d'a^^ 

Lqelff{x,)D'({z,J 



(3.38) 



where 



(pU^b)) = 



J ^ PTPffti^B^PT 



j d^K^K\D\{z,.Kl) 



(3.39) 



Assuming that all quark flavors have the same transverse momentum distribution we reach the 
result 

.2/„2^„ ^\ ' ^\ (3.40) 



{Pld^B^z,))=zl{pUxJ) + {K\{z 



3.3.2 Transversity measurements 

At the end of the previous chapter, we concluded that in inclusive deep inelastic scattering it is 
not possible to measure the transversity distribution function, h^. In one-particle inclusive deep 
inelastic scattering, we see from Eq. ( p.37dD that it is possible for the transversity distribution to 
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appear in an observable in connection with a chiral-odd fragmentation function, i.e. the Collins 
function. 

It is convenient to introduce the angle ^ = ^j^ + 0^. As specified before, we define the azimuthal 
angles with reference to the electron scattering plane. On the other side, it is possible to choose 
the transverse component of the target spin as the reference axis for the measurement of azimuthal 
angles. When expressing angles with respect to the target spin, we will use a superscript S. The 
relations between the angles in the two different systems is 0^ = — , = ^ ^^^^ 
particular, ^ = ^f^— 2^^. 

The relevant quantity to be measured is the azimuthal single transverse spin asymmetry 

(sind,) (x ^ . ^ l^^s^lhij^^^l^^ 

I^^s^ Ph± <^^uu 

( 1 Ag/) 5(3;) e\ J d^p,^ J \{kj ■ P^JiM^) wn^<A ^^-''^^ 

= 15 ' 



{l/xj,y^)A{y)j:^elf^ixj,)D^^iz, 



In this expression, the transverse momenta of and are entangled in a convolution inte- 
gral [ |151| ]. To simplify the situation, we have to make some assumptions on the transverse mo- 
mentum dependence of the distribution and fragmentation function. 

The simplest example is to suppose there is no intrinsic transverse momentum of the partons 
inside the target [|T5|,|3jl, i.e. 

h,{x,p\)^h,{x)^^. OA!) 
n 

Under this assumption, the pion transverse momentum with respect to the virtual photon is entirely 
due to the fragmentation process, i.e. Pj^^ = Kj, and the convolution can be disentangled 

where the approximation sign reminds us that the equality is assumption dependent. 

Another possibility is to assume that the transverse momentum distribution is Gaussian-like in 
both the distribution and fragmentation side, i.e. 

h^{x,pl)^h,ix)^——, Ht{z,K\)^Hi{z) ,\ • (3.44) 



The convolution becomes [ |141| ] 



(Kliz,)) \P, , I e-'P^x/(^L(-«.^/.)) 



where {P^^} is given by Eq. ( |3.40| ). Using the relation {\aj\)^ = {aj)n/A, valid for Gaussian 
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distributions, we can carry out the integration over Pj^^ and obtain 

(«l> (l''*J> 



/ 



(^r>v^/2 



{\Kt\) 



(3.46) 



where the Xg and Zj^ dependence of the average transverse momenta squared is understood. The 
asymmetry becomes 



{lMA{y)L,elff{x,)D1{z,) 



(3.47) 

where the approximation sign reminds us that the equality is assumption dependent. 

If we want to disentangle the convolution integral of Eq. ( |3.37d[ ) without making any assump- 
tion on the intrinsic transverse momentum distribution, we need to weight the integral with the 
magnitude of the pion transverse momentum [ |127| ]. This procedure results in the azimuthal trans- 
verse spin asymmetry 



sin(^)) (xg,);,zj = 



/ d d^PhL ( \PhL I sin ( d6 (T^ - d6 a, 



UT 



I^^S^^PhA 

{l/x,y^)Biy)z,L,elh1{xs)H;-WHz, 
{l/xj,y^)A{y) L,elff{x,)Dl{z,) 



(3.48) 



We achieved an assumption-free factorization of the Xg dependent transversity distribution and the 
Zfj dependent Collins function. The measurement of this asymmetry requires binning the cross 
section according to the magnitude of the pion transverse momentum. On the other side, this 
asymmetry represents potentially the cleanest method to measure the transversity distribution to- 
gether with the Collins function. Moreover, it turns out that it is possible to study the evolution of 
the moment //j^*^^^ with the energy scale [ |1U4| ], without incurring complications due to Sudakov 
factors [|54l]. However, the inclusion of transverse momentum raises delicate issues related to color 
gauge invariance, factorization and evolution [ |133[ , |124[ |. 

The Collins function can be measured also in e+e annihilation into two hadrons belonging 
to two different jets [|52|,|5T1]. We will only briefly mention this issue at the end of Chap. |^. The 
relevance of this measurement is clear, since an independent measurement of the Collins function 
would be extremely useful to pin down the transversity distribution, despite the problems with 
relating measurements at different energy scales [5^ . 
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3.4 Summary 

In this chapter we studied one-particle inclusive deep inelastic scattering at leading twist, includ- 
ing transverse momenta and neglecting the polarization of the final state hadron. To describe the 
process q* hY we introduced the fragmentation correlation function of Eq. ( p.l2D . At leading 



order in \/Q, the correlation function can be decomposed in two terms, containing the unpolar- 
ized fragmentation function Dj, and the chiral-odd, T-odd Collins fragmentation function, 
[Eq. dCT])]. 



Similarly to what was done in Chap. ^ we cast the correlation function in the form of a forward 
decay matrix in the quark chirality space [Eq. ( [3.31[ )]. Since the matrix is positive semidefinite, we 
were able to suggest for the first time a bound on the Collins function [Eq. ( [3.33bD ]. 



The Collins function is extremely important for possible transversity measurements. Being 
chiral odd, it appears in the single transverse spin asymmetry of Eq. ( |3.37dD in a convolution 



with the transversity distribution. The convolution can be disentangled by measuring the weighted 
azimuthal asymmetry of Eq. ( [3.48[ ). This is therefore a first way to observe the transversity distri- 



bution. 

This method to access transversity poses some problems. First of all, at present we have no 
convincing information on the magnitude of the Collins function. This is a problem common to all 
polarized fragmentation functions, and in particular to T-odd functions. We postpone this discus- 
sion to Chap. where we will mention the scarce information we have on the Collins function, 
we will attempt to estimate it in a consistent model, and we will address the question if T-odd 
fragmentation functions could offer a good chance to tackle the transversity. 

The second problem presented by the Collins function is the need of including transverse mo- 
mentum in the analysis. The study of transverse momentum dependent distribution and fragmenta- 
tion functions is per se a very intriguing subject, but one has to deal with theoretical subtleties (e.g. 
Sudakov factors, transverse gauge links, evolution equations), phenomenological difficulties (e.g. 
describing the transverse momentum dependence of the functions) and experimental challenges 
(e.g. measuring a weighted azimuthal asymmetry). 

In the next chapter, we will examine two-particle inclusive DIS and we will show that the 
transversity distribution can appear in connection with two different fragmentation functions, even 
if the cross section is integrated over the transverse momentum of the outgoing hadron. 
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Two-hadron 
fragmentation functions 

We often think that when we have com- 
pleted our study of one, we know all about 

two, because "two " is "one and one ". We 
forget that we have still to make a study of 

"and". 

Sir A. Eddington 



To observe transversity, an alternative process is represented by two-particle inclusive deep inelas- 
tic scattering, where we can introduce two-hadron fragmentation functions. The transverse spin 
of the target can be correlated via a transversely polarized quark to the relative transverse momen- 
tum of the hadronic pair instead of the transverse momentum of the outgoing hadron, as in the 
case of the Collins function. This provides a way in which the transversity can be probed without 
including partonic transverse momenta, thus avoiding several complications and subtleties. 

Two-hadron fragmentation functions and their relevance for transversity measurements have 
been partially analyzed in some articles [^,[TT7|,[TT^]. The most complete treatment of these func- 
tions has been carried out by Bianconi et al. in Ref. 0. Model calculations have been performed 
in Refs. ^^,[15^. Two-pion fragmentation functions have been studied also in the context of 
e^e^ annihilation with a somewhat different formalism [|27|]. 

In an apparently independent context, semi-inclusive production of spin-one hadrons (e.g. p, 
CO, ^) has been also studied and proposed as an alternative method to measure the transversity 
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distribution [92, 122, 19, 36]. However, to measure the polarization of the outgoing vector meson 
(e.g. p^) it is necessary to measure the four-momenta of the decay products (e.g. n^n^). Thus, 
the reaction ep ep^^{p^ n^n^) is just apart of the more general reaction ep en^n^ J^, 
namely the part where the total invariant mass of the pion pair is equal to the p mass. 

In this chapter, we will go along the same route presented in Ref. ^ and we will complement 
that treatment with the study of the partial wave expansion of two-hadron fragmentation functions. 
This step will prove to be essential to unravel the connection with spin-one fragmentation functions 
and interference fragmentation functions. 



4.1 Two-particle inclusive deep inelastic scattering 

In two-particle inclusive deep inelastic scattering, two of the hadrons belonging to the current frag- 
mentation region are detected in coincidence with the scattered electron, as shown schematically 
in Fig. 



Some experimental results on two pion production are already available [^, ^ . More 
recent data are available only on exclusive production of two pions [ ]57| , |5S| , |59| , [7], ^ and on two- 
hadron production in e+e annihilation [0,0,||,|5|,|T]]. 

If we consider only low invariant masses, hadron pairs are produced mainly in the s-wave 
channel, with a typically smooth distribution over the invariant mass, or in the p-wave channel, via 
a spin-one resonance, with its typical Breit-Wigner invariant mass distribution. This is the case of 
the production of two pions (which can proceed through a p (770) resonance), two kaons (0 (1020) 
resonance), a pion and a kaon (A'* (892) resonance). Fig. ^!2| shows the typical invariant-mass 
spectrum of pion pairs in two-particle inclusive DIS. 

The situation is different for proton-pion production, where on top of the smooth 5-wave con- 
tinuum the sharp spin-half A(1115) resonance appears. Since the A undergoes a parity-violating 
weak decay the formalism we will develop in this chapter will prove to be inadequate. Moreover, 
the resonance is so sharp that it is more appropriate and easier to study it in the framework of 
polarized single hadron production. For these reasons, although polarized A production represents 
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Figure 4.2. Typical invariant mass spectrum ofpion pairs produced in semi-inclusive DIS The p peak 
is visible. 



a very nice way to observe the transversity distribution, we will not touch this subject, for which it 
is possible to consult a vast literature [l83,||,[l|,|T9|,|0l||,|ll[l^,|l7l 0,161]. 



4.1.1 Kinematics 



We assume to detect two hadrons with masses and M2 and momenta and We denote the 
center-of-mass momentum as and the total invariant mass of the system as M^. We maintain the 
parametrization of the vectors P, Pj^ and q as given in Eqs. ( p.3[ ) and Eqs. ( p.lOD , but we introduce 
also the semi-difference of the two hadrons' momenta, R. We parametrize the new momenta 
according toQ 



1 



2 



2 



r {Ml-Ml)-^Ml 



— P 



2P 



h 



{l + Q {Ml-Ml) + ^^Ml 
2 ' IP- 

h 

{l-Q {Ml-Ml) + ^^Ml 
2 ' 



IP, 



h 



(4.1a) 
(4.1b) 
(4.1c) 



The total invariant mass squared, M|, has to be small compared to Q^. For later convenience, we 
compute the scalar quantities 



R' 



{M\+Ml) Ml 



(4.2a) 



'The connection with the variable t, used in Ref. Q is ^ = j (1 + 
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P,.R 



iL^Kll^Mf - (1 - QMf - (1 + QMi 



{Ml -Ml) 



Mi 



k ~\~ I ICj-' I 



2z 



kR = 



{Ml-Ml)-^Ml^^^k^ + \k^ 



Iz 



(4.2b) 

(4.2c) 
(4.2d) 

(4.2e) 



4.1.2 The hadronic tensor 

For a process with two outgoing hadrons we can define the cross section 



2^1 2^2 d'^f? 



which can be rewritten in terms of the different variables as 



; L^y (/ , l',Xe) 2MW^^ {q, P, S, P, , P^) , 



(4.3) 



d9(T 



a' 



-L^y{l,l',Xe) IMW^"" {q,P,S,P^,P^). (4.4) 



d C dMl d 0^ dz;^ '^ ^hi. ^^B d 05 322^^5^:52' 

The angle 0^ is the azimuthal angle of the vector Rj with respect to the lepton plane, measured 
in a plane perpendicular to the direction of the outgoing hadron. Neglecting l/Q corrections, it 
can be identified with the azimuthal angle measured in a plane perpendicular to the virtual photon 
direction or to the lepton beam. Notice the extra factor 8 in the denominator compared to the 
one-particle inclusive case. To obtain the previous formula, we made use of the steps 



d^Pj d^P^ 



d^P^ d^R 



dz, d'P,^ d C d'R^ _ dz, d^P, d C dMl d 



(4.5) 



The hadronic tensor for two-particle inclusive scattering is defined as 

IMW^^ia PS P. P.) = ^ y I 

X {Inf 5^^\q + P- P^^„ P^ P2) H^" {P, S, P^„ , Pi , P2) , 
H^\P,S,P^,,,P,,P2)^{P,S\J^{0)\P,,P2,:^''){P,,P2,^Vm\P,S). 



(4.6) 
(4.7) 



With this definition, it is possible to recover the hadronic tensor for one-particle inclusive scattering 
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by integration over the second hadron 



(4.8) 



where the state represents the sum of the states plus all possible states of the second 
hadron, and P^, = Pa^n -\-P2- 

In analogy to what we presented in Sec. [3.1.2| on page|3^, at leading order in 1/2 the hadronic 
tensor can be expressed in terms of correlation functions as 



2MW^''{q,P,S,P,,P2) = 32z, ^ Tr (4>(xg,P7.,5) A(z„*r, CM^, (^)^) /) 



(4.9) 



where the fragmentation correlation function has been generalized to include also the dependence 
on the vector R 

1 



A{zX,Ml<l>R^k^)^^J dk^A{k,P,,R] 



k-=p-/z 



(4.10) 



The fully differential cross section might be a little bit too complex for experimental measure- 
ments. To simplify the situation, we can perform the integration over the transverse part of the 
center-of-mass momentum, Pj^^. The integrated cross section is 



a' 



-L^y{l,l',K) 2MW^^ {q,P,S,P,,P2), (4.11) 



where 



2Miy^^(^,/',5,Pi,/'2) =32z,Tr[4>(x5,5) A{z,,.C.Mlh) /]• 



A(z,1;,m1^r)^^J dk+ d%Aik,P^,R) 



k-=P,/z 



(4.12a) 
(4.12b) 



To identify the leading twist part of the hadronic tensor and write down the cross section, we 
can trivially repeat the steps described in Sec. p.l.3| on page 
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4.2 The correlation function A 



The most general expansion of the quark-quark correlation matrix A, respecting Hermiticity and 
parity invariance, is [Hp ^ 



r c r r ('4 11') 



M} 

h 



The amplitudes C- are dimensionless real scalar functions Q = C-{k-Pi^,k^,R^,k-R). The last four 
terms are T-odd. 

We are going to consider first the case when no hadron transverse momentum, Pj^^, is de- 
tected and postpone the complete case to the second part of the chapter. We insert the general 
decomposition of Eq. ( |4.13| ) into Eq. ( |3.17c| ) and extract the leading-twist projection 

^_A(z,C,m2,(^,J7" = (^Di(z,C,Mf)+///<(z,C,M2)fc:^ (4.14) 

The prefactor has been chosen to have a better connection with one-hadron results, i.e. integrated 
over ^, M| and 0^. We introduced the two-hadron parton fragmentation functions 



D,{z. CM) = y d% dk^ d{2k.P,J 5 (k\ + ^ - ^ 



Hf{z, CMl) = \ J d% dk^ d{2k.P,) dU + ^+k^-^ 



-c, + -c. 



(4.15a) 
(4.15b) 



The fragmentation function is chiral-odd and T-odd. Like the Collins function, it can be used 
as a partner for the transversity distribution, as we will see in the next section. Notice that this 
function does not require the presence of partonic transverse momentum. Because of this, its 
evolution equations could be simpler than the ones of the Collins function. Since this function has 
the same operator structure as the transversity, it has been suggested that it could have the same 



evolution equations [ ^|160| , |56| ]. However, the situation could be complicated by the presence of 



the dependence on the variables ^ and Mr, which are not present in single particle functions [|161[]. 



^Note that we always use M^^ to render the amplitudes dimensionless, in contrast to what is done in Ref. 
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4.2.1 Correlation function in helicity formalism 



Once again, expressing the Dirac matrices of Eq. ( |4.14| ) in the chiral or Weyl representation, we 
obtain 



1 

ki %n 



/O 









-ie-'^«^//f(z,C,M^) 



0\ 

ie'^«^//<(z,C,M2) 

h 

DiizX^Ml) 




0/ 
(4.16) 

As already observed before, we can restrict ourselves to the subspace of good quark fields and 
adopt the chirality basis, and rewrite the correlation function as 



1 

87r 



ie^^«^//<(,,C,M2^^ 



(4.17) 



From the positivity of the previous matrix, we can derive bounds for the two-hadron fragmentation 
functions defined above: 



Di(z,C,m2)>0, 



\R 



T\ \lT< 



H<{z,CM)\<D,{zXM). 



(4.18a) 
(4.18b) 



4.3 Cross section and asymmetries 

In Fig. ^3] on the following page, we give a pictorial description of the vectors and angles involved 
in one-particle inclusive deep inelastic scattering. The differential cross section integrated over the 
center-of-mass transverse momentum reads 



(4.19) 

Inserting into the previous equation the formulae obtained for the distribution correlation matrix, 
Eq. ( [Z.48[ ), the elementary cross section, Eq. ( p.21D , and the two-hadron fragmentation matrix. 
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Figure 4.3. Description of the vectors and angles involved in two-particle inclusive deep inelastic scattering. 



Eq. ( |4.17| ), we obtain the following result^ 
dJo la} 



d C, dMf^ d dzf^ dx^ dy d Ansx^y^ 



A{y)f^{x,)D\{zXM] 



+ A.5^C(j)g^(x5)D^(z,C,M,2)+5(3;)|S^iy^sin((|,^ + <|,,)/.^(x5)//<^(z,C,M,2) 



(4.20) 



4.3.1 Transversity measurements 

To access transversity in two-particle inclusive deep inelastic scattering it is required to measure 
the azimuthal single transverse spin asymmetry [|150| 1 



/ d <^)5 d (^)^ d C dMl sin((^)^ + <p^) d^a^ 



UT 



jd^,d^i,di;dMld'a, 



uu 



Rn 



{l/xsy^)B{y) L,elh\{x^) j di; dMl^H<Hz,,^.Ml) 



(4.21) 

The most valuable characteristic of this asymmetry is that it does not require the measurement of 
the center of mass transverse momentum, so that the complications connected to the inclusion of 



^The different prefactor appearing in Ref. 150 is due to the different definition of the hadronic tensor and of the 
fragmentation functions. 
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partonic transverse momenta can be avoided. Apart from the usual variables Xg, y, z^, the only 
other variable to be measured is the angle + 0^. In case the transverse spin direction is used as 
a reference axis, instead of the lepton scattering plane, then the angle to be measured is 0^ — l^f . 



4.4 Partial wave expansion 

Up to now, we did not make any study of the inner structure of the two-hadron fragmentation 
functions. It is useful to expand them in partial waves, because if we restrict ourselves to systems 
with low invariant masses, the dominant contributions come only from the lowest harmonics, i.e. 
s and p waves. 



4.4.1 Center-of-mass parameters 



The partial-wave expansion can be performed only in the frame of reference of the center of mass 
of the hadron pair. As a first step, we need to express all vectors in this frame, i.e. 



h 



— — 00 



(4.22a) 



M] + |/?|2 - + |/?|2 - 2|/?| COS e 
2^2 



M2 + |/?|2- /m2 + |/?|2 + 2|/?|cos0 



(4.22b) 



2^2 



pM 
1 



pM 

2 



M2 + |/?|2- |/?|cos0 Jm2 + |/?|2 + |/?|cos0 



v/2 



/?| sinOcos^^, sin0 sin^^ 

i?| sinOcos^g, |/?| sin0 sin^^ 



Ml + |i?|2 - \R\ cos d J Ml + |/?|2 + \R\ cos 



V2 



V2 



(4.22c) 

|/?| sinQcos^^, — sin0 sin0^ 

(4.22d) 



where 



\R\ = 1, /m2 - 2 (m2 +m2) + (m2 -Mlf / Ml. 



(4.23) 



The polar angle 9 is illustrated in Fig. on the next page. The variable C, is connected to the 
center-of-mass variable cos in the following way: 



^^^ = W- (V^i + v/^2 + I^P - 21^1 cos 

h h 



(4.24) 
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'-on, 



p direction 
in the target rest-frame 





If 


e 


P2I 





Figure 4.4. Description of the polar angle 9, in the center-of-mass frame of the hadron couple. 



If M^—M2, notice that R is purely spacelike and 



\R\ = '\R\ = \sjMl-ml 



C = — -JmI-4M} cos 0. 



(4.25) 



4.4.2 The correlation function A with partial wave-expansion 

To perform the partial-wave expansion of the fragmentation functions, let us first express the cor- 
relation function A in terms of center-of-mass parameters 

^_A(z,C(cos0),M2,0^)7- = i-(^Di(z,C(cos0),M2)-Fi//<(z,C(cos0),M2) sinoM^^^^ 

(4.26) 

where n^j = [O, 0, cos 0^, sin 0^1 . In general, we can expand the fragmentation functions on the 
basis of the Legendre polynomials in the following way 

Z)i(z,C(cos0),m2) =£Di„(z,m2)P„(cos0) (4.27) 

n 

with 

^3 = d cos e P,, (cos ) (z, C (cos ) , Mf) . (4.28) 

Considering only two-hadron systems with a low invariant mass, we assume we can truncate 
the expansion in order to include only 5- wave and p-wave contributions to the correlation function. 
The connection between the correlation function A in terms of the variable C, and in terms of the 
variable cos 6 is 



2\R\ 

A(z,cos0,m2,<^^) = J-Ja(z,^,m2,<^^) 



(4.29) 



to take into account the fact that d^ = 2\R\/M, dcos0. Then we can write the partial-wave ex- 
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panded correlation function as 



^_A(z,cos0,Mf,0^)7-^i^ 



Di^uuizM^i'r D\P^^ (z, Mf) cos + D[ (z, Aff) - (3 cos^ - 1 ) 
+ i (h<^p^ (z, Ml) + H<IP (z, Mf) cos e) sin 1^ ^ ^_ , 

(4.30) 

The correlation function can be consequently written in matrix form as 

D(z,cos0,M^,(^^) ^ 

/C2X2 

1 + D\P^^ cos + D/'P ^ 1(3 cos2 - 1 ) +i (h<'J^ + //< cos ) sin |^ e'^« ^ 

h 

^ -i(^i^^^^+^i7^os0)sin0Me-<^« +D^j^^^cos (3cos2 - l) ^ 

(4.31) 

Here, it is not yet clear what are the motivations to assign such names to the functions. It will 
become clear after we distinguish the s- and p-wave contributions in the rest of the section and 
also after we study spin-one fragmentation functions in Chap. ^. 

The information encoded in the correlation function can be expressed in a different way. De- 
pending on the angular momentum of the system, the angular distribution of the two hadrons is 
characterized by a specific bilinear combination of spherical harmonics. If the angular momentum 
can be only or 1, as in the case of our truncated partial wave expansion, the angular distribution 
of the two hadrons can be fully described by the decay matrix 



1 



( 



\7m \7m * 



1 

An 



\ 





--y|sin0e'^« 


V^cos 


y^sin0e-'^« 


-y^sin0e"'^« 


1 sin^ 


— ^cos0sin0e^''^^ 


-|sin2 0e-2i^'« 


a/3 cos 


— 7= cos sin e"^« 


3cos2 


^cos0sin0e"'^« 


y|sin0e''^« 


-|sin2 0e2i'^'fi 


^ cos sin e^*^* 


\ sin^ 



(4.32) 



The angular momentum indices, j and /, run only from to 1. For convenience, we split the 
matrix into blocks: the upper- left block (j = / = 0) refers to the pure 5- wave component, the 
lower-right block (j = / = 1) describes the p-wave component, and it is subdivided according to 
the value of m,m' = +1,0, —1, Finally, the off-diagonal blocks describe the sp interference. 

The original correlation function can now be expressed as the trace of the decay matrix and a 
fragmentation matrix in the quark chirality space ® the hadronic system angular momentum space 



D(z,cos0,m2,<|,^) =D(z,m2)/-''^-^(0,<|,^) 

^2^2 ^2^2 



(4.33) 
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The solution for the fragmentation matrix is the 8x8 matrix 



I A, , . B, , . 

\ fm'.jm fm'jm/ 



(4.34) 



where the inner blocks, spanning the hadronic angular momentum space, read explicitly 



A, , . 

j'm,jm 



B ., , ■ 



( 












—D'^P 

l,UL 












DPP 

l,UU 




-DPP 












— D'^P 

V3 LUL 































DPP 

l,UU 




-DPP 



(4.35a) 



V 












;2v/2 |/f| u<i'P 
\/3 M^"l,UT 


= 2^/2 \R\ TT<sp 





•2^/2 t7<pp 
3 M^"l,LT 














3 M,,^l,Lr 















(4.35b) 



The choice of the indices in the names of the fragmentation functions is connected to their 
position in the matrix. The ss functions are typical of unpolarized two-hadron fragmentation. The 
sp sector describes the interference between s- and ;?-wave fragmentation; this is the sector studied 
by Jaffe et al. in Refs. |1 17| and |163| . The pp functions correspond to the spin-one fragmentation 
functions studied in Refs. |122| , p^|3^ . We will take a look at them from a different point of view in 
the next chapter and also the choice of the subscripts will become more transparent. 

Note that the matrix fulfills the properties of Hermiticity, conservation of angular momentum 
(requiring m + X2 = f^' + Xj) ^^id parity invariance [ |1 17[ ] 



X'2'^2 



—X2 ^^2 



(4.36) 



The imaginary entries of the matrix correspond to T-odd functions. 

In principle the fragmentation matrix could contain more functions, without violating any sym- 
metry, but after tracing it with the decay matrix, they would vanish. In this sense, only part of the 
full information contained in the fragmentation matrix can be analyzed through a parity-conserving 
process [|64|]. We shall come back to this issue after we studied spin-one fragmentation functions 
in the next chapter. 

Finally, if we trace the fragmentation matrix with the decay matrix, we obtain exactly the 
correlation function expanded in partial waves, Eq. ( |4.31D , except for the fact that the unpolarized 
fragmentation function turns out to be expressed in terms of pure s- and p-wave contributions, i.e 



1 



(4.37) 
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In any cross section, these two contributions are merged together and they are kinematically in- 
distinguishable, unless a specific behavior of the invariant mass dependence can be assumed. For 
instance, the p-wave contribution could be due to the existence of a resonance, emerging over a 
background of continuum 5-wave states [|79|] . 



4.4.3 Positivity bounds on partial-wave fragmentation functions 



The fragmentation matrix, Eq. ( [4.34[ ), has to be positive definite, thus allowing us to set positivity 
bounds on the two-hadron fragmentation functions. From the positivity of the diagonal matrix 
elements, we obtain the bounds 



(4.38a) 
(4.38b) 

(4.38c) 



while from the positivity of the two-dimensional minors we get 



\R\ 



D'P 

l,UL 



fj<sp 

\,UT 



IJ<PP 
l.LT 



,3 / 2 \ 3 

< , _ I npp J- -DPP < -D, 



< 



< __ / ( £)PP +-DPP ] { DPP --DPP ] <-D, 



(4.39a) 
(4.39b) 
(4.39c) 



4.4.4 Cross section and asymmetries with partial-wave expansion 



Thanks to the partial-wave expansion, the cross section at leading twist can be expressed in a 
factorized way [ ]112| , |117p 



d^CJ 



dcos 6 dM^ d dz^ dx^ dy d (^^ 
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The unpolarized and polarized parts of the cross section are Q 



uu 



d^CT, 



LL 



UT 



X 



Duuu i^hM) +coseD\P^^{z,M) + \ (3cos2 - l) DfJJz^M^ 



(4.41a) 



9 

X 



D,,uu i^hM) +COS 0D-;^^Jz„m2) + 1 (3cos2 - l) D^JJz^Mf) 



(4.41b) 



(X^e^ \R\ 



9 

X 



(4.41c) 



From the partial wave analysis we see in particular that the transversity distribution can be matched 
with two different chiral-odd, T-odd fragmentation functions, one pertaining to the interference 
between the s- and p-wave channels of two-hadron production, the second being a purely p- 
wave effect. This considerations agree with the past literature on interference fragmentation func- 
tions Ql 17| , |163| ] and spin-one fragmentation functions [ |122| , [T^ ,P^. A priori, we don't know what 
is the magnitude of these functions, nor their behavior with respect to the invariant mass of the 



system. Jaffe et al. Ql 17[ ] studied one of the possible mechanisms that could generate an sp inter- 
ference fragmentation function in production. They separated the production of the pion pair 
(which they did not evaluate) from a nn rescattering process, which determines the T-odd char- 
acter of the fragmentation function and implies a peculiar behavior with respect to the invariant 
mass, shown in Fig. (a). A different model was applied by Radici et al. [ ]150| ]: the expected 



invariant mass behavior turns out to be very different, as shown in Fig. (b), and the magnitude 
of the effect is estimated to be of the order of 1%. For what concerns pp fragmentation functions, 
at the moment there are no estimates of their magnitude and behavior. However, in general ;?-wave 
production of two hadrons becomes significant only when they come from the decay of a spin- 
one resonance. Because of this, we can expect that the invariant mass shape of these functions 
corresponds to a Breit-Wigner curve peaked at the resonance mass. 

If we integrate over cos 0, with the polar angle ranging from to n, the only surviving contri- 
butions to the cross sections are 



^'^uu=llir-\My)M^B)D,,uui^,M) 



Ttsx^y^ 



XeS^C{y)g,{xs)D,ij^{z,M) 



(4.42a) 
(4.42b) 



The distribution and fragmentation functions are understood to have a flavor index q. 
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(b) 



Figure 4.5. Two different models of the invariant mass behavior of the sp interference fragmentation func- 



tion: (a) Jaffe et al. [117], (b) Radici et al. [ 150]. 



UT 



Biy) I y sin(<^)^ + 0,) h,{xs)H<lf^ (z„ . (4.42c) 



It is possible to integrate over cos 6 with the polar angle going from —k/2 to n/2. While the first 
two cross sections vanish, the transverse cross section would turn out to be 



h 



(4.43) 

While the asymmetry in Eq. ( |4.42c[ ) is almost identical to the one discussed in Ref. |117| , the last 
asymmetry contains a contribution of pure p waves that, as already mentioned, will become par- 
ticularly relevant when the two hadrons are produced through a spin-one resonance. 



4.5 The correlation function A with transverse momentum 



We turn now to the correlation functions with the dependence on the transverse momenta, Eqs. ( [2.57p 
and ( |4.1Q ). Consequently, their general decomposition is richer than before. We already know what 
is the form of the correlation function 4> [Eq. ( |2.58[ )]. For what concerns the correlation function 
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A, the inclusion of transverse momentum results in 



+ mt{zXMlkl^j-Rj)^ + Gi{z,i;Mlk\^j-Rr) 



f^^^rs)^- (4.44) 



Once again, we can express the correlation function as a matrix in the chirality space of the quark 



I X .^^y I 



Ml 

n 



1 



Mi 

n 



(4.45) 



We can easily obtain positivity bounds on the transverse-momentum dependent functions 



/! 



<D 



1' 



Ml 

h 



Mt 

h 



2 

1 ; ' 



(4.46) 
(4.47) 



where it is understood that all the functions are dependent on the variables z, ^, M|, k\,kj-Rj. 

Note that if we integrate over kj, the fragmentation function will disappear, while parts of 
the functions Hf and will merge into a single function, Hf, and we will recover the results of 
Sec. ff^on page 15^. 



4.6 Cross section and asymmetries with transverse momentum 

Prior to integration over the center-of-mass transverse momentum we have the nine-fold cross 
section 



dCdM2d(^)^dz,d2p,^ dxgdyd^, 

Ep(^) , ^ 



XiX[ ^ dy / X2X2 



(4.48) 



To simplify the notation, we introduce the projection Uj /\bj = a-£'Jbj. Inserting the formulae 
obtained for the distribution correlation matrix, Eq. ( [2.48[ ), the elementary cross section, Eq. ( [3.21[ ), 
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and the two-hadron fragmentation matrix, Eq. ( |4.45| ), we obtain the following resu 



(X^e^ { \R I 



M ^ 1 



\R I 
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-Biy) cos{2(j)^)^ 
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'A 
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+ 5(y) 



+ 5(3;) sin(2(|),)^ 
+ B{y) cos(2(^J^ 



^(Pr • Phi) (^T ' ^hi) Pt ' j± 

MM,^ ^ 

Akj.) + 



Pt ■ Phi) (PhA. A fc^) + (^r • Ph±)iPh± A Pt) , ± 



{AA9) 



(4.50) 



(4.51) 



(4.52) 



^The UU and f/T cross sections correspond to the one calculated in Ref. 15C except for a difference in the overall 
coefficient - an extra factor 2 in the denominator, due to the use of the variable C, instead of ^, the lack of a factor 
[Inf in the denominator, due to the different definition of the hadronic tensor, an extra factor 2/7r due to the different 
definitions of the fragmentation functions- the use of Mf^ in the denominators instead of Mj, or Mj +M2, and a 
factor 2 difference in the definition of the coefficient C{y). 
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(4.53) 
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(4.54) 



4.7 Partial-wave expansion with transverse momentum 



To proceed with the partial wave expansion of transverse-momentum dependent functions, it turns 
out to be convenient to rewrite the correlation function in a somewhat different way, i.e. 



^1 + 



h 



M. ^ ^ ' 1 



h 



(4.55) 



The functions Hf and have been merged into the complex function 

Hi< (z,C,M2,4,e^(^^-'^^)) =Hi (z,C,m2,4,cos(0,-<^>^)) 

+ e'('^«-^^)//<(z,C,M2,4,cos((^),-(^)^) 



(4.56) 



Because of the presence of the azimuthal angles, this function can be expanded in spherical har- 
monics, provided we retain only those terms which are consistent with the definition in Eq. ( [4-.56[ ). 
If we wish to include only the s- and ;?-wave contributions, the partial wave expansion of the 
function H^^ takes the form 



= H^uy + H^^P^ cos e + Hfll ^ (3 cos2 - 1 ) + e^^-^R^ sin [h^^p^ + H^pp cos 
+ sin B {h^^;/^ + H^pp cos 6 + ^^ + H h<pp' cos 6 

+ ^m-^R) sin2 e H^PP + e-^'(^k-h) sin2 (^H^pp + M //< 



(4.57) 
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Note that the functions on the right-hand side depend only on the variables z, Mf^, kj. 

The expansion of the other two functions, Dj and G^, is considerably simpler, because they 
can depend only on cos(^y^ — 0^) 

^1 = ^wu + ^TuL ^ + ^i.LL \ (3 cos^ - 1) + cos(0;t - 0^) sin (ofjj'^ + Df PP cos 
+ cos(20^-20^) sin^eD^^^, (4.58a) 
= Gf'J^ + GfPPcos +cos((^^ - 0^) sine G^^^, (4.58b) 

where the functions on the right-hand side depend only on the variables z, Mf^, kj. 

As a conclusive remark, note that it would have been possible to expand the functions Hf and 
in a way similar to what we have done for and G^, but some care is required to properly 

treat the component H'^pp which is shared by the two functions. This is why we preferred to 

take an alternative way. In any case, for completeness we give also the expansions of and 
separately: 



^ H^^j^j +H^^P^cos +H^[[ - (Scos^ - l) + 2cos(0^ - 0^) sin [H^lf^ cos 



R\ 



+ 2cos(2(^^ - 2^j^) sin^ H^pp - sin^ jr^H^^P 



<pp 

TT' 



=Hf'J^ + HfPP'cos0+lcos{^^-(^n) sine//; 
As before, the functions on the right-hand side depend only on the variables z, M^, kj. 



(4.59a) 
(4.59b) 



As we have done in Sec. [4-.4.2| on page |5^, it is possible to rewrite the correlation function as a 
trace between the decay matrix, Eq. ( |4.32| ), and a fragmentation matrix in the quark chirality space 
® the hadronic system angular momentum space 



D{z,co^0,Mh^R^k.j^\ 
Once again, the result is a 8 x 8 matrix 
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(4.60) 
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(4.61) 



The inner blocks span the space of the orbital angular momentum of the hadronic system and read 
/ 
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The block C of the fragmentation matrix can be obtained from the block A by imposing parity 
invariance relations. 

The ss sector of the matrix contains functions analogous to the ones we discussed in the previ- 
ous chapter in the context of single-particle unpolarized fragmentation. The pp sector corresponds 
to the spin-one fragmentation functions studied in Ref. We will discuss them again in the 
next chapter. The sp interference sector has never been studied with the inclusion of transverse 
momentum. 

Similar considerations to the ones discussed after Eqs. ( [4-.35[ ) hold for the transverse momentum 
dependent fragmentation matrix. In particular, the matrix fulfills the properties of Hermiticity, 
conservation of angular momentum {m-\-X2 = ^' -^'Xi'^ and parity invariance 



£)] tn ,jm 

X2X2 



.1^1 DJ' -rn' ,j -m ^ 
^2 -^2 



(4.62) 



The imaginary parts of the matrix represent T-odd functions. Note that the off-diagonal blocks, 
which are chiral-odd, can only contain T-odd functions. The matrix could in principle contain 
other functions, but they are lost when tracing with the decay matrix and thus cannot be analyzed 
by a parity conserving decay process. 

From the fact that the matrix Z)(z,M|,^^) is positive semi-definite, it is possible to obtain 
positivity bounds on the fragmentation functions, as we have done already before. 



4.8 Summary 

In this chapter we analyzed two-particle inclusive deep inelastic scattering, at leading order in 
1/2, without and with partonic transverse momentum. For the description of the fragmentation 
side, we introduced a correlation function A that depends on the center-of-mass momentum of the 
two hadrons and on their relative momentum. 

In the first part of the chapter we neglected partonic transverse momentum. We distinguished 
one chiral-even T-even and one chiral-odd T-odd fragmentation function [Eq. 04-.14D ]. They are 



complicated objects depending on three variables, whereas single-particle fragmentation functions 
depend on one variable only. We derived positivity bounds on these functions [Eqs. ( [4.18D ] and we 
computed the cross section of two-particle inclusive scattering [Eq. ( 14.20[ )]. 

Since two-particle systems with a low invariant mass are usually produced in the s and p waves, 
we performed a partial wave expansion of the fragmentation functions taking into consideration 
only these two lowest modes. We split the chiral-even fragmentation function into three contribu- 
tions, pertaining to the pure s waves, pure p waves and sp interference. We split the chiral-odd 
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function into two contributions, typical of pure p waves and sp interference. These five new frag- 
mentation functions depend only on two variables, namely z and the invariant mass squared of the 
system [Eq. ( [4-.30D ]. Employing the usual helicity formalism, we discussed positivity bounds on 
the new fragmentation functions [Eqs. (|4.38D and ( |4.39| )1. We computed all non vanishing spin 
asymmetries of two-particle inclusive DIS at leading twist [Eqs. ( [4.41D ]. In particular, we showed 
that in the single transverse spin asymmetry of Eq. ( |4.41c[ ) the transversity distribution appears in 
connection with both chiral-odd fragmentation functions. Integrating the asymmetry over different 
ranges of the polar angle 0, we presented two distinct ways to access the transversity distribution. 
The first one [Eq. ( |4.42c| )1 involves only the sp interference term and corresponds to the one dis- 
cussed in Ref. |1 \1\ while the second one [Eq. ( |4.43[ )] involves also the pure p term, which could 
have an entirely different physical origin and should be particularly relevant in the presence of 
a spin-one resonance. These two asymmetries are in all respects two distinct ways to access the 
transversity distribution function. The second one has never been clearly indicated in the literature. 

In the rest of the chapter, we repeated the analysis of the correlation function including partonic 
transverse momentum and also in this case we applied a partial wave expansion. 
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Looks like we've made it, 
Look how far we've come, my baby. 
You are still the one that I love. 
The only one I dream of. 

S. Twain 



In the previous chapter we examined the production of two hadrons in the s and p waves. As 
we already mentioned, the p-wave sector of the analysis should overlap with the description of 
the fragmentation functions for spin-one mesons. The reason is that the polarization of the vector 
meson can be analyzed by measuring its decay in two other particles, as in the case of a p meson 
decaying into two pions. 

In this chapter, we will take a small detour from the mainstream of the thesis, and we shall deal 
with spin-one targets and spin-one fragments. In the first part, we will examine what is necessary 
for the description of a spin-one object and what are the differences with the spin-half case. We 
will introduce distribution functions for quarks inside spin-one objects, with and without transverse 
momentum. 

Then we will turn our attention to spin-one hadrons in the final state and we will examine 
the fragmentation functions for a quark producing a spin-one hadron. In the end, we will check 
whether the connection with two-hadron fragmentation can be established and if there are any 
differences between the two situations. 
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5.1 Spin-one targets 



In recent years some attention has been devoted to distribution functions characterizing spin-one 
targets, starting from the work of Hoodbhoy, Jaffe and Manohar [ |110| ]. Unfortunately, the only 
available spin-one target is the deuteron, which is essentially a weakly bound system of two spin- 
half hadrons. In fact, deuteron targets are often used in deep inelastic scattering with the major 
purpose of extracting the neutron distribution functions. However, the spin-one structure of the 
deuteron in itself can be very interesting [|144|, |5^^|169|]. A large amount of work is already 



present on deuteron structure functions, especially in Drell-Yan processes [106, 108]. 



5.1.1 Spin density matrix and spin tensor 



The description of particles with spin can be attained by using a spin density matrix p in the rest 
frame of the particle. The parametrization of the density matrix for a spin-7 particle is conveniently 
performed with the introduction of irreducible spin tensors up to rank 27. For example, we have 
already seen [cf. Eq. ( [2.42D ] that the density matrix of a spin-half particle can be decomposed on a 
Cartesian basis of 2 x 2 matrices, formed by the identity matrix and the three Pauli matrices. 



1 



p = -(l + S'a') 



(5.1) 



where we introduced the (rank-one) spin vector S'. 

To parametrize the density matrix of a spin-one particle we can choose a Cartesian basis of 
3x3 matrices, formed by the identity matrix, three spin matrices E' (generalization of the Pauli 
matrices to the three-dimensional case) and five extra matrices E'-'. These last ones can be built 
using bilinear combinations of the spin matrices. In three dimensions these combinations are no 
longer dependent on the spin matrices themselves, as would be for the Pauli matrices. We choose 



them to be (see Refs. [64| , |134| and ^ for a comparison) 



E'^' = ^ (E'E^' + E^E') - ^ 1 d'j. 
With these preliminaries, we can write the spin density matrix as 



(5.2) 



(5.3) 



where we introduced the symmetric traceless rank-two spin tensor T'j . We choose the following 
way of parametrizing the spin vector and tensor in the rest frame of the hadron. 
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(5.4a) 



(5.4b) 
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The parameter S^^ is often called alignment [ |5^|134| , pT] , |93| ]. 

Although we introduced the spin vector and tensor in the particle rest frame, it is possible 
to introduce covariant generalizations of them, satisfying the conditions P^S^ = and P^T^^ = 
0, where P^ is the momentum of the hadron. In a coUinear frame, the parametrizations of the 
covariant spin vector and tensor will be 
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(5.5a) 



(5.5b) 



In the particle rest frame (where P^ = Mj they correspond to Eqs. ( 15. 4| ) (note that here they 
are written in light-cone coordinates). 

Inserting in Eq. ( |5T3| ) the spin vector, Eq. ( |5.4aD , and the spin tensor, Eq. ( |5.4b| ), the explicit 
form of the spin density matrix p turns out to be 



3 1^1 
2^2 2V2 
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2V2 



2V2 



3 +'^LL 
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2^/2 
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(5.6) 



The treatment of spin-one particles can be done equivalently by introducing the complex polar- 
ization vector £, and its covariant generalization e^. The two formalisms can be related by means 
of the formulae [|134p 



S = Im (e* X £) 
or for the covariant generalizations 



Re 



£*£ 



Re 



8^v 



p2 



(5.7) 



(5.8) 



5.1.2 Interpretation of the components of the spin tensor 

A particular component of the spin tensor measures a combination of probabilities of finding the 
system in a certain spin state (defined in the particle rest frame). As "analyzing" spin states we can 
choose the eigenstates of the spin vector operator in a particular direction. We can write the spin 
vector operator in terms of polar and azimuthal angles. 



E'n ■ = cos ■& cos 9 -f cos '&sin^-\- sin i^, 



(5.9) 
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and we can denote its eigenstates as |m, 9)), m being their magnetic quantum number. The 
probability of finding one of these states can be calculated as 



^ [m, (p)] = Tr {p |m, (t^, (p))(m, (t?, (p)| } . 
From the explicit formula of the density matrix, Eq. (|5.6|), one can compute 



(5.10) 



^LL = -I <^ [1 ' (0. 0)] - I <^ [- 1 ' (0' 0)] + f [0, (0, 0)] , (5- 1 la) 

Sir = ^ [0, (-f ,0)] - ^ [0, (f ,0)] , = ^ [0, (-f , f )] - ^ [0, (f , f )] , (5.11b) 

5« = ^ [0, (f , -f )] - ^ [0, (f , f )] , 5^^ = ^ [0, (f , f )] - ^ [0, (f ,0)] . (5.1 Ic) 

Below, we suggest a diagrammatic interpretation of these probability combinations. Arrows repre- 
sent spin states m = +1 and m = — I'm the direction of the arrow itself, while dashed lines denote 
spin state m = in the direction of the line itself. 

transverse plane 



Q - a 



The probabilistic interpretations suggest straightforward bounds on the values the spin tensor 
parameters can achieve, namely 



1 2 



1 <5^y < 1, 



-l<^ir<l, 

where /, j = x,y. Finally, it is possible to define a total degree of polarization 

Is 3 

^[1[SI + (Sf )2 + (5^^)2] + ^ [35^ + (SIt)' + (Sir)' + (^¥1? + (S^r)'] ' ''^ 
whose value ranges between and 1. 



(5.12) 



(5.13) 



5.2 The correlation function <I> 
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5.2 The correlation function O 

The correlation function ^» has to fulfill the condition of Hermiticity and parity invariance. For 
spin-one hadrons, they are slightly different from the spin-half case, because of the presence of the 
spin tensor 

Hermiticity: ^{p,P,S) = f^^p,P,S,T)f, (5.14a) 

parity: ^{p,P,S) = )P^{p,P,-S,f))P (5.14b) 

where p^ = d'^^Pfx and so forth for the other vectors, and = 5^^ d^^Tfyp. 

The most general decomposition of the correlation function 4> for spin-one hadrons isQ 

4>(p, P, 5, r) = MAi 1 + A2 ^ + A3 + (^^ a^yP^p''^ + (iAg p • 5 75) 



M 



(5.15) 



+ ^P,PvT'^' i+^P,P.T^'H^P,P.T^'l^ 
+ (^^P.PvT^'Op.pPp-^ +A,,p^T^^yy+{^aypPPp^T^' 

The amplitudes A- are real functions A ■ = A-(p • P, p^) . The terms with the amplitudes A4, Ag, A^j, 
Ajg, Ajg, A^g and Ajg (included between parentheses) constitute the T-odd part of the correlation 
function, according to the definition 

^lo^^ipM = -irV ^T-odd(i^'^'^'3')irV- (5.i6) 

The leading-twist part of the correlation function, parametrized in terms of five distribution 
functions, is [P^,P7|] 

^+ Mx, T) 7+ = (/i {x) {x) 5^ 75 + ^1 (^) Jsh +b,{x) S^^ + ih^{x) ^^^) (5.17) 
The underlined function, /zj^^, is T-odd. Apart from the usual distribution functions defined in 



^In Ref. 



107 



a similar decomposition of the correlation function was attempted, but in an incorrect way. 
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Eq. ( [2.34] ), we introduced the spin-one parton distribution functions 

b^{x)= j d^p-^dp^d{2p-P)5{pj-+x^M^ + p^-2xp-P) 



4M2 



Mx ' 



(A +xA ) p' + Pt-M'x^ ^ 
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h^i^j-{x) = J d^pjdp^d{2p-P) 5 {p\+x^hf + p^-2xp-P) 



1 

AM 



Mx 



Mx) {-A,^-xA,^) 



(5.18a) 



(5.18b) 



Note that the distribution function b^, introduced in Ref. |1 10| , was called in a different way in 
Refs. ^ and ^ to follow a more systematic naming of the functions, especially in view of the 
inclusion of transverse momentum, as we will do in Sec. |5.2.2| . Finally, note that the function 
h^^j is T-odd, but does not require the presence of intrinsic transverse momentum, a feature that 
is absent in T-odd distribution functions for spin-half targets at leading order in l/Q. It could 
be interesting to study this object in order to clarify the role of intrinsic transverse momentum in 
generating T-odd effects in distribution functions [ [7^|124] ]. 



5.2.1 Correlation function in the helicity formalism 



The only difference from the analysis of Sec. [2.2. 1| on page |2^ is that we have to deal with a more 
complex spin density matrix and with a 3 x 3 target spin space. The connection with the correlation 
function and its matrix representation is 



where 
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Eventually, the 6 x 6 leading-twist scattering matrix turns out to be Q 
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(5.21) 



where the functions on the right-hand side depend on the variable x only. The matrix F{x) is 
Hermitean, parity invariant and conserves angular momentum (Aj -\-x[ = Aj + x^). 
From the positivity of the diagonal elements we obtain the bounds: 



/iW>o 

\g,{x)\<f,{x)-^-b,{x) < ^/iW, 
while positivity of two-dimensional minors gives the bound 

{h^{x)f+{h,^^{x)f < i (^f^{x) + ^b,ix)^ (^f^(x)+g,{x)-^-b^{x)^ 

<^^/iW+^iW)'<^(/iW)', 



(5.22a) 
(5.22b) 

(5.22c) 



(5.23) 



This bound is a generalization of the Soffer bound 01590 and must be fulfilled by any spin-one 
target. If we assume that T-odd distribution functions vanish due to time-reversal invariance, then 
the bound can be reduced to 



1/^1 W I < y ^ (/i W + ^ ^1 {x?j (/i {x)+g,ix)-^-b, (x)) 



(5.24) 



^Note the difference of sign in the imaginary components of the matrix with respect to Ref. ^ 
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5.2.2 Inclusion of transverse momentum 

We now consider the correlation function unintegrated over the parton transverse momentum, as 
defined in Eq. ( p.57D . For convenience, the correlation function can be decomposed in several terms 
in relation to the polarization state of the target, i.e. 4> = 4>jy + + + 4>^^ + 4>^y + ^tt- 
leading order in 1/2, these terms can be decomposed as (T-odd terms are underlined) 



^+4>j;(x,p^,5,r)7+= |/i+i^^|^+, (5.25a) 
^+ (x, , 5, r) 7+ = 5^ 1^ 1^ Tg + /zf^ Tg 1 1 , (5 .25b) 



(5.25c) 



^+4>^Jx,p^,5,r)7+ = 5^^ (5.25d) 

^+4>^^(x,p^,5, r)7+ = {/i,^ +^erVp^ 75 

+i^^,,+i^^|}^^, 



(5.25e) 



^ ix D ST)v+-'>f Pt-^tt-Pt pan PtPtp 



(5.25f) 



All the distribution functions on the right-hand side are understood to depend on x and p^. 
Similarly to the function /zj, defined in Eq. ( |2.60| ), it is convenient to define the functions 

KlA^.Pt) = h'lLTi^^pl) +hf^^\x,pl), (5.26a) 

h,TA^,p^T)=h[TA^.PT)+hfj^i^.PT)- (5.26b) 

In the rest of this Section, unless otherwise specified, all the functions are understood to depend 
on the variables x and pj. 

As we have already done in Eq. ( [2.62[ ), to simplify the notation it is preferable to introduce 
new complex functions with a real part, corresponding to a T-even function, and an imaginary part, 
corresponding to a T-odd function. To avoid the introduction of new names, we will simply call 
the new functions with the name of their real part, in the following way 

/?j -t-i/Zj^^ — > /Zj, f\TT~^ '^S\^TT^f\TTi (5-27) 



The correlation function ■ 
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Following steps analogous to the previous section, we can reconstruct the complete 6x6 scat- 
tering matrix. 



F{x,pj) 



A\A, 



B 



a;a, a;a, 
I C 

\ A'lA; A'jA;^ 



(5.28) 



The inner blocks span the target spin space and they are 



a;a, 
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(5.29) 



B 



a;a, 



^e>^.[/.^,-i(/.^-i/.^J] 



2 m2 



Mr 



y2/z* 

(5.30) 



\ 



This matrix is Hermitean. Conservation of angular momentum is guaranteed (A'^ +x[+l' = -^1 + 
X^)- The matrix is also parity invariant, i.e. 



F(..,p^)^'A = (-l)'V(x,/,^)-^i-^; 



(5.31) 



In fact, block C of the fragmentation matrix can be obtained from block A by imposing parity 
invariance relations. 

Since the matrix F{x,Pj) is positive semidefinite, we can extract bounds on the distribution 
functions. If we fully exploit the positivity of the scattering matrix, we can write several relations 
involving an increasing number of different functions. We feel this to be an excessive task if 
compared to the exiguity of information we have on the functions involved. Therefore, here we 
choose to focus only on the relations stemming from positivity of the two-dimensional minors of 
the matrix. 

Because of the symmetry properties of the matrix, only nine independent inequality relations 
between the different functions are produced :0 



\PtI 
2M2 



l^ir+/iLrl'< (/i + ^) (/i+^i-y 



(5.32a) 
(5.32b) 



Note the difference in Eq. (5.32g) compared to Eq. (38) of Ref B7L due to an error in the latter. 



80 



5. Spin one 



\Pt\' 
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4M2 
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< (/i+^i-^) (/i-^i-y 



3 



4M4 



l/irrl'<(/i+^i-y) (/i-^i-y) 



(5.32c) 
(5.32d) 
(5.32e) 
(5.32f) 
(5.32g) 
(5.32h) 
(5.32i) 



5.3 Cross sections and asymmetries for spin-one targets 

We will not dwell very much on the analysis of spin-one targets. The only available one is the 
deuteron, which is a system of two spin-half hadrons connected by nuclear interaction. In the 
approximation of independent scattering off the nucleons, the distribution function is expected 
to be small. However, calculations of b^ are available in the literature (see e.g. Ref. |169| ). Some 
calculations show that this function could be sizeable in the low-jc region [ |144| , |62| , p9|], due to 
nuclear shadowing effects. New precision measurements of the deuteron structure functions (e.g. 
from HERMES) may provide an experimental test of this expectation, as already observed in 
several papers pT| , |131| , p31 ]. 

We will focus only on totally inclusive deep inelastic scattering off a spin-one target at leading 
twist. As we have seen in Chap. ||, in this type of process we can access only chiral-even, transverse 
momentum independent functions. Spin-one targets are characterized by three such functions, /j, 
and b^. Therefore, we should be able to define three independent asymmetries. In addition to 
the indices ^ to denote the ± 1 polarization along the direction of the beam, with the index we 
will denote when the target is polarized along the direction of the beam but with magnetic number 
equal to 0. 

The total cross section integrated over the azimuthal angle 0^ is analogous to Eq. ( [2.7(JD 

^ I^J [^(3^) if^i^s) +^eS,B{y)gl{xs)] • (5.33) 



A few words on the value of the polarization coefficients are needed. When the target is polarized 
purely in the ±1 state along the beam direction, then 5^ = =f1 and = — 1/3 [cf. Eq. ( |5.1 la| )]. 
In case the polarization is not complete, then a possible notation could be 5^ = =f|5^| and 5^^ = 
— \Sii\/?>, with \Si\ = \Sii} < 1. When the target is polarized purely in the state along the beam. 



5.4 Spin-one hadrons in the final state 
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then 5*^ = and S^^ = 2/3. In case the polarization is not complete, we can use the notation 
5^^ = 2 1 5^^1/3, where the coefficient \S^j^\ is smaller than 1, and can be different from l^l^l. If we 
just sum the cross section with opposite vector polarization S^, in analogy to Eq. ( |2.71D , we will 
obtain 

i (d^CT^ + d^CT^) = ^^A{y) £4 (ffix,) - ^-^bfix,)) , (5.34) 

where a contamination of tensor polarization is still present. Thus, the unpolarized part of the cross 
section has to be properly defined as 

1 Atioc^ 
= :t {d^<y^^ + d^O^, + d^a^^) = -—2A{y) ^e^/f (x^). (5.35) 



The vector polarized part of the cross section corresponds exactly to Eq. ( [2.72p . We rewrite it here 
for convenience 

= ^ (d^^^-.^ - d^o^^) = --2 \S,\Biy) l^els^ix,). (5.36) 
Finally, we can define a tensor polarized part of the cross section 

d\, = - (2 d^a^, - d^a^^ - d^a^^) = \SlM{y) Y^e\b\{x,). (5.37) 

L D SX^y q 

Therefore, the appropriate asymmetry to measure the spin-one distribution function is 



5.4 Spin-one hadrons in the final state 

Instead of pointing the attention to spin-one targets, it is possible to analyze spin-one final state 
hadrons in semi-inclusive deep inelastic scattering or in e+ annihilation. This idea was first 
considered by Efremov and Teryaev [ |^ ] (and more recently rediscussed in Ref. [ |155| ]). A system- 



atic study was accomplished by Ji 0122| ], who singled out two new fragmentation functions, the 
function Sj, analogous to the distribution function b^, and the T-odd function ti-. These new frag- 
mentation functions can be observed in the production of vector mesons, e.g. p,K* , ^. However, 
these functions require polarimetry on the final-state meson, which can be done by studying the 
angular distribution of its decay products (e.g. in the case of meson). In this sense, vector 

meson production represents just a specific contribution to the more general case of two-particle 
production near the vector meson mass. 
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5.4.1 The decay of a spin-one hadron 



We want to describe the decay of a spin-one hadron with momentum and mass into two un- 
polarized particles with momenta P^^, Pj, and masses Mj, Mj. We assume the same parametrization 
of the momenta as presented in Chap. ^ in particular Eqs. ([4.1D. In general, the angular distribution 
of the decay products of a spin-one hadron into two unpolarized hadrons is 



W'(COS0,<^)^)=P^2^2^(0,0^ 



lAjJA^ 



(5.39) 



where G and 0^ are the polar and azimuthal angles of the vector R = {P^ — Pj)/! in the decaying 
particle rest frame, as defined in Chap. ^ Eqs. ( [4. 22] ), and they correspond to the angles usually 
measured in experiments [|69|,^[7l]. The decay matrix, has been defined in Eq. ( |4.32| ). In this 
chapter, we are only interested in the j = 1 sector, therefore for convenience we define a spin-one 
decay matrix 

i?(0,<^>^) ^^(0,<^)^) (5.40) 
As can be checked by explicit comparison, this matrix can be rewritten as 



1 

47t 



(5.41) 



where ^cm = ^cm/ l-^l • Notice that \R\ stands for the modulus of the vector R in the center-of-mass 
frame of reference. We know from Eq. ( t4.23| ) or Eq. ( [4-.25[ ) how to express it in terms of M^, M2 



and Mj^. 



In general, the decay matrix can be expressed in terms of analyzing powers: 



(5.42) 



and the decay distribution can be obtained accordingly as 

By comparing Eq. ( |5.41| ) with Eq. ( |5.42D we can identify 

A' = 0, 



(5.43) 



^ ~ 3 ^cm^cm' 



(5.44) 



from which we observe that the vector analyzing powers are zero [|5^].f*| 
The tensor analyzing power can be written in a covariant form: 

1 1 / P^P"^ 

i/?|2 3r M? 



(5.45) 



The vector analyzing powers can be different from zero in a decay that does not conserve parity. 
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We can use a parametrization of the tensor analyzing power analogous to that of the spin tensor, 
Eq. ( |5.5bD , provided we exchange the plus and minus components and we replace P and M with 



and Mj^. Then, we can write the parameters of the tensor in terms of the angles 9 and 0^: 

A^^ = i (Bcos^e-l) , (5.46a) 

= -sin20 cos^^, A^^ = -sin20 sin^^, (5.46b) 

Afj = - sin^ G cos2<^^, A^^, = - sin^ G sin2(^^. (5.46c) 



Substituting the explicit form of the tensor analyzing power in Eq. ( |5.43| ), we obtain the decay 
distribution (cf. Ref. p36| ) 

3/2 

W{cosG,^j^) = ^ -+S^j^{3cos^G-l)-SlTsm26 cos(j)j^-Sl^sm26 sin0^ 

^ . (5.47) 

- Sfj sin^ G cos 2(^^ - SJ^ sin^ G sin 20^ j . 



5.5 The correlation function A 

The decomposition of the correlation function A is analogous to that of <I>, just replacing p, P, 
A-{p ■ P,p^) with k, Pf^, Bj{k ■ Pj^.k^) and the spin vector and tensor, and T^^ , with the vector 
and tensor analyzing powers, A^ and A^^ . Therefore, the leading-twist part of the correlation 
functions, A(z, A) 7^, can be parametrized in terms of five fragmentation functions 

^_A(z, cos 0,(^)^)7- = ^ [D,{z) + G,{z)A^y^+H,{z)y,Jl.T+B,{z)A^^ + iH,^^T{z)ALT) ^-^ 

(5.48) 

The prefactor has been chosen to have a better connection to the unpolarized results, i.e. integrated 
over cos G and 0^. We retained the terms containing the vector analyzing powers, although we 
know that they vanish in the case of a parity-conserving decay. The function H^jjj- is T-odd and 



corresponds to the function hj of Ref. |122| . This function evolves with the energy scale in the same 



way as the transversity fragmentation function [ |160[ ] . 



The matrix representation of the leading part of A(z) can be obtained almost directly from the 
result of the matrix F{x), shown in Eq. ( |5.21| ), provided we apply the correct replacements. The 
result is 

D{zff2 = 
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D,+G, 


-\B. 








V2{H^+iH^^j.) 
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\ 



(5.49) 

Conservation of angular momentum in the fragmentation process requires A2 + = + ;f2- 

At this point, we can compare the matrix we obtained with the pure p-wave sector of the 
matrix D{z,M^) in Eq. ( |4.34D . We notice that the two matrix have an almost identical form, after 
we identify 

3 



B, 



H 



ILT 



DPP 

^ l,UU' 

3 

-DPP 

I [A\^H<PP 
4V 3M, 1'^^ 



(5.50a) 
(5.50b) 
(5.50c) 



There are, however, a couple of differences. The first difference is that two-hadron fragmentation 
functions depend also on the total invariant mass squared, while spin-one fragmentation functions 
do not, since they are supposed to be nonzero only at the mass of the hadron we are observing. 
More appropriately, since we are usually observing a resonance we should assume that spin-one 
fragmentation functions have to be multiplied by an invariant mass distribution, e.g. a Breit-Wigner 
curve peaked at the resonance mass. 

The second difference is that in the analysis of two-hadron fragmentation we missed the com- 
ponents Gj and H^, typical of vector polarization. The reason is that a parity violating decay is 
required to analyze these components. In fact, if we trace the matrix of Eq. ( |5.49D with a decay 
matrix such as the one defined in Eq. ( |5.40| ), the vector polarization functions, Gj and H^, would 
not appear in the decay distribution. When studying two-hadron fragmentation, we imposed parity 
invariance from the very beginning - from Eq. ( [4.1 3D - thus missing any contribution of vector 
polarization. 

We are not going to analyze the features of the correlation function A when partonic transverse 
momentum is included, as they essentially reproduce those of the correlation function studied 



in Sec. 5.2.2 on page 78, modulo the replacements we discussed at the beginning of this section. 



5.5.1 Positivity bounds on spin-one fragmentation functions 



The bounds stemming from the positivity of the fragmentation matrix are very similar to the ones 
we discussed for spin-one distribution functions in Sec. |5.2.1| on page From the diagonal 



5.6 Cross section and asymmetries for spin-one production 
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elements we get 



-^-D,{z)<B,{z) < 3D,{z), 
|Gi(z)|<Di(z)-i5i(z) < ^D,{z). 



(5.51a) 
(5.51b) 

(5.51c) 



Eqs. ( |5.51aD and ( |5.51bD correspond to the Eqs. ( |4.38b[ ) and ( |4.38cD . From positivity of the two- 
dimensional minors we get 

{H,LTiz)f+{H,{z)f<\(^D,{z) + lB,{z)^ (D,{z) + G,{z)-lB,{z)y (5.52) 

Because of the lack of information on vector polarized fragmentation functions, the bound in 
Eq. ( |5.51cD does not have a direct relevance for experiments, but it can be very useful to test 
the consistency of model calculations. The bound in Eq. ( p.52\ ) can be used in a less restrictive 



version [37] 



\HiLTiz)\<J{D,iz) + ^B,iz)^ (^D,{z)-lB,iz)^<^D,iz), (5.53) 



corresponding to Eq. ( [4.39cD . 



5.6 Cross section and asymmetries for spin-one production 

We consider one-particle inclusive DIS events where the target consists of a spin-half hadron and 
the fragment is a spin-one hadron. We assume the polarization of the final state hadron is ana- 
lyzed by means of a two-particle, parity-conserving decay. It is worthwhile to mention that the 
experimental analysis of spin-one final state polarization has been already performed for exclusive 
leptoproduction of vector mesons [|57|,§B|,§^[7|,||] and in e+e" annihilation [0,0,||,|,|T1]. Fig. ^ 
on the next page shows a typical distribution of the energy absorbed by the target in p production 
events. The peak represents the exclusive events, while the rest of the distribution contains deep 
inelastic events, which have always been excluded from data analyses. 

Let us first analyze the cross section integrated over the transverse momentum of the outgoing 
hadron. It will be differential in six variables 



d^CT 



dcos d(^)^ dz, dxg dyd(j)s ^^"^-"ajA; 



dy 



-) 



(5.54) 



where and are the decay angles discussed in Sec. |5.4.1| on page |S2. 
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Figure 5.1. Distribution of the energy absorbed by the target in p production events /|^. 



The unpolarized and polarized parts of the cross section are 

1 



My) /f (^fi) D\{z,) + - (3cos2 e - 1) Bliz, 



■h> ] ' 



d'^LL = I ^^^2 S.Ciy) glix,) (D\iz,) + \{3 cos2 - 1 ) Bliz,: 



UT 



Insxgy 



-B{y) \S^\ sinie sin(0^ + 05)/i^(xg)//j^^^(z,). 



(5.55a) 
(5.55b) 
(5.55c) 



The unpolarized cross section corresponds to Eq. (58) of Ref. ^ (there the unpolarized fragmenta- 
tion function has been neglected), while Eq. ( |5.55b| ) corresponds to Eq. (61) of Ref. ^ integrated 
over the hadron transverse momentum and finally Eq. ( |5.55c| ) corresponds to Eq. (63) of Ref. P^p| 
At the same time, note the correspondence with the pp sector of Eqs. ( |4.41| ), after we apply the 
identifications of Eqs. ( |5.50[ ). 

The asymmetry of Eq. ( |5.55cD contains the transversity distribution multiplied with the chiral- 
odd T-odd function H^^j. This observable is just a part of the one we presented in Eq. ( |4.41cD , but 
this is a good point to stress once more a couple of details. The function H^^^ is different from 
the sp interference fragmentation function H'^^j- Although a proper analysis of this function has 
to be done in the framework of two-particle fragmentation, it is important to realize that it enjoys 
the characteristics of a single-particle function. For instance, there is little doubt about the fact 
that it will be strongly peaked at the mass of some spin-one resonance. On the contrary, we don't 
know what is the invariant mass behavior oisp interference fragmentation functions. Moreover, its 
physical origin could be more similar to that of the Collins function than to that of sp interference 
fragmentation functions. Since we have some indications that the Collins function is possibly 
sizable, there is hope that also H^j^j, will turn out to be large, maybe larger than sp interference 
fragmentation functions. However, it must be remembered that H^^^j, does not depend on partonic 
transverse momentum. 



^ An overall minus sign is missing in Eq. (63) of Ref. 
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If we don't integrate over the transverse momentum of the outgoing hadron, the the cross 
section will be differential in eight variables 



dcos d(j)j^ dzj, d dxg dy d (j)g 



XiX'i ^ dy / X2X2 



Ri^^<l>R)j,^j,r (5-56) 



Obviously, the cross sections are much more complex than before. The following formulae cor- 
respond to the ones listed in Ref. ^ if we replace the tensor polarization components with the 
analyzing powers of Eqs. ( |5.46D . However, here we show also the unpolarized contributions and 
we examine also the terms with T-odd distribution functions 



Unpolarized lepton beam and unpolarized target 



Aiy) \ Jf [/i + i (3cos2 Q-\)J {f, 5 J 



^ Insxj^y'^ 
— sin 20 cos{^i^ — 0^) y 
2{kj ■Pf^^'^ — kj 



k P 



ILT 



sin^ 9 cos(2^y^ — 2^^) 



X 



X 



X 



M} 

n 



TT 



OL^e^ I 1 

+ E ^7;;:7^ 5 (y) <| - - (3 cos^ e - 1 ) cos 2(^), 



2{pj ■ Pj^j){kj^ ■ Pj^^) Pjkj 



MM,, 



+ sin20 cos(0;^ + 0^) 



Pt PhA 

M 



hiH 



1 ^^ILT 



+ sin 6 cos 2<^^ 



MM, ^11^^ 



+ sin20 cos(3(^^ — (^^) 



X y 



4(^r • PhA.fiPT -Phi)- ^i^T ■ Phi) (Pt -K)- kliPr ' Ka 

2MM} 

n 



■hiH 



ILT 



+ sin^ cos(40,^ - 2^^) y 

K^T-PhLfiPT-Kl 



kj - ^{kr ■ Ph±f [Pt ■ kr " ^(kj " Ph±) (Pt " Ka 



IMMl 

h 



IMMl 

n 



"1 ^irr 



(5.57) 



^Moreover, in Ref. ^ there are a few typos and there is a systematic error sign whenever the factors and 



1 5^^^ I are used. 
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Polarized lepton beam and unpolarized target 



^^^LU = - L C{y) I sin20 sin(0, - 0^) ^ 



fiG 



ILT 



+ sin^ 6 sin(20^ - 20^) 



h 



AG 



ITT 



(5.58) 



Unpolarized lepton beam and longitudinally polarized target 



+ sin^ e sin(20;j - 20^) 

-E 



5l|AW.|sm2esin((j.,-,/i„)^ 

2(^r '^hl)^ ~^T 



k p 



8ilG 



ITT 



X 



X 



2{pj ■ Pf^j){kj ■ Pf^^) Pt'^t IT 1. 



PT-PhJ 

M 



MM,, 



hl^lLT 



+ sin20 sin(0^ + 0^) 



+ sin^0 sin20^J^ 



~mm7 



+ sin20 sin(30;j — 0^) 



^{kT-PhLfiPT-PhJ-^i^T-PhJiPT-^T) -kliPr-PhA 

2MM2 

h 



'h\LH\LT 



+ sin^ e sin(40^ - 20^) 
8(^r-^±)'(/'r-^±) 



4 - 4(A:r " PhA.f [Pt ' ^t ^ "^{^t ' Ph±)(PT ' Ph±)] 



IMMl 

h 



2MMl 

h 



h\LH\TT 



(5.59) 



Polarized lepton beam and longitudinally polarized target 



1 



— sin20 cos(0^ — 0^) 

— sin^ Q cos (20^ — 20^) ^ 



k P 



Mu 



ILT 



^{^T'Ph±)'^ ~^T 



(5.60) 



Ml 

h 



ITT 
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Unpolarized lepton beam and transversely polarized target 



d^%r = E 2^^^ I ^(3^) j sin20 sin(0^ - 0^) ^ 



2MM. 



ILJ 



■ sin 20 sin(20^ ~^r~ ^s) ^ 

■ sin^ 6 sin(0^ — 2<^^ + 0^) 

■ sin^ 6 sin(3<^^ - 2</)^ - </)^) ^ 



2MM 



ILT 



h 



IMMl 

h 



X sin(0;j - 0^) 



SlT^lTT 

Pt Ph± 



IMMl 8lT^lTT 
h 



IMMl 

h 



fvrD, 
M •'^^ 1 

sin20 sin(0^ — 0^) 



- sin 20 sin(20^ ~^r~ 0^) -'^ 
+ sin^ e sin(0/j - 20^ + (j)^) J 



■ sin^ sin(30^ - 20^ - 0^) J 



M -^i^^i 



1 
3 



+ ^(3cos2 0-l) 



2MM 



2MM 



h 



2{kj ■ P^^i^Pj -kj) l^{pj-P^ 



2MMl 

h 



ITT 



UjkT ■ PhA_?{PT ■Ph±)- 2(fcr • Ph±)iPT ■ fer) 



2MM2 

h 



kj{Pj-Pf^_ 

2MMl 

h 



cxP'c^ I I 

+ E2^«W(3(3-s^e-l)-nW» + « 

sin20 sin(0^ + 0^) [h^ ^ilj] + sin^ 

— sin20 sin(20^ ~ 0/? + 0^) 



X 



xsin(0^- 20^-0^)^ 



(5.61) 



' ^ h-\- 

^^^^ h^H^j.j 



X 
X 
X 



2{kj -Py^ 



\2 t2 



2M2 

h 



— sin^ sin(30^ — 20^ + 0^) 



'4ik^-P,J'^-3klik,^-P, 



2Ml 

h 



+ ^(3cos2 0-l)sin(30;j-05) 



4(/>r ■ PhA.f{kT ■Ph±)- ^(Pt ■ Pha) (Pt " ^r) " PUh " Ph±) ,± „± ' 
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— sin 20 sin(20^ + ~ ^s) 

- sin^ Q sin((^^ + 2(^^ - (/)^) 
+ sin 20 sin(40^ -^r- 0^) 



AM^Ml 

n 



X 



AM^Ml 

n 



h\Tti\LT 



+ sin^ e sin(50;j - 20^ - 0^) 



'24(fc^-P,^)[2(p^./>,^)2-p2] 



4M2m3 

h 



+ 



X 



h 
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Polarized lepton beam and transversely polarized target 

Pt PhA 



M 



+ -(3cos2 0-l) 



M 



— sin 20 cos(2(^y^ ~^r~ ^s) 

— sin^ cos(^^ - 2^^ + y 

— sin^ 9 cos(3^y^ — 2^^ — 0^) 



sin20cos(0^ — 



(Pt ' ^T) „ r) 

2MM, ^ir^i^r 



2MM, 



ILT 



h 



2{kj ■Ph^){PT '^t) ~^^{Pt 'PjiA 

IMMl 

h 



ITT 



Hk^ ■ PhA.f{PT -PhD- 2(^r • Ph±) (Pt ■ ^t) 
2MM} 

h 



2MMl 

h 



SitD 



ITT 



sin 20 COS (0^ — ^g) J' 



~2MM7 



IitG 



ILT 



+ sin 20 cos(2^y^ ~^r~ ^5) 
- sin^ cos(^^ - 2(j)j^ + J 
+ sin^ cos(30^ - 20^ - J 



"^^P T P hL>^T ' P hL> Pj kj 



2MM, 



\LT 



h 



■flT^ 



IT ^ ITT 



2{kj -Pfij^jPT '^t) ~k^{PT 'PjiA 

2MMl 

h 

" '^ {kj ■ Pi,^f{PT -PhJ- ^(kT ■ Phi) (Pt ■ ^t) 
2MMl 

h 



krjPT'PhA 

2MMl 

11 



flT^i 



TT 



(5.62) 



The previous formulae describe the pp sector of two-hadron production with transverse mo- 
mentum and partial- wave expansion, which we did not compute in Chap. H. 



5.7 Summary 

In this chapter, we examined the description of spin-one hadrons, both in the role of targets and of 
final-state fragments. 

To describe spin-one targets, we introduced a rank- two spin tensor in the decomposition of the 
spin density matrix. The presence of the spin tensor made the decomposition of the correlation 
function 4> more complex than the spin-half case. Neglecting partonic transverse momentum, 
we could describe the leading twist part of the correlation function by means of five distribution 



functions [Eq. ( |5.17[ )]. Besides the familiar three functions of the spin-half case (/j, and h^), we 



needed to introduce the function and the T-odd function h^jjj.. The latter was never discussed 
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in the literature before and it is the first leading-twist T-odd distribution function that survives the 
integration over partonic transverse momentum. 

We studied the correlation function in the framework of the helicity formalism and we obtained 
for the first time positivity bounds on spin-one distribution functions [Eqs. ( |5.22| ) and ( |5.23| )1. In 
particular, we demonstrated that the bounds on the helicity and transversity distribution functions 
are different from the spin-half case. We carried out the analysis of the correlation function includ- 
ing also transverse momentum, leading to the decomposition of the correlation function shown in 
Eqs. ( |5.25D , together with the positivity bounds of Eqs. ( |5.32[ ). 

We briefly discussed the cross section for inclusive deep inelastic scattering on spin-one targets, 
emphasizing the contribution of the distribution function [Eqs. ( |5.34{ ), ( |5.35D and ( |5.37| )]. 

We devoted the second half of the chapter to the analysis of spin-one hadron as current frag- 
ments. We explained how it is possible to probe the polarization of an outgoing hadron if it un- 
dergoes a two-particle decay. After studying the correlation function A, in Sec. we listed all 
the spin asymmetries occurring at leading order in 1/2 in spin-one deep inelastic leptoproduction, 
with or without integration over the outgoing hadron's transverse momentum. In particular, we 
highlighted the single transverse spin asymmetry of Eq. ( |5.55c[ ). Even though this asymmetry was 
already included in the discussion of two-hadron fragmentation functions, in this chapter we made 
clear that it has the characteristics of a single-particle fragmentation function and it is not related 
to the two-hadron sp interference function H^y'j- 

In the next chapter, we are going to study how it is possible to generate T-odd single-particle 
fragmentation functions in the framework of a field theoretical approach. 
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A model estimate 
of the Collins function 

Models are to be used, not believed. 

H. Teil 



In the previous chapters, we extensively discussed how it is possible to observe the transversity 
distribution function in connection with T-odd chiral-odd fragmentation function, in particular the 
Collins function H^, the two-hadron interference function H^^j or the spin-one fragmentation 
function H^^j. Unfortunately, at the moment we have scarce or no information about these func- 
tions. Therefore, we need to address the important issue of estimating them and check whether they 
could be measurable. Beside giving an indication of the magnitude of unknown fuctions, model 
evaluations, however rough or naive, serve some useful purposes: they show whether a nonzero 
function can be obtained in the framework of known theories, they pave the way for future im- 
proved estimates, they shed light on some crucial properties, and they analyze the consequences 
of the assumptions of the model. In this chapter, we focus on the Collins function, which could be 
regarded as the prototype of a T-odd function, and present an estimate of it. 

In spite of the apparent difficulty in modeling T-odd effects, a nonvanishing Collins function 
can be obtained through a consistent one-loop calculation, in a description where massive con- 
stituent quarks and pions are the only effective degrees of freedom and interact via a simple pseu- 
doscalar coupling, as we have shown in Ref. [B^. We point out that Collins himself suggested the 
idea of dressing the quark propagator as a possible mechanism to produce a nonzero [ [77| ] and a 
more specific way to achieve this goal was mentioned by Suzuky in Ref. |162| . The model discussed 
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in our early work was admittedly unfit to reproduce the phenomenology of the Collins function. 
Here, we show the result of a recent calculation [^, which gives a more reasonable estimate. 
At present, only one attempt to theoretically estimate the Collins function for pions exists [E^, 



and little phenomenological information is available from experiments. The HERMES Collabora- 
tion reported the first measurements of single spin asymmetries in semi-inclusive DIS [[T^,[ll]], giv- 



ing an indication of a possibly nonzero Collins function. Preliminary results have been presented 



by the SMC collaboration as well [|65|]. The Collins function has also been invoked to explain large 
azimuthal asymmetries in pp^ nX [OFDp. However, all these analyses are plagued by large 



uncertainties, and hence do not allow any conclusive statement yet. Recently, a phenomenologi- 
cal estimate of the Collins function has been proposed [PO], combining results from the DELPHI, 



SMC and HERMES experiments. However, in spite of all the efforts to pin down the Collins 
function, the knowledge we have at present is still insufficient. 

In this chapter we calculate the Collins function for pions in a chiral invariant approach at a low 
energy scale, as we have done in Ref. We use the model of Manohar and Georgi [ |139| ], which 



incorporates chiral symmetry and its spontaneous breaking, two important aspects of QCD at low 
energies. The spontaneous breaking of chiral symmetry leads to the existence of (almost massless) 
Goldstone bosons, which are included as effective degrees of freedom in the model. Quarks appear 
as further degrees of freedom as well. However, in contrast with the current quarks of the QCD 
Lagrangian, the model uses massive constituent quarks - a concept that has been proven very 
successful in many phenomenological models at hadronic scales. With the exception of Ref. |121|, 



the implications of a chiral invariant interaction for fragmentation functions into Goldstone bosons 
at low energy scales remain essentially unexplored. 

Although the applicability of the Manohar-Georgi model is restricted to energies below the 
scale of chiral symmetry breaking ^ I GeV, this might be sufficient to calculate soft objects. 
In this kinematical regime, the chiral power counting allows setting up a consistent perturbation 
theory [ |170| ]. The relevant expansion parameter is given by l/A^, where / is a generic external 



momentum of a particle participating in the fragmentation. To guarantee the convergence of the 
perturbation theory, we restrict the maximum virtuality ji^ of the decaying quark to a soft value. 
We mostly consider the case jU^ = 1 GeV^. 

The outline of the chapter is as follows. We first give the details of our model and present the 
analytical results of our calculation. Next, we discuss our results and compare them with known 
observables, indicating the choice of the parameters of our model. Then, we present the features of 
our prediction for the Collins function and its moments. Finally, using the outcome of our model, 
we estimate the leading order asymmetries containing the Collins function in semi-inclusive DIS 
and in e+e^ annihilation into two hadrons. 



6.1 Calculation of the Collins function 

Considering the fragmentation process of a quark into a pion, q*{k) 7t{p)Y, we use the expres 
sions of the unpolarized fragmentation function Z)^ and the Collins function given in Eq. ( [3.28[ ) 



6.1 Calculation of the Collins function 
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For convenience, we reproduce that definition in a more explicit way:[] 

D,{z,z^k'r)= Ak+ l:x[^{k,p)y 
e'-'k 1 r 



(6.1a) 



_ ; 



(6.1b) 



with niy^ denoting the pion mass and A being the correlation function defined in Eq. ( [3.12D . 



We now use the chiral invariant model of Manohar and Georgi [ |139| ] to calculate the matrix 
elements in the correlation function. Neglecting the part that describes free Goldstone bosons, the 
Lagrangian of the model reads 

J^ = f{i^ + Y-m + g^4Y^)xj/. (6.2) 
In Eq. ( |6.2| ) the pion field enters through the vector and axial vector combinations 

Vf'='-[u\d^u], A^ = ^{wt<9M, (6.3) 

with u = exp(iT • k/2Ft^), where the T,- are the generators of the SU(2) flavor group and Fj^ = 
93MeV represents the pion decay constant. In absence of resonances, the pion decay constant 
determines the scale of chiral symmetry breaking via = AkF^. The quark mass m and the axial 
coupling constant are free parameters of the model that are not constrained by chiral symmetry. 
The values of these parameters will be specified in Sec. |6]^. Although we limit ourselves here to 
the SU(2) flavor sector of the model, the extension to strange quarks is straightforward, allowing in 
particular the calculation of kaon fragmentation functions. For convenience we write down explic- 
itly those terms of the interaction part of the Lagrangian (p^ that are relevant for our calculation. 



To be specific we need both the interaction of a single pion with a quark and the two-pion contact 
interaction, which can easily be obtained by expanding the nonlinear representation u in terms of 
the pion field: 

^7iqq = -j^^^yi^y5'^-^^^^^ (6.4a) 

■^K7tqq = -^Wyii'C-{jt^di'jt)\lf. (6.4b) 

Performing the numerical calculation of the Collins function, it turns out that the contact interaction 
( |6.4b[ ), which is a direct consequence of chiral symmetry, plays a dominant role. 

At tree level, the fragmentation of a quark is modeled through the process q* nq, where 
Fig. ^?T]on the next page represents the corresponding unitarity diagram. Using the Lagrangian in 
Eq. ( |6.4a| ), the correlation function at lowest order reads 



^0, M = -§2(^ im-t^'n)* 2;r 5 ((* - p)^ - m^) . (6.5) 



'Note that this definition of slightly differs from the original one given by Collins [|77[]. 
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Figure 6.1. Lowest-order cut diagram describing the fragmentation of a quark into a pion. 



This correlation function allows to compute the unpolarized fragmentation function by means 
of Eq. ( |6.1a| ), leading to 



Note that the expression in Eq. ( |6.6| ) is only weakly dependent on the transverse momentum of 
the quark. In fact, is constant as a function of kj, if iriji = and (or) m = 0. Because our 
approach is limited to the soft regime, we will impose an upper cutoff on the kj integration, as will 



be discussed in more detail in Sec. |6|^. This in turn leads to a finite D^{z) after integration over 
the transverse momentum. 

The SU(2) flavor structure of our approach implies the relations 

j^U-.K+ ^ j^d-^K+ ^ jyU-^K- ^ ^d-^K- ^ , (6.7b) 



where is the result given in Eq. (|6.6|). In the case of unfavored fragmentation processes Dj 
vanishes at tree level, but will be nonzero as soon as one-loop corrections are included. Accord- 
ing to the chiral power counting, one-loop contributions to are suppressed by a factor /^/A^ 
compared to the tree level result. The maximum momentum up to which the chiral perturbation 
expansion converges numerically can be determined only by an explicit calculation of the one-loop 
corrections. 

As in the case of a pseudoscalar quark-pion coupling [Q], the Collins function turns out to 



be zero in Bom approximation. To obtain a nonzero result, we have to resort to the one-loop level. 



In Fig. |63 the corresponding diagrams are shown, where we have displayed only those graphs 



that contribute to the Collins function. The explicit calculation of Hj- is similar to our previous 



work 034] ] . The relevant ingredients of the calculation are the self-energy and the vertex correction 



diagrams. These ingredients are sketched in Fig. ^3] and can be expressed analytically as 

^ ^ 4F2 7 [iTif \[k - Z)2 - m2] [/2 -m\y 

1^ '^^ 8^3 ^^J [iKf [{k-p-lY-m^] [{k-lY-m^][l^-miy ^ ^ ^ 

Tikn^- i r / (/ + i^)(/-^ + m)/ 
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Figure 6.2. One-loop corrections to the fragmentation of a quark into a pion contributing to the Collins 
function. The Hermitean conjugate diagrams (h.c.) are not shown explicitly. 



P 



P ^ 



Figure 6.3. One-loop self-energy, and vertex corrections. 



where flavor factors have been suppressed. For later purpose, we give here the most general 
parametrization of the functions E, Fj and Fj, 

Z{k)=A^ + Bm, (6.9a) 
r,{k,p) = (q +Dj + E,^ + Fjjt^ , (6.9b) 

^lif'^P) = ^^5 (C^+Dj + E^ll: + F^^i) . (6.9c) 

The real parts of the functions A, B, Cj, Dj etc. could be UV divergent and require in principle 
a proper renormalization. Here, we do not need to deal with the question of renormalization at 
all, since only the imaginary parts of the loop diagrams are important when calculating the Collins 
function [|4|]. 

Taking now flavor factors properly into account, the contributions to the correlation function 
generated by the diagrams (a), (b) and (c) in Fig. 15^ are given by 



X J* 75 i±4 ^(*) 2,t S{{k-i,)--m^), (6. 10a) 

1 (^ + m 
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(6.10b) 



(6.10c) 



The correlation functions of the Hermitean conjugate diagrams follow from the Hermiticity condi- 

Summing the contributions of all diagrams and inserting the resulting correlation function in 
Eq. (OH), we eventually obtain the result 



1^ - 3m Im (a + 5) - Im (q - mE^ + {k^ - m^)F^) 



+ 2 Im (Cj - mE^ + {k^ - tn^)F^ ^ 



(6.11) 



Thus, the Collins function is entirely given by the imaginary parts of the coefficients defined in 
Eqs. (|6.9|). We can compute these imaginary parts by applying Cutkosky rules to the self-energy 
and vertex diagrams of Fig. on the page before. Explicit calculation leads to 



2m^ 



k^ — np- / — m 

1 - 



2 



Im (Ci - mE, + ik^ - m^)F,) = m {k^ - n?) 



2 V k^ 

2 2n /^3fc^+"^^-"^l 



'1 ' 



Ik^ 



(6.12a) 
(6.12b) 



, . 2 k^-m^ + ml ( 2 2 o 2w 



,111 ^•••'H 



Im (C2 - mE^ + [k^ - m^)F^) = ^^^jpi ^ 



m') 1 



2 2 
m —m 



k^ 



(6.12c) 



where we have introduced the so-called Kallen function, X{k^ ,np- ,m\) = [k^ — {m -\- m^Y]\k?- — 
(m — m^Y], and the factors 

= j d"^ I d -ml) d {{k-lf - m2) = ^ (yt2, m2, m2 ) {k^ - (m + m^f) , (6.13a) 

2 y ' (fc_^_/)2_^2 



2A/A(F,m2,m2; 



In 



1 + 



fc2m2 — (m2 — m2 )2 



(r-(m + m;j)2) . 



(6.13b) 



These integrals are finite and vanish below the threshold of quark-pion production, where the self- 
energy and vertex diagrams do not possess an imaginary part. 

Thus, Eq. ( |6.1 ID in combination with Eqs. ( |6.12| ) gives the explicit result for the Collins func- 
tion in the Manohar-Georgi model to lowest possible order. Because of its chiral-odd nature, the 
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Collins function would vanish in this model if we set the mass of the quark to zero. The same phe- 
nomenon has been observed in the calculation of a chiral-odd twist-3 fragmentation function [ |121| ]. 
The result in Eq. ( |6.1 1| ) corresponds, e.g., to the fragmentation u — > n^. The expressions for the 
remaining favored transitions are obtained in analogy to Eqs. (|6.7|). Unfavored fragmentation pro- 
cesses in the case of the Collins function appear only at the two-loop level. 



6.2 Estimates and phenomenology 

6.2.1 Unpolarized fragmentation function and the choice of parameters 

We now present our numerical estimates, where all results for the fragmentation functions in this 
subsection refer to the transition u^n^. To begin with we calculate the unpolarized fragmentation 
function (z), which is given by 

D,{z) = nj^ dKlD^{z,Kj-), (6.14) 

where Kj = —zk^ denotes the transverse momentum of the outgoing hadron with respect to the 
quark direction. The upper limit on the Kj integration is set by the cutoff on the fragmenting quark 
virtuality, /i^, and corresponds to 

4m..=zil-z)^^-zm^-il-z)ml. (6.15) 

In addition to m and g^, the cutoff jl^ is the third parameter of our approach that is not fixed a 
priori. However, as will be explained below, the possible values of /i^ can be restricted when 
comparing our results to experimental data. Unless otherwise specified, we always use the values 



m 



0.3GeV, ,?A = 1' Ai^ = lGeV2. (6.16) 



At the relevant places, the dependence of our results on possible variations of these parameters 
will be discussed. A few remarks concerning the choice in Eq. ( |6.16| ) are in order. The value of 
m is a typical mass of a constituent quark. The choice for the axial coupling can be seen as a kind 
of average number of what has been proposed in the literature. For instance, in a simple SU(6) 
spin-flavor model for the proton one finds ^ 0.75 in order to obtain the correct value for the 



axial charge of the nucleon [ ]1390 . On the other hand, large Nc arguments favor a value of the order 
of 1 [ [171^ , while, according to a recent calculation in a relativistic point-form approach [|59[], a 
g^ slightly above 1 seems to be required for describing the axial charge of the nucleon. Finally, 
our choice for /i^ ensures that the momenta of the outgoing pion and quark, in the rest frame of 
the fragmenting quark, remain below values of the order 0.5 GeV. In this region we believe chiral 
perturbation theory to be applicable, meaning that our leading order result can provide a reliable 
estimate. 



In Fig. 5.4 on page 101 we show the result for the unpolarized fragmentation function D'f 



Notice that in general the fragmentation functions vanish outside the kinematical limits, which in 
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our model are given by 

_ 1 

■^max,min 2 




7 ? \ ^ 9 



^2 ; 



(6.17) 



corresponding to the situation when the upper limit of the Kj integration becomes equal to zero. 
We consider our tree level result as a pure valence-type part of D"^^^ . The sea-type (unfavored) 
transition a ^ ;r+ is strictly zero at leading order. Therefore, we compare the model result to the 
valence-type quantity D"^^ — D"^^ , where the fragmentation functions have been taken from 
the parametrization of KretzeiQ [ |129| ] at a scale = I GeV^. Obviously, the z dependence of both 



curves is in nice agreement, which is a nontrivial result. For example, in the pseudoscalar model 
that we used in our previous work []34l], Z)j behaves quite differently and peaks at an intermediate 



z value. 

On the other hand, we underestimate the parametrization of Ref. |129| by about a factor of 2. 
Some remarks are in order at this point. Although a part of the discrepancy might be attributed 
to the uncertainty in the value of g^, the most important point is to address the question as to 
what extent we can compare our estimate with existing parametrizations. The parametrization 
of [ |129| ] serves basically as input function of the perturbative QCD evolution equations, used to 
describe high-energy e^e~ data, and displays the typical logarithmic dependence on the scale Q^. 
A value of 2^ = 1 GeV^ is believed to be already beyond the limit of applicability of perturbative 
QCD calculations. On the other hand, our approach displays, to a first approximation, a linear 
dependence on the cutoff /i^. It is supposed to be valid at low scales and it is also stretched to 
the limit of its applicability for ji^ = \ GeV^. In this context it should also be investigated to 
what extent the inclusion of one-loop corrections, which allow for the additional decay channel 
q* — > TtTtq, will increase the result for at /x^ = 1 GeV^. Finally, we want to remark that to our 
knowledge there exists no strict one-to-one correspondence between the quark virtuality /i^ and 
the scale used in the evolution equation of fragmentation functions, which in semi-inclusive DIS is 
typically identified with the photon virtuality Q^. For all these reasons, a smooth matching of our 
calculation and the parametrization of [ ]1290 cannot necessarily be expected. Despite these caveats, 
the correct z behavior displayed by our result for suggests that the calculation can well be used 
as an input for evolution equations at a low scale. In the next subsection we will elaborate more on 
this point in connection with the Collins function. 

The best indication of the appropriate value of the cutoff jU^ may be obtained when comparing 
our calculation to experimental data of the average transverse momentum of the outgoing hadron 
with respect to the quark, which we evaluate according to 

{\K^\{z)) = f ""'''dKl\K^\D,{z,K^T)- (6.18) 
(zj Jo 



2r 



^Other parametrizations [|125L B3l [1301] use a starting energy scale > 2 GeV^, which is too high to allow a 
comparison with our results. Moreover, the valence parts of these other parametrizations would have to be obtained 
with some extra assumptions and seems too arbitrary. 
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Figure 6.4. Model result for the unpolarized quark fragmentation function D"^'' ( solid line) and compari- 
son with the parametrization ofRef. (gray line). 
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Figure 6.5. Model result for the average hadron transverse momentum for different choices of the cutoff: 
jj,^ =0.5 GeV^ ( dotted line ), jX^ — 1 GeV^ ( solid line j, jU^ = 1 .5 GeV^ f dashed line ) and comparison with a 
fit to experimental results from DELPHI mU (gray line). 



102 



6. A model estimate of the Collins function 




In Fig. ^3] on the page before we show the resuk of this observable as a function of z for three 
different choices of the parameter /i^. As a comparison, we also show a fit (taken from Ref. |2^) to 
experimental data obtained by the DELPHI Collaboration As in the case of D^{z), the shape of 
our result is very similar to the experimental one, which we consider as an encouraging result. For 
^2 = 1 GeV^ our curve is about 30% below the data. Such a disagreement is not surprising, keeping 
in mind that at LEP energies higher order perturbative QCD effects (e.g. gluon bremsstrahlung, 
unfavored fragmentations, etc.) play an important role, leading in general to a broadening of the Kj. 
distribution. For experiments at lower energies, however, where perturbative QCD contributions 
can be neglected in a first approximation, it may be possible to exhaust the experimental value for 
{\Kj\{z)) with genuine soft contributions as described in our model. This in turn would determine 
the appropriate value of the cutoff /i^. For example, such a method of matching our calculation 
with experimental conditions could be applied at HERMES kinematics, even though the method 
is somewhat hampered since Kj is not directly measured in semi-inclusive DIS. In this case, one 
rather observes the transverse momentum of the outgoing hadron with respect to the virtual photon, 
which depends on both Kj and the transverse momentum of the partons inside the target pj, 
as shown in Eq. ( |3.40D . 



6.2.2 Collins function 



We now turn to the description of our model result for the Collins function. In Fig. //^ is 
plotted for three different values of the constituent quark mass, m = 0.2, 0.3, 0.4 GeV. In a large z 
range, the function does not depend strongly on the precise value of the quark mass, if we choose it 
within reasonable limits. That is why we can confidently fix m = 0.3 GeV for our numerical studies. 
It is very interesting to observe that the behavior of the unpolarized fragmentation function is 
quite distinct from that of the Collins function: while the former is decreasing as z increases, the 
latter is growing. The different behavior of the two functions becomes even more evident when 



looking at their ratio, shown in Fig. At present, there exists some evidence of z behavior of the 



Collins function and it is in agreement with our results. We briefly discuss this subject in Sec. ^ 
on page |109| . 
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Figure 6.7. Model result for /D^. 



The ratios of the Collins function or any of its moments with are almost independent of the 
coupling constant gj^. The reason is that the one-loop correction containing the contact interaction 
is only proportional to g\, as is, and is dominating over the others. Furthermore, the ratio 
H^/D^ is nearly independent of the cutoff /i^. In conclusion, the prediction shown in Fig. is 
almost independent of the choice of parameters in our approach. 

At this point we would like to add some general remarks concerning the z behavior of our 
results. It turns out that the shape of all the results does not vary much when changing the pa- 
rameters within reasonable limits. In particular, variations of and of the cutoff jl^ only change 
the normalization of the curves but not their shape. In this sense our calculation of fragmentation 
functions has a good predictive power for the z behavior of the function. This has a direct practical 
consequence if one uses, for instance, our result of the Collins function as input in an evolution 
equation: the z dependence of the input function can be adjusted to the shape of our H^, while its 
normalization can be kept free in order to account for uncertainties in the values of g^ and jU^. 

In Fig. 0| on the next page we plot the ratio 



n 



H^^'l^){z) ^ 
D,{z) ^D;{zj 



Izm.^ 



(6.19) 



which enters the transverse single spin asymmetry of Eq. ( |3.43D . This quantity rises roughly lin- 
early within a large z range, leading to a similar z behavior of the transverse spin asymmetry. 

is no longer independent of the cutoff jU^, but rather the same dependence as in the 
case of {\Kj\) (Fig. on page |101D can be assumed. In Fig. |6]^ on the next page, this ratio is 
compared to the expression 



{\Kr\{z)) Htiz) 



H 



Iznij^ D^{z) 



n 



D\{z) 



dK4 



Iznij^ 



(6.20) 



A very close agreement between the two different curves can be observed, indicating that the 
model predicts a quite similar transverse momentum dependence of both the Collins function and 
Z)j. In the literature, this feature is sometimes assumed in phenomenological parametrizations of 
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Figure 6.8. Model result for H'^^^/'^'^ /Dy (solid line) and comparison with the product 
{{\Kj\) /Izniji) [Hy /Dy) (dashed line). Note that the positivity bound requires the ratio to be smaller than 
0.5. 



H^. Note, however, that in our approach deviations from this simple behavior can be expected, if 
Z)j is also calculated consistently to the one-loop order. 

The Collins function has to fulfill the positivity bound of Eq. ( p.33b[ ), which can be rewritten 



as 



1^. 



Ht{z,K^)<-D,{z,4). 



1 



(6.21) 



(6.22) 



Integration over gives the simplified expression 

D,{z) 

which is satisfied by our model calculation. It is clear, however, that increasing the value of /i^ will 
eventually result in a violation of the positivity condition. To avoid such a violation, we should 
calculate and H^- consistently at the same order, i.e., the one-loop corrections to should be 
included. By doing so, the positivity bound will be fulfilled even at larger values of /i^, for which 
our numerical results are no longer trustworthy. 

From our results, we expect an increasing behavior of the azimuthal asymmetry in p ^ nX 
as function of Xp, qualitatively similar to what has been predicted in Ref. ^in the context of the 
Lund fragmentation model. At this point, it is also interesting to discuss the comparison of our 
results with the ones obtained using the so-called "Collins guess". On the basis of very general 
assumptions, Collins suggested a possible behavior for the transverse spin asymmetry containing 
[|77|]. This suggestion has been used in the literature (see, e.g., Refs. |145| , |128| , ^|137D to 
propose the following shape for the Collins function: 



H^^'/^\z)^nJdKj 



M. 



c 



■D,{z,4) 



(6.23) 



2z M^ + Kj/z^- 

with the parameter ranging between 0.3 and O.VGeV.^ Using our model outcome for the 
unpolarized fragmentation function, we apply Eq. ( |6.23| ) to estimate H^'^^/'^\ and in Fig. we 

^Note that even this particular form does not correspond precisely to what proposed in Ref. even if it is often 
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show how this compares to Eq. ( |6.19[ ). There is a rough qualitative agreement with the Collins 



ansatz for the lowest value of the parameter M^, although it is not growing fast enough compared 
to Eq. ( |6.19| ). On the other hand, in the Manohar-Georgi model there is no agreement with the 
Collins ansatz for high values of the parameter M^, which might indicate that the relation suggested 



in Eq. (6.23) should be handled with care. 



Finally, in Fig. |6.10| on the following page we display the quantity 



because this ratio appears in the weighted asymmetries of Eq. ( |3.48| ). In Fig. |6.10| on the next page, 
the expression 

is also shown for comparison. Once again, there is a remarkable agreement between the two 
different expressions, confirming the quite similar behavior of and . 



6.2.3 Asymmetries in semi-inclusive DIS and e^e annihilation 

We turn now to estimates of possible observables containing the Collins function. We will take into 
consideration one-particle inclusive DIS, where the Collins function appears in connection with the 
transversity distribution of the nucleon, and e+e^ annihilation into two hadrons belonging to two 
different jets. 

In the first case, we show predictions for both transverse spin asymmetries defined in Eqs. (p.43p 



and ( |3.48D . Different calculations can be found in the literature, e.g., in Refs. ^|126| , p37| . To esti 



mate the magnitude of the asymmetries, we need inputs for the distribution functions, in particular 
for the transversity distribution. Several model calculations of this function are available at present 



referred to as "Collins ansatz". 
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Figure 6.10. Model result H^^^^ /D^ (solid line) and comparison with the product {{Kj) /2z^ml.){H^ /Df) 
( dashed line ). 



(see [ PT] ] for a comprehensive review). We refrain from considering many different examples and 
rather restrict the analysis to two limiting situations. In the first case we adopt the "nonrelativistic" 
assumption = g^, while in the second case we exhaust the upper bound on the transversity dis- 
tribution, i.e., = j{fi+gi) 0159| ]. We use the simple parametrization of and /j suggested in 
Ref. At the moment, more sophisticated parametrizations are available, taking scale evolution 
into account also. However, all these parametrizations are compatible with each other to the extent 
of our purpose here, which is to give an estimate of the asymmetries for a low scale. We focus on 
the production of n^, where the contribution of down quarks is negligible, not only because of the 
presence of unfavored fragmentation functions, but also because the transversity distribution for 
down quarks appears to be much smaller than for up quarks in model calculations. 

In Fig. |6.1 1| we present the azimuthal asymmetry defined in Eq. ( p.43D as a function of Xg, after 
integrating numerator and denominator over the variables y and Zj^, for the two cases described 
above. In Fig. |6.12| , we present the same asymmetry as a function of z^, after integrating over y 
and Xg. As already mentioned before, our prediction is supposed to be valid at a low energy scale of 
about 1 GeV^. Neglecting evolution effects, it could be utilized for comparison with experiments 
at a scale of a few GeV^. We assume the value of the transverse polarization to be = 0.75. 
In performing the integrations, we apply the kinematical cuts typical of the HERMES experiment, 
as described in [[T^]. Therefore, our prediction is particularly significant for HERMES, which 
is expected to be the first experiment to measure this asymmetry. In principle, the simultaneous 
study of the Xg and Z/^ dependence of the asymmetry yields separate information on the distribution 
and fragmentation parts and allows one to extract both up to a normalization factor [ |126| ]. Note, 
however, that this procedure relies on the assumption of up-quark dominance and is valid only if 
the Xg dependence of the asymmetry can be ascribed entirely to the distribution functions and the 
Zj^ dependence entirely to the fragmentation functions. Kinematical cuts could partially spoil this 
situation. We would like to stress that our calculation predicts an asymmetry up to the order of 
10%, which should be within experimental reach, and suggests the possibility of distinguishing 
between different assumptions on the transversity distribution. 

Using the same procedure as before, we have estimated the asymmetry defined in Eq. ( p.48[ ). 



6.2 Estimates and phenomenology 



107 





0.25 






Ed 


0.2 








0.15 








0.1 




0.05 




Figure 6.11. Azimuthal transverse spin asymmetry (sin0)j^j, as a function of x^. Solid line: assuming 



/ij = gi. Dashed line: assuming h^ — \{fi +gi)- The functions f^ and g^ are takenfrom /|^ 
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Figure 6.12. Azimuthal transverse spin asymmetry {sm<j))jjj as a function of Zf^. Solid line: assuming 



/ij = gy Dashed line: assuming h^ — ^(/j The functions f^ and g^ are takenfrom /|7( 
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Figure 6.13. Azimuthal spin asymmetry (iP^jj^l/m^j sin0)yj- as a function of Xg. Solid line: assuming 



/ij = gy Dashed line: assuming h^ — j{fi + gi)- The functions and g^ are taken from /|7( 
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Figure 6.14. Azimuthal spin asymmetry (jP^j^l/m^j sin^)jyj- as a function of Zf,. Solid line: assuming 
/ij — gy Dashed line: assuming /ij = ^(/j + gi). The functions /j and g^ are taken from ^ 



containing the weighting with l/m^j;. The results are shown in Fig. |6.13| as a function of and 
in Fig. |6.14| as a function of z^. The magnitude of this asymmetry is higher than in the unweighted 
case, which is partially due to the fact that the weighting enhances the asymmetry by about a factor 
of 2. 

In addition to appearing in semi-inclusive DIS in connection with the transversity distribution 
of the nucleon, the Collins function can be independently extracted from another process, i.e. 



electron-positron annihilation into two hadrons belonging to two back-to-back jets [pT|,p2|]. We 



restrict ourselves to the case of photon exchange only. In this process, one of the two hadrons (say 
hadron 2) defines the scattering plane together with the leptons and determines the direction with 
respect to which the azimuthal angles must be measured. The cross section is differential in five 
variables, e.g. Zi,Z2,y, |^/^^|, 0- The variables and are the longitudinal fractional momenta of 
the two hadrons. In the center of mass frame y — {1+ cos 0)/2, where 6 is the angle of hadron 
2 with respect to the momentum of the incoming leptons. The vector Pj^^ denotes the transverse 
component of the momentum of hadron 1 and ^ is its azimuthal angle with respect to the scattering 
plane. For a more detailed description of the kinematical variables we refer to [ pT| , |5^ . 
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Figure 6.15. Azimuthal asymmetry {P^j_ cos20)^^^_ /or e+e annihilation into two hadrons, integrated over 
the range 0.2 <Z2< 0.8 ( solid line), and over the range 0.5 < Z2 ^ 0.8 f dashed line). 



We define the azimuthal asymmetry 



(P^^COS20>^+^_(0,Z1,Z2) 







P cos 20 ^e+g- 






2d5a+ _ 



2 sin^ 



(6.26) 



where summations over quark flavors are understood. The weighting with a second power of 
\Pf^^ I in the numerator is necessary to obtain a deconvoluted expression. We prefer to use the same 
weighting in the denominator as well, to avoid a modification of the asymmetry just caused by the 
weighting. 

In Fig. |6.15| we present the estimate of the asymmetry defined above, entirely based on our 
model. The asymmetry has been integrated over and 9, leaving the dependence on alone. 
We have extended the integration interval all the way to [0, Tz], to obtain a conservative estimate. 
In fact, limiting the interval to [7r/4, 3;r/4] will enhance the asymmetry by a factor of 2, approxi- 
mately. Because the Collins function increases with increasing z, we also get a larger asymmetry 
by restricting the integration range for Zj- As an illustration of this feature, in Fig. |6.15| we present 
two results, obtained from two different integration ranges. Our prediction is supposed to be valid 
only at low energy scales and should be evolved for comparison with higher energy experiments. 
It is important to note that we estimate the asymmetry to be of the order of about 5%, and thus it 
should observable in experiments. 



6.3 Comparison with existing data 

The HERMES collaboration recently measured the single longitudinal spin asymmetry 

1 



UL' 



Q 



c^hj^{x) + C2h^{x)\H^^^/^\z) + other terms 



(6.27) 
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On purpose, we avoid entering the details of the formula. It is difficult to extract information on the 
Collins function and on the transversity distribution from this asymmetry [|6T|]. In a recent analysis 
of this asymmetry, Efremov et al. extracted a behavior H^/D^ z for z < 0.7, although some 
questionable assumptions were used to obtain this result. 

If we assume that the other terms in ( |6.27[ ) are small, then the dependence of the asymmetry 
should be almost entirely due to the Collins function. To compare this behavior with our model 
estimate, as a first step we parametrize our result for the ratio /D^ with a simple analytic 

form 

Hf^"'\z) 1 1 

— -— ^0.316z + 0.0345 0.00359- (6.28) 

£>i(z) l-z (l-z) 

Fig. |6.16| shows the result of our model together with this parametrization (note that we extended 
the plot to higher values of z compared to Fig. ^]8| on page |104| ). 

In Fig. |6.17| we compare our parametrization with HERMES data on (sin^y^)^;^ [ |TT1 ] and pre- 
liminary data on the same asymmetry from the CLAS collaboration at JLAB [^. Note that we 
arbitrarily normalized our curve to take into account the unknown distribution functions and pref- 
actors. The agreement of the Zj^ shape is remarkable. 



6.4 Summary 

We have estimated the Collins fragmentation function for pions at a low energy scale by means of 
the Manohar-Georgi model. This model contains three essential features of nonperturbative QCD: 
massive quark degrees of freedom, chiral symmetry and its spontaneous breaking (with pions as 
Goldstone bosons). Because of the chiral invariant interaction between pions and quarks, the frag- 
mentation process can be evaluated in a perturbative expansion. The constituent quark mass, the 
axial pion-quark coupling and the maximum virtuality /i^ of the fragmenting quark are free 
parameters of our approach. The quark mass and are constrained within reasonable limits. To 
ensure the convergence of the chiral perturbation expansion, jU^ cannot exceed a typical hadronic 
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Figure 6.17. Comparison between the results of our model and data from the HERMES and CLAS experi- 
ments. 



scale. We have mostly considered the value jU^ = 1 GeV^, which guarantees that the momenta of 
the particles produced in the fragmentation process stay well below the scale of chiral symmetry 
breaking, ^ 1 GeV. To determine the appropriate value of /i^, the average transverse momen- 
tum of a data set could be used. In any case, we observed that variations of the free parameters 
within reasonable limits have only a minor influence on the shape of the results, implying that our 
approach has a good predictive power for the z behavior of the various functions. 

We have found that the Manohar-Georgi model reproduces reasonably well the unpolarized 
pion fragmentation function and the average transverse momentum of a produced hadron as a 
function of z, supporting the idea of describing the fragmentation process by such a chiral invariant 
approach. 

Compared to the unpolarized fragmentation function, modeling the Collins function is con- 
siderably more difficult, mainly because of its chiral-odd and T-odd nature. In our approach, the 
helicity flip required to generate a chiral-odd object is caused by the mass of the constituent quark, 
while the T-odd behavior is produced via one-loop corrections. The Collins function exhibits a 
quite distinct behavior from the unpolarized fragmentation function. In particular, the ratio H^/D^ 
is strongly increasing with increasing z. 

On the basis of our results, we have calculated the transverse single-spin asymmetry in semi- 
inclusive DIS where the Collins function appears in combination with the transversity of the nu- 
cleon. This observable will be measured in the near future at HERMES and could also be inves- 
tigated at COMPASS, Jlab (upgraded) and EIC. For typical HERMES kinematics the asymmetry 
is of the order of 10%, giving support to the intention of extracting the nucleon transversity in 
this way. We believe that our estimate of the Collins function, despite its uncertainties, can be 
very useful for this extraction. Finally, we found an encouraging agreement between the shape 
of our estimate and the trend of the single spin asymmetry measured by the HERMES coUabora- 
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tion [jT^,[TT]], although it is not clear if this asymmetry is originated solely by the Collins function. 

More information on the Collins function from the experimental side is required. In this re- 
spect, the most promising experiment seems to be e+e^ annihilation into two hadrons, where 
appears squared in an azimuthal cos 20 asymmetry. According to our calculation, an asymmetry 
of the order of 5% can be expected, which should be measurable at high luminosity accelerators, 
such as BABAR and BELLE [ |146| ] . Dedicated measurements of the Collins function would be ex- 
tremely important for the extraction of the transversity distribution. Moreover, they could answer 
the question whether a chiral invariant Lagrangian can be used to model the Collins function. 



7 



Conclusions and outlook 



/ was born not knowing and 
have had only a little time to 
change that here and there. 

R. Feynman 



The main focus of this thesis was to present some ways of accessing the transversity distribution 
of quarks inside hadrons, hy This was the main motivation to perform a thorough analysis of 
one-particle and two-particle inclusive deep inelastic scattering at leading order in l/Q, with a 
particular attention to T-odd fragmentation functions. This thesis is certainly not the first work on 
this important subject, however it complements the existing literature in many respects. 



7.1 Conclusions 

Several observables containing the transversity distribution appear all along the thesis, but three of 
them are particularly promising and have been highlighted. The first method presented to access 
transversity is the measurement of the azimuthal single spin asymmetry described in Eq. ( [3.48D , or 



the similar asymmetry described in Eq. ( p.43[ ). In these asymmetries, the transversity distribution 



appears in combination with the Collins fragmentation function, which provides the chirality flip 
necessary to compensate the chiral-odd nature of /ij. 

The Collins function was introduced in Ref. ^ as early as 1993. However, it was looked at 
with some skepticism because it is a single-particle T-odd fragmentation function. It was conjec- 
tured to be small or to vanish altogether [ |117[ ] . In this thesis, we tackled the question by calculating 
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directly the Collins function in a consistent and time-reversal invariant field theoretical approach. 
Our model calculation dispels the doubts about the vanishing of the Collins function and, more 
in general, of T-odd fragmentation functions. Moreover, our calculation aspires to describe the 
Collins function in a qualitative and quantitative way. At the moment, there is no unquestionable 
experimental information on the Collins function. However, the single longitudinal spin asymme- 
try measured by the HERMES collaboration [jT^,[Tl|] can be originated by the Collins function. We 
tested our model on this experimental results and we found an encouraging agreement between the 
shape of our estimate and the trend of the data. 

It would be of remarkable importance to collect extra information on the Collins function, not 
only to allow the extraction of the transversity distribution, but also to open a new window on 
polarization in hadronization processes. New measurements could consider different asymmetries 
in deep inelastic scattering (HERMES, COMPASS, CLAS, EIC), pp scattering (RHIC) and e+e^ 
annihilation into two hadrons belonging to two jets (BABAR, BELLE). 

The second observable containing the transversity presented in this thesis is the asymmetry 
of Eq. (|4.42cD, involving a chiral-odd T-odd fragmentation function generated by the interference 



between the production of two hadrons in the s and p waves. This asymmetry was indicated for the 



first time by Jaffe, Jin and Tang in 1998 [ |1 17p . However, their analysis suffered two limitations: 
it did not present a complete treatment of two-hadron fragmentation functions and it relied just 
on one specific mechanism to give some phenomenological indications. An improved analysis of 
two-hadron fragmentation was presented by Bianconi, Boffi, Jakob and Radici p^, together with 
different models of the interference fragmentation function [ p3| , |150| ]. In this thesis, we comple- 



mented their formalism by performing a partial- wave decomposition of two-hadron fragmentation 
functions, with and without including partonic transverse momentum.^ We were able to separate 
the contribution of the s and p waves and we recovered the results of Jaffe, Jin and Tang in the 
transverse momentum integrated sp interference sector. 

Besides rediscussing the interference fragmentation function in a more exhaustive framework, 
another interesting asymmetry came out of the analysis, as shown in Eq. ( j4.43D . This asymme- 



try contains a chiral-odd T-odd function pertaining to the pure p-wave sector of two-hadron pro- 
duction. Such a contribution was already mentioned by Collins, Heppelmann and Ladinsky in 
1994 [[79]] and identified more formally by Ji in the same year [ |122| ]. However, our work is the first 
one to clearly distinguish the p-wave from the ^/^-interference part in the single transverse spin 
asymmetry. It is useful to stress once more that the p-wave contribution should not be overlooked. 
It vanishes when integrating over the full range of the polar angle 6, unlike the sp interference 
term. On the other hand, it might be larger than the latter - especially in the presence of a spin-one 
resonance -thus offering an excellent way to access the transversity distribution. 

We looked more carefully at p-wave two-hadron fragmentation functions from a different point 
of view, i.e. comparing them with spin-one fragmentation functions, introduced for the first time 



by Ji 01221 1. In the thesis, we carried out the analysis of spin-one functions in a more detailed 
way compared to the original work of Ji. We showed that they overlap almost exactly with pure 
p-wave two-hadron fragmentation functions, except for the lack of a dependence on the invariant 



The results of the partial-wave expansion are presented here for the first time and still await publication. 
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mass squared. We included partonic transverse momentum and we presented the complete cross 
section of spin-one leptoproduction at leading order inl/Q. As in the case of the Collins function, 
spin-one functions could provide us with fresh information on spin effects in hadronization. 

An important achievement of the thesis is the connection between the formalism of correlation 
functions and of helicity matrices. The two approaches were known to be equivalent at leading 
order in l/Q, but the connection had to be worked out more explicitly. We expressed in the he- 
licity matrix formalism all the partonic functions we considered - spin-half and spin-one distribu- 
tion functions, without and with transverse momentum, and one-hadron, two-hadron and spin-one 
fragmentation functions, without and with transverse momentum. A significant outcome of this 
examination was the derivation of positivity bounds for all of them. Positivity bounds can be valu- 
able tools, as they provide guidance to estimate unknown functions and they test the consistency 
of model calculations. 



7.2 Outlook 

The formalism of two-hadron and spin-one fragmentation functions has to be completed by analyz- 
ing all the azimuthal asymmetries containing the transversity distribution function, in case partonic 



transverse momentum is included. The only step in this direction was presented in Ref. |T]and was 
limited to the spin-one contributions. Subleading twist contributions have been neglected in the 
thesis and they deserve further survey. 

Another very useful extension would be to calculate the cross section of e~^e~ annihilation into 
two couples of hadrons belonging to two different jets, including in the analysis the distinction of 
the s- and p-wave contributions. A measurement of two-hadron fragmentation functions in this 
process would be extremely interesting first of all because - as already mentioned - they contain 
valuable information on the role of polarization in fragmentation processes, and secondly because 
- as in the case of the Collins function - an independent measurement of these functions would be 
very important for a clear extraction of the transversity distribution. 

For what concerns our model calculation of T-odd functions, an immediate application would 
be to use our results to estimate other asymmetries containing the Collins function. Furthermore, 
the shape of our results could be used to choose a specific analytical form to parametrize the Collins 
function and fit future experimental data. 

Another possible development is to use a similar approach to estimate T-odd interference frag- 
mentation functions and T-odd spin-one fragmentation functions. The analysis is complicated by 
the presence of a vector meson, which is beyond the reach of our model as formulated at present. 
In principle, it is possible to extend the model. In any case, even without any extension it would be 
already possible to explore the effect of one-loop corrections on "incoherent" two-pion production, 
and compare the relative importance of sp interference and pure p-wave contributions, even in the 
absence of a resonance (e.g. in tt^tt^ production). 

From a more formal point of view, the use of a chiral invariant approach for describing fragmen- 
tation functions at low energy scales should be investigated further. First of all, we should check 
if the chiral perturbation expansion is reliable. Secondly, we should try to describe fragmentation 
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functions entirely, and not their "valence" part only. Finally, it would be interesting to explore the 
connection between our one-loop corrections and the recent models of T-odd effects [[72|], which 
are also based on one-loop effects. 
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Summary 



This thesis discusses three different ways to observe the transverse spin of quarks inside the nucle- 
ons. This problem relates to one of the unanswered questions of present day physics: what is the 
internal structure of protons and neutrons? 

One of the ways to investigate such a question is to study the process of deep inelastic scatter- 
ing, in which a nucleon is bombarded with a focused beam of electrons (or other leptons) at very 
high energies. The electrons penetrate inside the nucleon and the way they are scattered yields 
information on its inner structure. A down-to-earth analogy would be shooting at a car with a ma- 
chine gun, penetrating its hood and recording the way bullets bounce on its internal components - 
engine, gear box and all. 

Deep inelastic scattering provides us with snapshots of the nucleon at high resolution, a few 
percents of the proton size. Extracting the relevant information on the structure of the nucleon 
is a complex task and requires a lot of theoretical effort. Due to the high energies and high mo- 
menta involved, we resort to the techniques of quantum field theory, which incorporates relativity 
and quantum mechanics. At low resolution - in low energy experiments - nucleons appear to be 
unbreakable particles with a certain mass, electric charge and spin. At high resolution - in deep in- 
elastic scattering - nucleons display an extremely intricate structure, involving a very high number 
of smaller particles, quarks and gluons. At present, we are not able to match these two pictures, 
in other words we are not able to explain how quarks and gluons interact to give origin to the 
characteristics of the nucleons. 

A key issue in understanding the nucleon is to explain the composition of its spin in terms of 
the underlying quark and gluon structure. So far, we acquired a good knowledge of the helicity 
distribution of quarks inside targets with a longitudinal spin. The exploration of the transverse 
spin distribution (or transversity) of quarks in targets with a transverse spin would give us a new 
perspective on the dynamics inside the nucleon. To a certain extent, the helicity and the transverse 
spin distributions represent a front and a side view of the nucleon spin. 

Observing the transversity of quarks is an arduous task. Several approaches have been sug- 
gested in theory, but none of them has been put into practice yet. In this thesis, three ways to probe 
the transversity are highlighted and analyzed. They are far from exhausting all possible methods, 
but they seem to be promising and could be adopted in experiments in the near future. All of them 
require that in the deep inelastic process we keep track not only of the scattered electron, but also 
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of one or two of the fragments coming from the disintegration of the nucleon. The detection of the 
fragments offers the opportunity to unravel more details of the structure of the target, but the price 
to pay is that we have to deal not only with the way the quarks are arranged in the nucleon, but 
also with the way they produce the final fragments, i.e. with fragmentation functions. 

From a more technical point of view, this thesis starts with reviewing the formalism to deal 
with deep inelastic scattering. The parton distribution functions are introduced, describing the 
way quarks and gluons are disposed inside the nucleon. Totally inclusive deep inelastic scattering, 
where only the scattered electron is detected, is discussed as the first and simplest way to observe 
distribution functions. To identify some way to measure the quark transversity, this thesis analyzes 
one-particle and two-particle inclusive deep inelastic scattering, where one and two of the outgoing 
hadrons are detected in coincidence with the electron. These processes require the introduction of 
parton fragmentation functions, describing the way a quark evolves into final state hadrons. This 
thesis shows that transversity can be measured in the above processes, in connection with three 
different fragmentation functions: the first one requires the observation of an unpolarized final 
state hadron with transverse momentum, the second requires the observation of the interference 
between the s- and p-wave production of two hadrons, the third requires the observation of pure 
p-wave two-hadron production, or equivalently of a spin-one resonance. All these fragmentation 
functions fall in the category of T-odd fragmentation functions: they require the presence of final 
state interactions, or else they are forbidden by time-reversal invariance. The last part of the thesis 
looks at the possibility of modeling this kind of fragmentation functions and investigates whether 
they can be large enough to allow the extraction of transversity. 



